Changeset 7656

Timestamp:

2017-02-08T14:55:40+01:00 (7 years ago)

Author:

isabella

Message:

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

File:

• branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90 (modified) (4 diffs)

branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

@@ -771 +771 @@
 
                         IF ( ln_incpc ) THEN
+                         IF (lwp) WRITE(numout,*) 'after end of IAU - IPC - saving s/tbias_asm_stscale'
                            DO jk = 1, jpkm1
                               tbias_asm_stscale_out(:,:,jk) =  tbias_asm_stscale(:,:,jk) * zfrac
@@ -785 +786 @@
 
             ELSE ! no assimilation increments, simply decay pressure correction (e.g for forecasts)
+                 ! this occurs at obsoper step as well ( ln_asmiau is false ) 
 
                DO jk = 1, jpkm1
@@ -806 +808 @@
 
                IF (ln_incpc) THEN
-                  
+                   IF (lwp) WRITE(numout,*) 'obsoper or forecast mode - IPC - computing tbias_i and sbias_i'
                    DO jk = 1, jpkm1
-                      tbias_i(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) )
-                      sbias_i(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) )
+                      tbias_i(:,:,jk) = tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) )
+                      sbias_i(:,:,jk) = sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) )
                    ENDDO
-
-                   IF ( kt == nn_bias_itwrt ) THEN
-                      DO jk = 1, jpkm1
-                         tbias_asm_stscale_out(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 )
-                         sbias_asm_stscale_out(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 )
-                      ENDDO
-                   ENDIF
-
                ENDIF
 
@@ -912 +906 @@
 
       IF ( kt == nn_bias_itwrt ) THEN
-
-
-         IF ( ln_bias_pc_app ) THEN
             tbias_p_out(:,:,:)   = tbias_p(:,:,:)
             sbias_p_out(:,:,:)   = sbias_p(:,:,:)
             IF (ln_incpc) THEN
                tbias_i_out(:,:,:)   = tbias_i(:,:,:)
-               sbias_i_out(:,:,:)   = sbias_i(:,:,:)
-                 IF (ln_incpc_only) THEN
-                  IF(lwp) WRITE(numout,*) ' t/p bias_p_out and tbias_asm_out are set to zero if incpc only'
-                   tbias_p_out(:,:,:)   = 0.0
-                   sbias_p_out(:,:,:)   = 0.0
-                   tbias_asm_out(:,:,:) = 0.0
-                   sbias_asm_out(:,:,:) = 0.0
-                 ENDIF
+               sbias_i_out(:,:,:)   = sbias_i(:,:,:)            
             ENDIF
-
-         ELSE
-
-            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'
-            tbias_p_out(:,:,:)   = 0.0
-            sbias_p_out(:,:,:)   = 0.0
-            tbias_asm_out(:,:,:) = 0.0
-            sbias_asm_out(:,:,:) = 0.0
-            IF (ln_incpc) THEN
-            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'
-               tbias_i_out(:,:,:)   = 0.0
-               sbias_i_out(:,:,:)   = 0.0
-               tbias_asm_stscale_out(:,:,:) = 0.0
-               sbias_asm_stscale_out(:,:,:) = 0.0
-            ENDIF
-         ENDIF
       ENDIF