Changeset 7656


Ignore:
Timestamp:
2017-02-08T14:55:40+01:00 (4 years ago)
Author:
isabella
Message:

added other output messages and changed one formula to sum up from tbias_i .

File:
1 edited

Legend:

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

    r7653 r7656  
    771771 
    772772                        IF ( ln_incpc ) THEN 
     773                         IF (lwp) WRITE(numout,*) 'after end of IAU - IPC - saving s/tbias_asm_stscale' 
    773774                           DO jk = 1, jpkm1 
    774775                              tbias_asm_stscale_out(:,:,jk) =  tbias_asm_stscale(:,:,jk) * zfrac 
     
    785786 
    786787            ELSE ! no assimilation increments, simply decay pressure correction (e.g for forecasts) 
     788                 ! this occurs at obsoper step as well ( ln_asmiau is false )  
    787789 
    788790               DO jk = 1, jpkm1 
     
    806808 
    807809               IF (ln_incpc) THEN 
    808                    
     810                   IF (lwp) WRITE(numout,*) 'obsoper or forecast mode - IPC - computing tbias_i and sbias_i' 
    809811                   DO jk = 1, jpkm1 
    810                       tbias_i(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    811                       sbias_i(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
     812                      tbias_i(:,:,jk) = tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
     813                      sbias_i(:,:,jk) = sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    812814                   ENDDO 
    813  
    814                    IF ( kt == nn_bias_itwrt ) THEN 
    815                       DO jk = 1, jpkm1 
    816                          tbias_asm_stscale_out(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 ) 
    817                          sbias_asm_stscale_out(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 ) 
    818                       ENDDO 
    819                    ENDIF 
    820  
    821815               ENDIF 
    822816 
     
    912906 
    913907      IF ( kt == nn_bias_itwrt ) THEN 
    914  
    915  
    916          IF ( ln_bias_pc_app ) THEN 
    917908            tbias_p_out(:,:,:)   = tbias_p(:,:,:) 
    918909            sbias_p_out(:,:,:)   = sbias_p(:,:,:) 
    919910            IF (ln_incpc) THEN 
    920911               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 
     912               sbias_i_out(:,:,:)   = sbias_i(:,:,:)             
    929913            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 
    945          ENDIF 
    946914      ENDIF 
    947915 
Note: See TracChangeset for help on using the changeset viewer.