Changeset 7653 for branches

Timestamp:

2017-02-07T16:37:29+01:00 (7 years ago)

Author:

isabella

Message:

adding changes to reset variables following options

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

@@ -404 +404 @@
            tbias_asm(:,:,:) = 0.0_wp
            sbias_asm(:,:,:) = 0.0_wp
+           tbias_asm_out(:,:,:) = 0.0_wp
+           sbias_asm_out(:,:,:) = 0.0_wp
          ENDIF
 
@@ -409 +411 @@
            tbias_asm_stscale(:,:,:) = 0.0_wp
            sbias_asm_stscale(:,:,:) = 0.0_wp
+           tbias_asm_stscale_out(:,:,:) = 0.0_wp
+           sbias_asm_stscale_out(:,:,:) = 0.0_wp
          ENDIF
 
@@ -645 +649 @@
                      &               ( t_asm_mem**(real(kt)/ztsday) * sbias_asm(:,:,jk)  +                     &
                      &               ztfrac * t_asm_upd * s_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof1(:,:)
-                     tbias_p(:,:,jk) = tbias_p(:,:,jk) +                        &
+
+                    tbias_p(:,:,jk) = tbias_p(:,:,jk) +                        &
                      &               ( t_asm_mem**(real(kt)/ztsday) * tbias_asm(:,:,jk)  +                     &
                      &               ztfrac * t_asm_upd * t_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof2(:,:)
@@ -907 +912 @@
 
       IF ( kt == nn_bias_itwrt ) 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_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
+            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