I used climate data (monthly mean value)for NO3, O2, and so on (COADS 05 and WOA2009). I used modis chlorophyll monthly mean value in 2016 as the initial data and forcing condition. phytoplankton was set by chlorophyll.
now, the model blow up.
i was wondering if the forcing condition (U,V, pair, tair) are not match with NH4 (constants for all grids)? should they match with the changes of the chlorophyll and phytoplankton?
what should i do when the model blow up? should i try to track the fortran code and find which step is something wrong?
thank you so much for your help in advance, best wishes
Code: Select all
326 0001-01-03 00:54:00.00 1.805362E-02 1.009324E+02 1.009504E+02 4.729840E+11 (004,028,18) 3.956102E-01 1.365641E-01 1.130226E+01 4.059272E+00 327 0001-01-03 01:03:00.00 1.707953E-02 1.009146E+02 1.009317E+02 4.728986E+11 (004,027,19) 3.526031E-01 1.868201E-01 1.182530E+01 5.002697E+00 328 0001-01-03 01:12:00.00 1.635206E-02 1.008860E+02 1.009024E+02 4.727603E+11 (006,025,18) 2.752542E-01 4.995576E-01 1.422893E+01 7.306116E+00 329 0001-01-03 01:21:00.00 1.656034E-02 1.008469E+02 1.008635E+02 4.725701E+11 (006,024,19) 2.383016E-01 1.426992E+00 2.045112E+01 1.410777E+01 330 0001-01-03 01:30:00.00 2.925909E-02 1.007976E+02 1.008268E+02 4.723289E+11 (006,023,20) 1.342662E-01 5.816475E+00 4.073108E+01 4.860259E+01 Blowing-up: Saving latest model state into RESTART file
Code: Select all
GET_NGFLD - ammonium southern boundary condition, 0001-01-16 00:00:00.00 (Grid= 01, Rec=0000001, Index=1, File: roms_bry.nc) (Tmin= 15.0000 Tmax= 345.0000) t = 15.0000 (Min = 5.00000000E+00 Max = 5.00000000E+00) GET_NGFLD - chlorophyll southern boundary condition, 0001-01-16 00:00:00.00 (Grid= 01, Rec=0000001, Index=1, File: roms_bry.nc) (Tmin= 15.0000 Tmax= 345.0000) t = 15.0000 (Min = 2.69459943E-02 Max = 1.90680868E+00)