Hi, everyone, I run ROMS and got the empty results (the NaN was shown in the first step and was attached below). I used Woapisces (have been provided by Olivier Aumont. These data have been computed by a coupled NEMO-PISCES global model) to compute the initial and boundary conditions of NO3 and O2.Then, I used WOA2009 ( World Ocean Atlas directory 2009 )to compute the initial and boundary conditions of temperature and salty. Then, I used coads dataset to compute forcing and boundary condition of Wind stress, Heat fluxes w3, Fresh water fluxes (evaporation - precipitation), Sea surface temperature and heat flux sensitivity to the sea surface temperature (dQdSST),Sea surface salinity,Short wave radiation, pair air pressure. the chlorophyll initial conditions was computed by using modis data in 2016,and NH4, LdetritusN, and et al were set as a constant with the time and loction in the boundary files. finally, i got a empty result and showd as below. Is there any way I can figure out what is the problem? the grid files? the water depth problem? or the parameter in the forcing file , boundary file ,and initial file are not match? or show i check some scripts to find out it?Thank you so much for your help in advance. I have long struggled with the problem TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed

for Sea surface temperature and salinity in forcing files were computed by using the coads dataset (coads05) ,but the salinity and temperature in the boundary files were computed by using woapisces dataset. Would that be a problem?

in that all the winds are negative. I would ask the model to write out surface fluxes and stresses to see if those look reasonable. I would try a much shorter timestep to see if you can get it to run long enough to output some useful fields so you can figure out where in the domain it is going bad.

I would not worry about starting from rest. The model spins up to match the density fields pretty quickly.

You had said you were applying the stresses from COADS, but yet the model is reading winds, air pressure and the rest. This looks odd:in that all the winds are negative. I would ask the model to write out surface fluxes and stresses to see if those look reasonable. I would try a much shorter timestep to see if you can get it to run long enough to output some useful fields so you can figure out where in the domain it is going bad.

thx,I really appreciated the suggestions you gave. it is important clue for me who only have biological background.I changed the timestep.

Code:

100 ntimes Number of timesteps for 3-D equations.90.000 dt Timestep size (s) for 3-D equations.5 ndtfast Number of timesteps for 2-D equations between

(100 steps were finished and I did not copy them all)the parameter only at first time step existed in the avg.nc file. I attached the nc file here.I will try to turnoff the biological mode part and check the forcing, grid files again.

i did my best to figure out why latent,sensible,and shlux is so wield in history file "roms_his.nc" (results were attached below). however, I have no clue. "lwrad" seems ok in history file(roms_his.nc).I go back to check the forcing file (roms_frc.nc) which I used "coads" and "NCEP/NCAR Reanalysis Monthly Means" data (attached below).

it is only one effective value in my region (the region is too small), and this value was obtained from coads data. acutally, i set the value of TIC,alkalinity by myself. NO3 which from coads data was multiplied by a different factor to get TIC,alkalinity, and so on.

Who is online

Users browsing this forum: No registered users and 1 guest

You cannot post new topics in this forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum