Changeset 8036


Ignore:
Timestamp:
2017-05-17T19:39:22+02:00 (3 years ago)
Author:
jwhile
Message:

Further merge of r5518_pcbias

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  
    13081308   rn_htrmax         =  200.0   ! max. depth of transition range 
    13091309/ 
     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  
    4444   USE wrk_nemo        ! Memory Allocation 
    4545   USE timing          ! Timing 
     46   USE biaspar         ! bias correction variables 
    4647 
    4748   IMPLICIT NONE 
     
    8485      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    8586      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 
    8690      !!---------------------------------------------------------------------- 
    8791      ! 
     
    9296         ztrdu(:,:,:) = ua(:,:,:) 
    9397         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(:,:) 
    94107      ENDIF 
    95108      ! 
     
    112125      IF(ln_ctl)   CALL prt_ctl( tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask,   & 
    113126         &                       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 
    114133      ! 
    115134      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  
    466466                            CALL dia_hsb_init   ! heat content, salt content and volume budgets 
    467467                            CALL     trd_init   ! Mixed-layer/Vorticity/Integral constraints trends 
     468                            CALL     bias_init  ! Pressure correction bias 
    468469      IF( lk_diaobs     ) THEN                  ! Observation & model comparison 
    469470                            CALL dia_obs_init            ! Initialize observational data 
  • branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7950 r8036  
    103103      IF( ln_crs     )       CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" )   ! tell iom we are at time step kstp 
    104104 
     105      IF( ln_bias )          CALL bias_opn( kstp ) 
     106 
    105107      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    106108      ! Update data, open boundaries, surface boundary condition (including sea-ice) 
     
    269271      IF( lk_trabbl      )   CALL tra_bbl    ( kstp )       ! advective (and/or diffusive) bottom boundary layer scheme 
    270272      IF( ln_tradmp      )   CALL tra_dmp    ( kstp )       ! internal damping trends 
     273      IF( ln_bias        )   CALL tra_bias   ( kstp ) 
    271274      IF( lk_bdy         )   CALL bdy_tra_dmp( kstp )       ! bdy damping trends 
    272275                             CALL tra_adv    ( kstp )       ! horizontal & vertical advection 
     
    293296               &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    294297               &                                    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 ) 
    295299      ELSE                                                  ! centered hpg  (eos then time stepping) 
    296300         IF ( .NOT. lk_dynspg_ts ) THEN                     ! eos already called in time-split case 
     
    307311         IF( ln_zdfnpc   )   CALL tra_npc( kstp )                ! update after fields by non-penetrative convection 
    308312                             CALL tra_nxt( kstp )                ! tracer fields at next time step 
     313         IF( ln_bias )       CALL dyn_bias( kstp ) 
    309314      ENDIF 
    310315 
     
    381386      ENDIF 
    382387 
     388 
     389      IF( lrst_bias )          CALL bias_wrt     ( kstp ) 
     390 
    383391      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    384392      ! Coupled mode 
  • branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r6491 r8036  
    100100 
    101101   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) 
    103105   USE asminc           ! assimilation increments      (tra_asm_inc routine) 
    104106   !                                                   (dyn_asm_inc routine) 
Note: See TracChangeset for help on using the changeset viewer.