nmathewa / MOM6dev

MOM6 development
0 stars 0 forks source link

the simulation shows results on land and not on sea #32

Closed Parag0206 closed 2 years ago

Parag0206 commented 2 years ago

The simulation produces results(all fields like SST,SSS,etc.) on land and not on the sea.

Screenshot from 2022-04-04 11-55-26

This maybe due to improper mask

Parag0206 commented 2 years ago

Yes, the msks were the problem, the land and ocean masks both were swapped. But, after correcting, still the error persists

Parag0206 commented 2 years ago

The topo file generation method is not the reason of the problem(whether topog file is generated from ETOPO or GEBCO does not matter).

Parag0206 commented 2 years ago

further the ocean domain shows that all values are missing

nmathewa commented 2 years ago
nmathewa commented 2 years ago

bath = etopo.depth.values

n_bath = -1*(np.where(bath > 0,0,bath))

Take a look, values greater than 0 will be 0 and less than zero will be multplied by -1

nmathewa commented 2 years ago
Parag0206 commented 2 years ago

As mentioned above, the problem is that the ocean values for bathymetry are below 0(sign convention). Convert them to positive as mentioned here

bath = etopo.depth.values

n_bath = -1*(np.where(bath > 0,0,bath))

Take a look, values greater than 0 will be 0 and less than zero will be multplied by -1

or use cdo expr to solve this