Changeset 8036
- Timestamp:
- 2017-05-17T19:39:22+02:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM
- Files:
-
- 5 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/CONFIG/SHARED/namelist_ref
r7941 r8036 1308 1308 rn_htrmax = 200.0 ! max. depth of transition range 1309 1309 / 1310 !----------------------------------------------------------------------- 1311 &nambias ! Bias pressure correctiom 1312 !----------------------------------------------------------------------- 1313 ln_bias = .false. 1314 ln_bias_asm = .false. 1315 ln_bias_rlx = .false. 1316 ln_bias_ofl = .false. 1317 ln_bias_ts_app = .false. 1318 ln_bias_pc_app = .false. 1319 fb_t_asm = 0.0 1320 fb_t_rlx = 0.0 1321 fb_t_ofl = 1.0 1322 fb_p_asm = 1.0 1323 fb_p_rlx = 1.0 1324 fb_p_ofl = 0.0 1325 eft_rlx = 365.0 1326 eft_asm = 365.0 1327 t_rlx_upd = 0.1 1328 t_asm_upd = 0.1 1329 nn_lat_ramp = 0 1330 bias_time_unit_asm = 86400.0 1331 bias_time_unit_rlx = 1.0 1332 bias_time_unit_ofl = 1.0 1333 cn_bias_tot = "bias_tot.nc" 1334 cn_bias_asm = "bias_asm.nc" 1335 cn_dir = './' 1336 ln_bsyncro = .FALSE. 1337 fctamp = 1. 1338 rn_maxlat_bias = 23.0 1339 rn_minlat_bias = 10.0 1340 nn_bias_itwrt = 15 1341 ln_itdecay = .FALSE. 1342 / -
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r6486 r8036 44 44 USE wrk_nemo ! Memory Allocation 45 45 USE timing ! Timing 46 USE biaspar ! bias correction variables 46 47 47 48 IMPLICIT NONE … … 84 85 INTEGER, INTENT(in) :: kt ! ocean time-step index 85 86 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdu, ztrdv 87 REAL(wp), DIMENSION(jpi,jpj,jpk) :: z_rhd_st ! tmp density storage for pressure corr 88 REAL(wp), DIMENSION(jpi,jpj) :: z_gru_st ! tmp ua trends storage for pressure corr 89 REAL(wp), DIMENSION(jpi,jpj) :: z_grv_st ! tmp va trends storage for pressure corr 86 90 !!---------------------------------------------------------------------- 87 91 ! … … 92 96 ztrdu(:,:,:) = ua(:,:,:) 93 97 ztrdv(:,:,:) = va(:,:,:) 98 ENDIF 99 ! 100 IF ( ln_bias .AND. ln_bias_pc_app ) THEN 101 z_rhd_st(:,:,:) = rhd(:,:,:) ! store orig density 102 rhd(:,:,:) = rhd_pc(:,:,:) ! use pressure corrected density 103 z_gru_st(:,:) = gru(:,:) 104 gru(:,:) = gru_pc(:,:) 105 z_grv_st(:,:) = grv(:,:) 106 grv(:,:) = grv_pc(:,:) 94 107 ENDIF 95 108 ! … … 112 125 IF(ln_ctl) CALL prt_ctl( tab3d_1=ua, clinfo1=' hpg - Ua: ', mask1=umask, & 113 126 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 127 ! 128 IF ( ln_bias .AND. ln_bias_pc_app ) THEN 129 rhd(:,:,:) = z_rhd_st(:,:,:) ! restore original density 130 gru(:,:) = z_gru_st(:,:) 131 grv(:,:) = z_grv_st(:,:) 132 ENDIF 114 133 ! 115 134 IF( nn_timing == 1 ) CALL timing_stop('dyn_hpg') -
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r7563 r8036 466 466 CALL dia_hsb_init ! heat content, salt content and volume budgets 467 467 CALL trd_init ! Mixed-layer/Vorticity/Integral constraints trends 468 CALL bias_init ! Pressure correction bias 468 469 IF( lk_diaobs ) THEN ! Observation & model comparison 469 470 CALL dia_obs_init ! Initialize observational data -
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/step.F90
r7950 r8036 103 103 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell iom we are at time step kstp 104 104 105 IF( ln_bias ) CALL bias_opn( kstp ) 106 105 107 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 106 108 ! Update data, open boundaries, surface boundary condition (including sea-ice) … … 269 271 IF( lk_trabbl ) CALL tra_bbl ( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 270 272 IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends 273 IF( ln_bias ) CALL tra_bias ( kstp ) 271 274 IF( lk_bdy ) CALL bdy_tra_dmp( kstp ) ! bdy damping trends 272 275 CALL tra_adv ( kstp ) ! horizontal & vertical advection … … 293 296 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 294 297 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 298 IF( ln_bias ) CALL dyn_bias( kstp ) 295 299 ELSE ! centered hpg (eos then time stepping) 296 300 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case … … 307 311 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection 308 312 CALL tra_nxt( kstp ) ! tracer fields at next time step 313 IF( ln_bias ) CALL dyn_bias( kstp ) 309 314 ENDIF 310 315 … … 381 386 ENDIF 382 387 388 389 IF( lrst_bias ) CALL bias_wrt ( kstp ) 390 383 391 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 384 392 ! Coupled mode -
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r6491 r8036 100 100 101 101 USE crsfld ! Standard output on coarse grid (crs_fld routine) 102 102 USE biaspar ! bias param 103 USE bias ! bias routines (tra_bias routine 104 ! dyn_bias routine) 103 105 USE asminc ! assimilation increments (tra_asm_inc routine) 104 106 ! (dyn_asm_inc routine)
Note: See TracChangeset
for help on using the changeset viewer.