Changeset 7089


Ignore:
Timestamp:
2016-10-25T16:39:26+02:00 (4 years ago)
Author:
isabella
Message:

simplifiyng the code

Location:
branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6328 r7089  
    13311331   nn_bias_itwrt  = 15 
    13321332   ln_itdecay     = .FALSE. 
    1333 / 
     1333   ln_incpc       = .FALSE. 
     1334/ 
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

    r7082 r7089  
    147147      & sbias_rlx_out, &   !: Output salinity bias field 
    148148      & tbias_p_out,   &   !: Output temperature bias field for P correction 
    149       & sbias_p_out        !: Output salinity bias field for P correction 
     149      & sbias_p_out,   &   !: Output salinity bias field for P correction 
    150150 
    151151   INTEGER, PRIVATE :: nn_lat_ramp     ! choice of latitude dependent ramp 
     
    155155   LOGICAL, PRIVATE :: ln_itdecay      ! evolve bias correction at every time step.   
    156156   LOGICAL, PRIVATE :: ln_incpc        ! incremental pressure correction 
    157    LOGICAL, PRIVATE :: lat_ramp_incpc  ! ramp for incremental pressure correction 
    158157 
    159158   REAL(wp), PRIVATE, DIMENSION(:,:), ALLOCATABLE :: fbcoef 
     
    217216         & cn_bias_tot, cn_bias_asm, cn_dir, sn_tbias_ofl, sn_sbias_ofl,  & 
    218217         & ln_bsyncro, fctamp, rn_maxlat_bias, rn_minlat_bias,            & 
    219          & nn_bias_itwrt, ln_itdecay, ln_incpc, lat_ramp_incpc 
     218         & nn_bias_itwrt, ln_itdecay, ln_incpc 
    220219  
    221220 
     
    277276    WRITE(numout,*) '     evolve pcbias at each timestep ln_itdecay     = ',ln_itdecay 
    278277    WRITE(numout,*) '     incremental press. correction  ln_incpc       = ',ln_incpc 
    279          WRITE(numout,*) '     lat ramp for inc.press.corr    lat_ramp_incpc = ',lat_ramp_incpc 
    280278         WRITE(numout,*) '  Parameters for parition of bias term ' 
    281279         WRITE(numout,*) '                                    fb_t_rlx       = ',fb_t_rlx 
     
    309307            &      tbias_p(jpi,jpj,jpk), & 
    310308            &      sbias_p(jpi,jpj,jpk), & 
    311             &      tbias_i(jpi,jpj,jpk), & 
    312             &      sbias_i(jpi,jpj,jpk), & 
    313309            &      rhd_pc(jpi,jpj,jpk) , & 
    314310            &      gru_pc(jpi,jpj)     , & 
     
    520516 
    521517 
    522       IF ( lat_ramp_incpc) THEN 
     518      IF ( ln_incpc) THEN 
    523519       ! not sure if this should be a special case of nn_lat_ramp == 2 
    524520         minlat_bias = 3.0_wp 
     
    582578      ztscale = 0.1_wp 
    583579      zwgt    = 1.0_wp 
    584       tbias_i(:,:,:) = 0.0_wp 
    585       sbias_i(:,:,:) = 0.0_wp 
    586580 
    587581      IF ( ln_bias_asm ) THEN 
     
    635629 
    636630                     DO jk = 1, jpkm1 
    637  
    638                         tbias_i(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     631                        tbias_p(:,:,jk) = tbias_p(:,:,jk) + ( t_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    639632                        &                - ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    640                         sbias_i(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     633                        sbias_p(:,:,jk) = sbias_p(:,:,jk) + ( s_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    641634                        &                - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    642  
    643635                     ENDDO 
    644636 
     
    688680                    
    689681                   DO jk = 1, jpkm1 
    690                       tbias_i(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     682                      tbias_p(:,:,jk) = tbias_p(:,:,jk) + ( t_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    691683                      &                - ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    692                       sbias_i(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     684                      sbias_p(:,:,jk) = sbias_p(:,:,jk) + ( s_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    693685                      &                - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    694                     ENDDO 
     686                   ENDDO 
    695687 
    696688                 ENDIF 
     
    864856         tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_p(:,:,:) 
    865857         tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_p(:,:,:) 
    866  
    867          IF (ln_incpc) THEN 
    868             tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) -tbias_i(:,:,:) - tbias_p(:,:,:) 
    869             tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) -sbias_i(:,:,:) - sbias_p(:,:,:) 
    870          ENDIF 
    871  
    872858      ELSE 
    873859         tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_p(:,:,:) 
    874860         tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) 
    875  
    876          IF (ln_incpc) THEN 
    877             tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_i(:,:,:) - tbias_p(:,:,:) 
    878             tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_i(:,:,:) - sbias_p(:,:,:) 
    879          ENDIF 
    880  
    881861      ENDIF 
    882862 
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/biaspar.F90

    r6328 r7089  
    2323   LOGICAL, PUBLIC :: ln_bias_pc_app = .FALSE. !: estimate bias from assim incr 
    2424   LOGICAL, PUBLIC :: lrst_bias      = .FALSE. !: estimate bias from assim incr 
     25   LOGICAL, PUBLIC :: ln_incpc       = .FALSE. !: incr pressure correction 
    2526 
    2627   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: & 
Note: See TracChangeset for help on using the changeset viewer.