Changeset 12609
- Timestamp:
- 2020-03-26T11:53:23+01:00 (5 years ago)
- Location:
- branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 3 edited
-
ASM/asminc.F90 (modified) (5 diffs)
-
OBS/diaobs.F90 (modified) (1 diff)
-
OBS/obs_oper.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r12594 r12609 60 60 #endif 61 61 LOGICAL, PUBLIC :: ln_bkgwri = .FALSE. !: No output of the background state fields 62 LOGICAL, PUBLIC :: ln_avgbkg = .FALSE. !: No output of the mean background state fields 62 63 LOGICAL, PUBLIC :: ln_asmiau = .FALSE. !: No applying forcing with an assimilation increment 63 64 LOGICAL, PUBLIC :: ln_asmdin = .FALSE. !: No direct initialization … … 133 134 REAL(wp), POINTER, DIMENSION(:,:) :: hdiv ! 2D workspace 134 135 !! 135 NAMELIST/nam_asminc/ ln_bkgwri, ln_ balwri,&136 NAMELIST/nam_asminc/ ln_bkgwri, ln_avgbkg, ln_balwri, & 136 137 & ln_trainc, ln_dyninc, ln_sshinc, & 137 & ln_phytobal, ln_slchltotinc, ln_slchldiainc, &138 & ln_slchlnoninc, ln_schltotinc, ln_slphytotinc, &139 & ln_slphydiainc, ln_slphynoninc, ln_spco2inc, &140 & ln_sfco2inc, ln_plchltotinc, ln_pchltotinc, &141 & ln_pno3inc, ln_psi4inc, ln_pdicinc, ln_palkinc, &142 & ln_pphinc, ln_po2inc, &143 138 & ln_asmdin, ln_asmiau, & 144 139 & nitbkg, nitdin, nitiaustr, nitiaufin, niaufn, & 145 & ln_salfix, salfixmin, nn_divdmp, & 146 & ln_seaiceinc, ln_sitinc, ln_temnofreeze, & 147 & mld_choice_bgc, rn_maxchlinc 140 & ln_salfix, salfixmin, nn_divdmp 148 141 !!---------------------------------------------------------------------- 149 142 … … 171 164 WRITE(numout,*) ' Namelist nam_asminc : set assimilation increment parameters' 172 165 WRITE(numout,*) ' Logical switch for writing out background state ln_bkgwri = ', ln_bkgwri 166 WRITE(numout,*) ' Logical switch for writing mean background state ln_avgbkg = ', ln_avgbkg 173 167 WRITE(numout,*) ' Logical switch for writing out balancing increments ln_balwri = ', ln_balwri 174 168 WRITE(numout,*) ' Logical switch for applying tracer increments ln_trainc = ', ln_trainc 175 169 WRITE(numout,*) ' Logical switch for applying velocity increments ln_dyninc = ', ln_dyninc 176 170 WRITE(numout,*) ' Logical switch for applying SSH increments ln_sshinc = ', ln_sshinc 177 WRITE(numout,*) ' Logical switch for Direct Initialization (DI) ln_asmdin = ', ln_asmdin178 171 WRITE(numout,*) ' Logical switch for applying SIC increments ln_seaiceinc = ', ln_seaiceinc 179 172 WRITE(numout,*) ' Logical switch for applying SIT increments ln_sitinc = ', ln_sitinc 180 WRITE(numout,*) ' Logical switch for phytoplankton balancing ln_phytobal = ', ln_phytobal 181 WRITE(numout,*) ' Logical switch for applying slchltot increments ln_slchltotinc = ', ln_slchltotinc 182 WRITE(numout,*) ' Logical switch for applying slchldia increments ln_slchldiainc = ', ln_slchldiainc 183 WRITE(numout,*) ' Logical switch for applying slchlnon increments ln_slchlnoninc = ', ln_slchlnoninc 184 WRITE(numout,*) ' Logical switch for applying schltot increments ln_schltotinc = ', ln_schltotinc 185 WRITE(numout,*) ' Logical switch for applying slphytot increments ln_slphytotinc = ', ln_slphytotinc 186 WRITE(numout,*) ' Logical switch for applying slphydia increments ln_slphydiainc = ', ln_slphydiainc 187 WRITE(numout,*) ' Logical switch for applying slphynon increments ln_slphynoninc = ', ln_slphynoninc 188 WRITE(numout,*) ' Logical switch for applying spco2 increments ln_spco2inc = ', ln_spco2inc 189 WRITE(numout,*) ' Logical switch for applying sfco2 increments ln_sfco2inc = ', ln_sfco2inc 190 WRITE(numout,*) ' Logical switch for applying plchltot increments ln_plchltotinc = ', ln_plchltotinc 191 WRITE(numout,*) ' Logical switch for applying pchltot increments ln_pchltotinc = ', ln_pchltotinc 192 WRITE(numout,*) ' Logical switch for applying pno3 increments ln_pno3inc = ', ln_pno3inc 193 WRITE(numout,*) ' Logical switch for applying psi4 increments ln_psi4inc = ', ln_psi4inc 194 WRITE(numout,*) ' Logical switch for applying pdic increments ln_pdicinc = ', ln_pdicinc 195 WRITE(numout,*) ' Logical switch for applying palk increments ln_palkinc = ', ln_palkinc 196 WRITE(numout,*) ' Logical switch for applying pph increments ln_pphinc = ', ln_pphinc 197 WRITE(numout,*) ' Logical switch for applying po2 increments ln_po2inc = ', ln_po2inc 173 WRITE(numout,*) ' Logical switch for Direct Initialization (DI) ln_asmdin = ', ln_asmdin 198 174 WRITE(numout,*) ' Logical switch for Incremental Analysis Updating (IAU) ln_asmiau = ', ln_asmiau 199 175 WRITE(numout,*) ' Timestep of background in [0,nitend-nit000-1] nitbkg = ', nitbkg … … 204 180 WRITE(numout,*) ' Logical switch for ensuring that the sa > salfixmin ln_salfix = ', ln_salfix 205 181 WRITE(numout,*) ' Minimum salinity after applying the increments salfixmin = ', salfixmin 206 WRITE(numout,*) ' Choice of MLD for BGC assimilation mld_choice_bgc = ', mld_choice_bgc207 WRITE(numout,*) ' Maximum absolute chlorophyll increment (<=0 = off) rn_maxchlinc = ', rn_maxchlinc208 182 ENDIF 209 183 … … 365 339 !-------------------------------------------------------------------- 366 340 367 ALLOCATE( t_bkginc(jpi,jpj,jpk) ) 368 ALLOCATE( s_bkginc(jpi,jpj,jpk) ) 369 ALLOCATE( u_bkginc(jpi,jpj,jpk) ) 370 ALLOCATE( v_bkginc(jpi,jpj,jpk) ) 371 ALLOCATE( ssh_bkginc(jpi,jpj) ) 372 ALLOCATE( seaice_bkginc(jpi,jpj)) 373 ALLOCATE( sit_bkginc(jpi,jpj) ) 341 IF ( ln_trainc ) THEN 342 ALLOCATE( t_bkginc(jpi,jpj,jpk) ) 343 ALLOCATE( s_bkginc(jpi,jpj,jpk) ) 344 t_bkginc(:,:,:) = 0.0 345 s_bkginc(:,:,:) = 0.0 346 ENDIF 347 IF ( ln_dyninc ) THEN 348 ALLOCATE( u_bkginc(jpi,jpj,jpk) ) 349 ALLOCATE( v_bkginc(jpi,jpj,jpk) ) 350 u_bkginc(:,:,:) = 0.0 351 v_bkginc(:,:,:) = 0.0 352 ENDIF 353 IF ( ln_sshinc ) THEN 354 ALLOCATE( ssh_bkginc(jpi,jpj) ) 355 ssh_bkginc(:,:) = 0.0 356 ENDIF 357 IF ( ln_seaiceinc ) THEN 358 ALLOCATE( seaice_bkginc(jpi,jpj)) 359 seaice_bkginc(:,:) = 0.0 360 ENDIF 361 IF ( ln_sitinc ) THEN 362 ALLOCATE( sit_bkginc(jpi,jpj)) 363 sit_bkginc(:,:) = 0.0 364 ENDIF 374 365 #if defined key_asminc 375 366 ALLOCATE( ssh_iau(jpi,jpj) ) 376 #endif377 t_bkginc(:,:,:) = 0.0378 s_bkginc(:,:,:) = 0.0379 u_bkginc(:,:,:) = 0.0380 v_bkginc(:,:,:) = 0.0381 ssh_bkginc(:,:) = 0.0382 seaice_bkginc(:,:) = 0.0383 sit_bkginc(:,:) = 0.0384 #if defined key_asminc385 367 ssh_iau(:,:) = 0.0 386 368 #endif -
branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r12380 r12609 968 968 ENDIF 969 969 970 #if defined key_cice 970 971 IF ( TRIM(cobstypessurf(jtype)) == 'sit' ) THEN 971 972 CALL obs_rea_snowdepth( surfdataqc(jtype), n2dintsurf(jtype), thick_s(:,:) ) 972 973 ENDIF 974 #endif 973 975 974 976 IF ( TRIM(cobstypessurf(jtype)) == 'sst' .AND. ln_sstbias ) THEN -
branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90
r12380 r12609 950 950 IF ( surfdataqc%lclim ) surfdataqc%rclm(jobs,1) = zclm(1) 951 951 952 #if defined key_cice 952 953 IF ( TRIM(surfdataqc%cvars(1)) == 'FBD' ) THEN 953 954 ! Convert radar freeboard to true freeboard (add 1/4 snow depth; 1/4 based on ratio of speed of light in vacuum compared to snow (3.0e8 vs 2.4e8 m/s)) … … 971 972 surfdataqc%nqc(jobs) = 4 972 973 ENDIF 974 #endif defined key_cice 973 975 974 976 END DO
Note: See TracChangeset
for help on using the changeset viewer.
