Changeset 6278 for branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO
- Timestamp:
- 2016-02-03T11:05:58+01:00 (8 years ago)
- Location:
- branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r6277 r6278 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 ! for pressure correction 88 REAL(wp), DIMENSION(jpi,jpj) :: z_gru_st ! for pressure correction 89 REAL(wp), DIMENSION(jpi,jpj) :: z_grv_st ! for pressure correction 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_pcbias/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90
r6277 r6278 121 121 INTEGER :: numout = 6 !: logical unit for output print; Set to stdout to ensure any early 122 122 ! output can be collected; do not change 123 INTEGER :: numnam = -1 !: logical unit for namelist 123 124 INTEGER :: numnam_ref = -1 !: logical unit for reference namelist 124 125 INTEGER :: numnam_cfg = -1 !: logical unit for configuration specific namelist -
branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r6277 r6278 452 452 IF( lk_diadct ) CALL dia_dct_init ! Sections tranports 453 453 CALL dia_hsb_init ! heat content, salt content and volume budgets 454 CALL trd_init ! Mixed-layer/Vorticity/Integral constraints trends 454 CALL trd_init ! Mixed-layer/Vorticity/Integral constraints trend 455 IF (lk_bias ) CALL bias_init ! Pressure correction bias 455 456 IF( lk_diaobs ) THEN ! Observation & model comparison 456 457 CALL dia_obs_init ! Initialize observational data -
branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO/OPA_SRC/step.F90
r6277 r6278 97 97 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell iom we are at time step kstp 98 98 99 IF( lk_bias ) THEN 100 IF( ln_bias ) CALL bias_opn( kstp ) 101 ENDIF 102 99 103 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 100 104 ! Update data, open boundaries, surface boundary condition (including sea-ice) … … 255 259 IF( lk_trabbl ) CALL tra_bbl ( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 256 260 IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends 261 IF( lk_bias ) THEN 262 IF( ln_bias ) CALL tra_bias ( kstp ) 263 ENDIF 257 264 IF( lk_bdy ) CALL bdy_tra_dmp( kstp ) ! bdy damping trends 258 265 CALL tra_adv ( kstp ) ! horizontal & vertical advection … … 279 286 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 280 287 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 288 IF( lk_bias ) THEN 289 IF( ln_bias ) CALL dyn_bias( kstp ) 290 ENDIF 281 291 ELSE ! centered hpg (eos then time stepping) 282 292 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case … … 293 303 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection 294 304 CALL tra_nxt( kstp ) ! tracer fields at next time step 305 IF( lk_bias ) THEN 306 IF( ln_bias ) CALL dyn_bias( kstp ) 307 ENDIF 295 308 ENDIF 296 309 … … 353 366 ENDIF 354 367 IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file 368 IF( lk_bias ) THEN 369 IF( lrst_bias ) CALL bias_wrt ( kstp ) 370 ENDIF 355 371 356 372 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> -
branches/UKMO/dev_r5518_pcbias/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r6277 r6278 99 99 100 100 USE crsfld ! Standard output on coarse grid (crs_fld routine) 101 101 USE biaspar ! bias param 102 USE bias ! bias routines (tra_bias routine 103 ! dyn_bias routine) 102 104 USE asminc ! assimilation increments (tra_asm_inc routine) 103 105 ! (dyn_asm_inc routine)
Note: See TracChangeset
for help on using the changeset viewer.