- Timestamp:
- 2017-08-09T14:56:31+02:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM
- Files:
-
- 1 added
- 5 edited
- 1 copied
-
CONFIG/SHARED/namelist_ref (modified) (2 diffs)
-
NEMO/OPA_SRC/ASM/asmbal.F90 (copied) (copied from branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/ASM/asmbal.F90) (2 diffs)
-
NEMO/OPA_SRC/ASM/asminc.F90 (modified) (16 diffs)
-
NEMO/OPA_SRC/ASM/asmlogchlbal_hadocc.F90 (added)
-
NEMO/OPA_SRC/ASM/asmpar.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/nemogcm.F90 (modified) (3 diffs)
-
NEMO/OPA_SRC/step.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/CONFIG/SHARED/namelist_ref
r8401 r8428 1293 1293 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1294 1294 ln_seaiceinc = .false. ! Logical switch for applying sea ice increments 1295 ln_logchlinc = .false. ! Logical switch for applying logchl increments 1296 ln_logchlbal = .false. ! Logical switch for applying logchl multivariate balancing 1295 1297 ln_temnofreeze = .false. ! Logical to not add increments if temperature would fall below freezing 1296 1298 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] … … 1302 1304 salfixmin = -9999 ! Minimum salinity after applying the increments 1303 1305 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1306 mld_choice_bgc = 5 ! MLD criterion to use for biogeochemistry assimilation 1307 rn_maxchlinc = -999.0 ! maximum absolute non-log chlorophyll increment from logchl assimilation 1308 ! <= 0 implies no maximum applied (switch turned off) 1309 ! > 0 implies maximum absolute chl increment capped at this value 1304 1310 / 1305 1311 !----------------------------------------------------------------------- -
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbal.F90
r7726 r8428 16 16 USE asminc ! Main assimilation increments module 17 17 USE asmpar ! Parameters for the assimilation interface 18 #if defined key_fabm 19 USE par_fabm 20 #elif defined key_medusa && defined key_foam_medusa 18 #if defined key_medusa && defined key_foam_medusa 21 19 USE par_medusa 22 20 #elif defined key_hadocc … … 75 73 CALL iom_rstput( kt, kt, inum, 'rdastp' , zdate ) 76 74 77 IF ( ln_logchltotinc .OR. ln_logchlpftinc ) THEN 78 #if defined key_fabm 79 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl1', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 80 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl2', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 81 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl3', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 82 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl4', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 83 #elif defined key_medusa && defined key_foam_medusa 75 IF ( ln_logchlinc ) THEN 76 #if defined key_medusa && defined key_foam_medusa 84 77 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chn', logchl_balinc(:,:,:,jpchn) ) 85 78 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chd', logchl_balinc(:,:,:,jpchd) ) 86 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phn', logchl_balinc(:,:,:,jpphn) ) 87 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phd', logchl_balinc(:,:,:,jpphd) ) 88 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_pds', logchl_balinc(:,:,:,jppds) ) 89 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zmi', logchl_balinc(:,:,:,jpzmi) ) 90 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zme', logchl_balinc(:,:,:,jpzme) ) 91 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_din', logchl_balinc(:,:,:,jpdin) ) 92 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_sil', logchl_balinc(:,:,:,jpsil) ) 93 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_fer', logchl_balinc(:,:,:,jpfer) ) 94 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jpdet) ) 95 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dtc', logchl_balinc(:,:,:,jpdtc) ) 96 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jpdic) ) 97 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jpalk) ) 98 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_oxy', logchl_balinc(:,:,:,jpoxy) ) 79 IF ( ln_logchlbal ) THEN 80 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phn', logchl_balinc(:,:,:,jpphn) ) 81 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phd', logchl_balinc(:,:,:,jpphd) ) 82 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_pds', logchl_balinc(:,:,:,jppds) ) 83 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zmi', logchl_balinc(:,:,:,jpzmi) ) 84 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zme', logchl_balinc(:,:,:,jpzme) ) 85 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_din', logchl_balinc(:,:,:,jpdin) ) 86 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_sil', logchl_balinc(:,:,:,jpsil) ) 87 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_fer', logchl_balinc(:,:,:,jpfer) ) 88 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jpdet) ) 89 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dtc', logchl_balinc(:,:,:,jpdtc) ) 90 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jpdic) ) 91 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jpalk) ) 92 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_oxy', logchl_balinc(:,:,:,jpoxy) ) 93 ENDIF 99 94 #elif defined key_hadocc 100 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_nut', logchl_balinc(:,:,:,jp_had_nut) )101 95 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phy', logchl_balinc(:,:,:,jp_had_phy) ) 102 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zoo', logchl_balinc(:,:,:,jp_had_zoo) ) 103 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jp_had_det) ) 104 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jp_had_dic) ) 105 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jp_had_alk) ) 96 IF ( ln_logchlbal ) THEN 97 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_nut', logchl_balinc(:,:,:,jp_had_nut) ) 98 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zoo', logchl_balinc(:,:,:,jp_had_zoo) ) 99 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jp_had_pdn) ) 100 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jp_had_dic) ) 101 CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jp_had_alk) ) 102 ENDIF 106 103 #endif 107 104 ENDIF -
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r8400 r8428 22 22 !! ssh_asm_inc : Apply the SSH increment 23 23 !! seaice_asm_inc : Apply the seaice increment 24 !! logchl_asm_inc : Apply the logchl increment 24 25 !!---------------------------------------------------------------------- 25 26 USE wrk_nemo ! Memory Allocation … … 44 45 #endif 45 46 USE sbc_oce ! Surface boundary condition variables. 47 USE zdfmxl, ONLY : & 48 ! & hmld_tref, & 49 #if defined key_karaml 50 & hmld_kara, & 51 & ln_kara, & 52 #endif 53 & hmld, & 54 & hmlp, & 55 & hmlpt 56 #if defined key_top 57 USE trc, ONLY: & 58 & trn, & 59 & trb 60 USE par_trc, ONLY: & 61 & jptra 62 #endif 63 #if defined key_medusa && defined key_foam_medusa 64 USE asmlogchlbal_medusa, ONLY: & 65 & asm_logchl_bal_medusa 66 USE par_medusa 67 #elif defined key_hadocc 68 USE asmlogchlbal_hadocc, ONLY: & 69 & asm_logchl_bal_hadocc 70 USE par_hadocc 71 #endif 46 72 47 73 IMPLICIT NONE … … 54 80 PUBLIC ssh_asm_inc !: Apply the SSH increment 55 81 PUBLIC seaice_asm_inc !: Apply the seaice increment 82 PUBLIC logchl_asm_inc !: Apply the logchl increment 56 83 57 84 #if defined key_asminc … … 61 88 #endif 62 89 LOGICAL, PUBLIC :: ln_bkgwri = .FALSE. !: No output of the background state fields 90 LOGICAL, PUBLIC :: ln_balwri = .FALSE. !: No output of the assimilation balancing increments 63 91 LOGICAL, PUBLIC :: ln_asmiau = .FALSE. !: No applying forcing with an assimilation increment 64 92 LOGICAL, PUBLIC :: ln_asmdin = .FALSE. !: No direct initialization … … 67 95 LOGICAL, PUBLIC :: ln_sshinc = .FALSE. !: No sea surface height assimilation increment 68 96 LOGICAL, PUBLIC :: ln_seaiceinc !: No sea ice concentration increment 97 LOGICAL, PUBLIC :: ln_logchlinc = .FALSE. !: No log10(chlorophyll) increment 98 LOGICAL, PUBLIC :: ln_logchlbal = .FALSE. !: Don't apply multivariate log10(chlorophyll) balancing 69 99 LOGICAL, PUBLIC :: ln_salfix = .FALSE. !: Apply minimum salinity check 70 100 LOGICAL, PUBLIC :: ln_temnofreeze = .FALSE. !: Don't allow the temperature to drop below freezing … … 91 121 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ssh_bkg, ssh_bkginc ! Background sea surface height and its increment 92 122 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: seaice_bkginc ! Increment to the background sea ice conc 123 124 REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE :: logchl_bkginc !: Increment to background logchl 125 #if defined key_top 126 REAL(wp), PUBLIC, DIMENSION(:,:,:,:), ALLOCATABLE :: logchl_balinc !: Increment to BGC variables from logchl assim 127 #endif 128 REAL(wp) :: rn_maxchlinc = -999.0 !: maximum absolute non-log chlorophyll increment from logchl assimilation 129 !: <= 0 implies no maximum applied (switch turned off) 130 !: > 0 implies maximum absolute chl increment capped at this value 131 132 INTEGER :: mld_choice_bgc = 5 !: choice of mld criteria to use for biogeochemistry assimilation 133 !: 1) hmld - Turbocline/mixing depth [W points] 134 !: 2) hmlp - Density criterion (0.01 kg/m^3 change from 10m) [W points] 135 !: 3) hmld_kara - Kara MLD [Interpolated] 136 !: 4) hmld_tref - Temperature criterion (0.2 K change from surface) [T points] 137 !: 5) hmlpt - Density criterion (0.01 kg/m^3 change from 10m) [T points] 93 138 94 139 !! * Substitutions … … 133 178 REAL(wp), POINTER, DIMENSION(:,:) :: hdiv ! 2D workspace 134 179 !! 135 NAMELIST/nam_asminc/ ln_bkgwri, 180 NAMELIST/nam_asminc/ ln_bkgwri, ln_balwri, & 136 181 & ln_trainc, ln_dyninc, ln_sshinc, & 182 & ln_logchlinc, ln_logchlbal, & 137 183 & ln_asmdin, ln_asmiau, & 138 184 & nitbkg, nitdin, nitiaustr, nitiaufin, niaufn, & 139 185 & ln_salfix, salfixmin, nn_divdmp, & 140 & ln_seaiceinc, ln_temnofreeze 186 & ln_seaiceinc, ln_temnofreeze, & 187 & mld_choice_bgc, rn_maxchlinc 141 188 !!---------------------------------------------------------------------- 142 189 … … 163 210 WRITE(numout,*) ' Namelist namasm : set assimilation increment parameters' 164 211 WRITE(numout,*) ' Logical switch for writing out background state ln_bkgwri = ', ln_bkgwri 212 WRITE(numout,*) ' Logical switch for writing out balancing increments ln_balwri = ', ln_balwri 165 213 WRITE(numout,*) ' Logical switch for applying tracer increments ln_trainc = ', ln_trainc 166 214 WRITE(numout,*) ' Logical switch for applying velocity increments ln_dyninc = ', ln_dyninc … … 168 216 WRITE(numout,*) ' Logical switch for Direct Initialization (DI) ln_asmdin = ', ln_asmdin 169 217 WRITE(numout,*) ' Logical switch for applying sea ice increments ln_seaiceinc = ', ln_seaiceinc 218 WRITE(numout,*) ' Logical switch for applying logchl increments ln_logchlinc = ', ln_logchlinc 219 WRITE(numout,*) ' Logical switch for logchl multivariate balancing ln_logchlbal = ', ln_logchlbal 170 220 WRITE(numout,*) ' Logical switch for Incremental Analysis Updating (IAU) ln_asmiau = ', ln_asmiau 171 221 WRITE(numout,*) ' Timestep of background in [0,nitend-nit000-1] nitbkg = ', nitbkg … … 176 226 WRITE(numout,*) ' Logical switch for ensuring that the sa > salfixmin ln_salfix = ', ln_salfix 177 227 WRITE(numout,*) ' Minimum salinity after applying the increments salfixmin = ', salfixmin 228 WRITE(numout,*) ' Choice of MLD for BGC assimilation mld_choice_bgc = ', mld_choice_bgc 229 WRITE(numout,*) ' Maximum absolute chlorophyll increment (<=0 = off) rn_maxchlinc = ', rn_maxchlinc 178 230 ENDIF 179 231 … … 223 275 224 276 IF ( ( ( .NOT. ln_asmdin ).AND.( .NOT. ln_asmiau ) ) & 225 .AND.( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ) .OR. ( ln_seaiceinc) )) & 226 & CALL ctl_stop( ' One or more of ln_trainc, ln_dyninc, ln_sshinc and ln_seaiceinc is set to .true.', & 277 & .AND.( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ).OR.( ln_seaiceinc ).OR. & 278 & ( ln_logchlinc ) )) & 279 & CALL ctl_stop( ' One or more of ln_trainc, ln_dyninc, ln_sshinc, ln_seaiceinc,', & 280 & ' and ln_logchlinc is set to .true.', & 227 281 & ' but ln_asmdin and ln_asmiau are both set to .false. :', & 228 282 & ' Inconsistent options') … … 233 287 234 288 IF ( ( .NOT. ln_trainc ).AND.( .NOT. ln_dyninc ).AND.( .NOT. ln_sshinc ).AND.( .NOT. ln_seaiceinc ) & 235 & ) & 236 & CALL ctl_warn( ' ln_trainc, ln_dyninc, ln_sshinc and ln_seaiceinc are set to .false. :', & 289 & .AND.( .NOT. ln_logchlinc ) ) & 290 & CALL ctl_warn( ' ln_trainc, ln_dyninc, ln_sshinc, ln_seaiceinc,', & 291 & ' and ln_logchlinc are set to .false. :', & 237 292 & ' The assimilation increments are not applied') 238 293 … … 254 309 & ' Background time step for Direct Initialization is outside', & 255 310 & ' the cycle interval') 311 312 IF ( ( ln_balwri ).AND.( .NOT. ( ln_logchlinc ) ) ) THEN 313 CALL ctl_warn( ' Balancing increments are only calculated for logchl', & 314 & ' Not assimilating logchl, so ln_balwri will be set to .false.') 315 ln_balwri = .FALSE. 316 ENDIF 317 318 IF ( ( ln_logchlbal ).AND.( .NOT. ( ln_logchlinc ) ) ) THEN 319 CALL ctl_warn( ' Cannot calculate logchl balancing increments if not assimilating logchl', & 320 & ' Not assimilating logchl, so ln_logchlbal will be set to .false.') 321 ln_logchlbal = .FALSE. 322 ENDIF 256 323 257 324 IF ( nstop > 0 ) RETURN ! if there are any errors then go no further … … 350 417 ssh_iau(:,:) = 0.0 351 418 #endif 352 IF ( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ).OR.( ln_seaiceinc ) ) THEN 419 IF ( ln_logchlinc ) THEN 420 ALLOCATE( logchl_bkginc(jpi,jpj) ) 421 logchl_bkginc(:,:) = 0.0 422 #if defined key_top 423 ALLOCATE( logchl_balinc(jpi,jpj,jpk,jptra) ) 424 logchl_balinc(:,:,:,:) = 0.0 425 #endif 426 ENDIF 427 IF ( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ).OR.( ln_seaiceinc ) & 428 & .OR.( ln_logchlinc ) ) THEN 353 429 354 430 !-------------------------------------------------------------------- … … 422 498 ! to allow for differences in masks 423 499 WHERE( ABS( seaice_bkginc(:,:) ) > 1.0e+10 ) seaice_bkginc(:,:) = 0.0 500 ENDIF 501 502 IF ( ln_logchlinc ) THEN 503 CALL iom_get( inum, jpdom_autoglo, 'bckinlogchl', logchl_bkginc(:,:), 1 ) 504 ! Apply the masks 505 logchl_bkginc(:,:) = logchl_bkginc(:,:) * tmask(:,:,1) 506 ! Set missing increments to 0.0 rather than 1e+20 507 ! to allow for differences in masks 508 WHERE( ABS( logchl_bkginc(:,:) ) > 1.0e+10 ) logchl_bkginc(:,:) = 0.0 424 509 ENDIF 425 510 … … 1155 1240 1156 1241 END SUBROUTINE seaice_asm_inc 1242 1243 SUBROUTINE logchl_asm_inc( kt ) 1244 !!---------------------------------------------------------------------- 1245 !! *** ROUTINE logchl_asm_inc *** 1246 !! 1247 !! ** Purpose : Apply the chlorophyll assimilation increments. 1248 !! 1249 !! ** Method : Calculate increments to state variables using nitrogen 1250 !! balancing. 1251 !! Direct initialization or Incremental Analysis Updating. 1252 !! 1253 !! ** Action : 1254 !!---------------------------------------------------------------------- 1255 INTEGER, INTENT(IN) :: kt ! Current time step 1256 ! 1257 INTEGER :: jk ! Loop counter 1258 INTEGER :: it ! Index 1259 REAL(wp) :: zincwgt ! IAU weight for current time step 1260 REAL(wp) :: zincper ! IAU interval in seconds 1261 !!---------------------------------------------------------------------- 1262 1263 IF ( kt <= nit000 ) THEN 1264 1265 zincper = (nitiaufin_r - nitiaustr_r + 1) * rdt 1266 1267 #if defined key_medusa && defined key_foam_medusa 1268 !CALL asm_logchl_bal_medusa() 1269 CALL ctl_stop( 'Attempting to assimilate logchl into MEDUSA, ', & 1270 & 'but not fully implemented yet' ) 1271 #elif defined key_hadocc 1272 CALL asm_logchl_bal_hadocc( logchl_bkginc, zincper, mld_choice_bgc, & 1273 & rn_maxchlinc, ln_logchlbal, logchl_balinc ) 1274 #else 1275 CALL ctl_stop( 'Attempting to assimilate logchl, ', & 1276 & 'but not defined a biogeochemical model' ) 1277 #endif 1278 1279 ENDIF 1280 1281 IF ( ln_asmiau ) THEN 1282 1283 !-------------------------------------------------------------------- 1284 ! Incremental Analysis Updating 1285 !-------------------------------------------------------------------- 1286 1287 IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 1288 1289 it = kt - nit000 + 1 1290 zincwgt = wgtiau(it) ! IAU weight for the current time step 1291 ! note this is not a tendency so should not be divided by rdt 1292 1293 IF(lwp) THEN 1294 WRITE(numout,*) 1295 WRITE(numout,*) 'logchl_asm_inc : logchl IAU at time step = ', & 1296 & kt,' with IAU weight = ', wgtiau(it) 1297 WRITE(numout,*) '~~~~~~~~~~~~' 1298 ENDIF 1299 1300 ! Update the biogeochemical variables 1301 ! Add directly to trn and trb, rather than to tra, as not a tendency 1302 #if defined key_medusa && defined key_foam_medusa 1303 DO jk = 1, jpkm1 1304 trn(:,:,jk,jp_msa0:jp_msa1) = trn(:,:,jk,jp_msa0:jp_msa1) + & 1305 & logchl_balinc(:,:,jk,jp_msa0:jp_msa1) * zincwgt 1306 trb(:,:,jk,jp_msa0:jp_msa1) = trb(:,:,jk,jp_msa0:jp_msa1) + & 1307 & logchl_balinc(:,:,jk,jp_msa0:jp_msa1) * zincwgt 1308 END DO 1309 #elif defined key_hadocc 1310 DO jk = 1, jpkm1 1311 trn(:,:,jk,jp_had0:jp_had1) = trn(:,:,jk,jp_had0:jp_had1) + & 1312 & logchl_balinc(:,:,jk,jp_had0:jp_had1) * zincwgt 1313 trb(:,:,jk,jp_had0:jp_had1) = trb(:,:,jk,jp_had0:jp_had1) + & 1314 & logchl_balinc(:,:,jk,jp_had0:jp_had1) * zincwgt 1315 END DO 1316 #endif 1317 1318 ! Do not deallocate arrays - needed by asm_bal_wri 1319 ! which is called at end of model run 1320 ENDIF 1321 1322 ELSEIF ( ln_asmdin ) THEN 1323 1324 !-------------------------------------------------------------------- 1325 ! Direct Initialization 1326 !-------------------------------------------------------------------- 1327 1328 IF ( kt == nitdin_r ) THEN 1329 1330 neuler = 0 ! Force Euler forward step 1331 1332 #if defined key_medusa && defined key_foam_medusa 1333 ! Initialize the now fields with the background + increment 1334 ! Background currently is what the model is initialised with 1335 CALL ctl_warn( ' Doing direct initialisation of MEDUSA with chlorophyll assimilation', & 1336 & ' Background state is taken from model rather than background file' ) 1337 trn(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) + & 1338 & logchl_balinc(:,:,:,jp_msa0:jp_msa1) 1339 trb(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) 1340 #elif defined key_hadocc 1341 ! Initialize the now fields with the background + increment 1342 ! Background currently is what the model is initialised with 1343 CALL ctl_warn( ' Doing direct initialisation of HadOCC with chlorophyll assimilation', & 1344 & ' Background state is taken from model rather than background file' ) 1345 trn(:,:,:,jp_had0:jp_had1) = trn(:,:,:,jp_had0:jp_had1) + & 1346 & logchl_balinc(:,:,:,jp_had0:jp_had1) 1347 trb(:,:,:,jp_had0:jp_had1) = trn(:,:,:,jp_had0:jp_had1) 1348 #endif 1349 1350 ! Do not deallocate arrays - needed by asm_bal_wri 1351 ! which is called at end of model run 1352 ENDIF 1353 ! 1354 ENDIF 1355 ! 1356 END SUBROUTINE logchl_asm_inc 1157 1357 1158 1358 !!====================================================================== -
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmpar.F90
r6486 r8428 20 20 & c_asmtrj = 'assim_trj', & !: Filename for storing the 21 21 !: reference trajectory 22 & c_asminc = 'assim_background_increments' 22 & c_asminc = 'assim_background_increments', & !: Filename for storing the 23 23 !: increments to the background 24 24 !: state 25 & c_asmbal = 'assim.balincs' !: Filename for storing the 26 !: balancing increments calculated 27 !: for biogeochemistry 25 28 26 29 INTEGER, PUBLIC :: nitbkg_r !: Background time step referenced to nit000 -
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r8400 r8428 61 61 USE asminc ! assimilation increments 62 62 USE asmbkg ! writing out state trajectory 63 USE asmbal ! writing out assimilation balancing increments 63 64 USE diaptr ! poleward transports (dia_ptr_init routine) 64 65 USE diadct ! sections transports (dia_dct_init routine) … … 160 161 IF( ln_dyninc ) CALL dyn_asm_inc( nit000 - 1 ) ! Dynamics 161 162 IF( ln_sshinc ) CALL ssh_asm_inc( nit000 - 1 ) ! SSH 163 IF( ln_logchlinc ) CALL logchl_asm_inc( nit000 - 1 ) 162 164 ENDIF 163 165 ENDIF … … 191 193 192 194 IF( lk_diaobs ) CALL dia_obs_wri 195 ! 196 IF( ( lk_asminc ).AND.( ln_balwri ) ) CALL asm_bal_wri( nitend ) ! Output balancing increments 193 197 ! 194 198 IF( ln_icebergs ) CALL icb_end( nitend ) -
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/step.F90
r8400 r8428 253 253 ! Passive Tracer Model 254 254 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 255 IF( lk_asminc .AND. ln_asmiau .AND. ln_logchlinc ) & 256 & CALL logchl_asm_inc( kstp ) ! logchl assimilation 255 257 CALL trc_stp( kstp ) ! time-stepping 256 258 #endif
Note: See TracChangeset
for help on using the changeset viewer.
