Changeset 7653


Ignore:
Timestamp:
2017-02-07T16:37:29+01:00 (4 years ago)
Author:
isabella
Message:

adding changes to reset variables following options

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

    r7652 r7653  
    404404           tbias_asm(:,:,:) = 0.0_wp 
    405405           sbias_asm(:,:,:) = 0.0_wp 
     406           tbias_asm_out(:,:,:) = 0.0_wp 
     407           sbias_asm_out(:,:,:) = 0.0_wp 
    406408         ENDIF 
    407409 
     
    409411           tbias_asm_stscale(:,:,:) = 0.0_wp 
    410412           sbias_asm_stscale(:,:,:) = 0.0_wp 
     413           tbias_asm_stscale_out(:,:,:) = 0.0_wp 
     414           sbias_asm_stscale_out(:,:,:) = 0.0_wp 
    411415         ENDIF 
    412416 
     
    645649                     &               ( t_asm_mem**(real(kt)/ztsday) * sbias_asm(:,:,jk)  +                     & 
    646650                     &               ztfrac * t_asm_upd * s_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof1(:,:) 
    647                      tbias_p(:,:,jk) = tbias_p(:,:,jk) +                        & 
     651 
     652                    tbias_p(:,:,jk) = tbias_p(:,:,jk) +                        & 
    648653                     &               ( t_asm_mem**(real(kt)/ztsday) * tbias_asm(:,:,jk)  +                     & 
    649654                     &               ztfrac * t_asm_upd * t_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof2(:,:) 
     
    907912 
    908913      IF ( kt == nn_bias_itwrt ) THEN 
    909          tbias_p_out(:,:,:)   = tbias_p(:,:,:) 
    910          sbias_p_out(:,:,:)   = sbias_p(:,:,:) 
    911          IF (ln_incpc) THEN 
    912             tbias_i_out(:,:,:)   = tbias_i(:,:,:) 
    913             sbias_i_out(:,:,:)   = sbias_i(:,:,:) 
     914 
     915 
     916         IF ( ln_bias_pc_app ) THEN 
     917            tbias_p_out(:,:,:)   = tbias_p(:,:,:) 
     918            sbias_p_out(:,:,:)   = sbias_p(:,:,:) 
     919            IF (ln_incpc) THEN 
     920               tbias_i_out(:,:,:)   = tbias_i(:,:,:) 
     921               sbias_i_out(:,:,:)   = sbias_i(:,:,:) 
     922                 IF (ln_incpc_only) THEN 
     923                  IF(lwp) WRITE(numout,*) ' t/p bias_p_out and tbias_asm_out are set to zero if incpc only' 
     924                   tbias_p_out(:,:,:)   = 0.0 
     925                   sbias_p_out(:,:,:)   = 0.0 
     926                   tbias_asm_out(:,:,:) = 0.0 
     927                   sbias_asm_out(:,:,:) = 0.0 
     928                 ENDIF 
     929            ENDIF 
     930 
     931         ELSE 
     932 
     933            IF(lwp) WRITE(numout,*) ' ln_bias_pc_app is =',ln_bias_pc_app, 'so t/p bias_p_out and tbias_asm_out are set to zero' 
     934            tbias_p_out(:,:,:)   = 0.0 
     935            sbias_p_out(:,:,:)   = 0.0 
     936            tbias_asm_out(:,:,:) = 0.0 
     937            sbias_asm_out(:,:,:) = 0.0 
     938            IF (ln_incpc) THEN 
     939            IF(lwp) WRITE(numout,*) ' ln_bias_pc_app is =',ln_bias_pc_app,' ln_incpc =',ln_incpc,'so t/p bias_i_out and tbias_asm_stscale_out are set to zero' 
     940               tbias_i_out(:,:,:)   = 0.0 
     941               sbias_i_out(:,:,:)   = 0.0 
     942               tbias_asm_stscale_out(:,:,:) = 0.0 
     943               sbias_asm_stscale_out(:,:,:) = 0.0 
     944            ENDIF 
    914945         ENDIF 
    915946      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.