Hello,
I am writing in connection with ROMS_3.7_NPZD_IRON model set-up.
The model is running but it is not reading dissolved iron eastern boundary condition.
I prescribed open boundary for NO3 and iron and closed boundary for the rest variables in npzd_iron.in file.
Initially, the model was not running as it was asking for phytFe_east boundary.
Hereafter I prescribed a small value (0.001 millimole/m3) for the phytFe_east boundary.
Then, the model started running but it is not reading dissolved iron eastern boundary condition instead it is reading phytoplankton, iron concentration eastern boundary condition.
**********
GET_NGFLD - temperature eastern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 3.01271591E+01)
GET_NGFLD - salinity eastern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 3.54105110E+01)
GET_NGFLD - nitrate eastern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 3.33553314E+01)
GET_NGFLD - phytoplankton, iron concentration eastern bt = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 1.00000000E-03 Max = 1.00000000E-03)
GET_NGFLD - temperature southern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 2.35944996E+01)
GET_NGFLD - salinity southern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 3.59412613E+01)
GET_NGFLD - nitrate southern boundary condition, t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 1.92311306E+01)
GET_NGFLD - dissolved iron southern boundary condition,t = 3986 19:30:00
(Rec=0000012, Index=2, File: roms_bry_1x12_ocean_time.nc)
(Tmin= 3652.0000 Tmax= 3986.8125)
(Min = 0.00000000E+00 Max = 7.90963888E-01)
**********
I am not sure but it seems that the model code is not reading iron eastern boundary condition properly.
Best,
Kunal