Changeset 7552 for branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

Timestamp:

2017-01-13T09:14:12+01:00 (7 years ago)

Author:

isabella

Message:

added one switch ln_incpc_only to enable selection of IPC scheme only

File:

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

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

@@ -226 +226 @@
          & cn_bias_tot, cn_bias_asm, cn_dir, sn_tbias_ofl, sn_sbias_ofl,  &
          & ln_bsyncro, fctamp, rn_maxlat_bias, rn_minlat_bias,            &
-         & nn_bias_itwrt, ln_itdecay, ln_incpc, ztscale, zwgt
+         & nn_bias_itwrt, ln_itdecay, ln_incpc,ln_incpc_only, ztscale, zwgt
  
 
@@ -286 +286 @@
     WRITE(numout,*) '     evolve pcbias at each timestep ln_itdecay     = ',ln_itdecay
     WRITE(numout,*) '     incremental press. correction  ln_incpc       = ',ln_incpc
+    WRITE(numout,*) '     incremental press. correction only  ln_incpc_only       = ',ln_incpc_only
     WRITE(numout,*) '     incremental press. correction  zwgt           = ',zwgt
     WRITE(numout,*) '     incremental press. correction  ztscale        = ',ztscale
@@ -312 +313 @@
          IF ( (.NOT. ln_itdecay) .AND. ln_incpc) &   
             &   CALL ctl_stop (' if you set ln_incpc to .true. then you need to set ln_itdecay to .true. as well' )
+
+         IF ( (.NOT. ln_incpc) .AND. ln_incpc_only) &   
+            &   CALL ctl_stop (' if you set ln_incpc_only to .true. then you need to set ln_incpc to .true. as well' )
 
       ENDIF
@@ -951 +955 @@
             tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_p(:,:,:) - tbias_i(:,:,:)
             tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:)
+            IF ( ln_incpc_only ) THEN
+               tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_i(:,:,:)
+               tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_i(:,:,:)
+            ENDIF
          ENDIF
       ELSE
@@ -958 +966 @@
             tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_p(:,:,:) - tbias_i(:,:,:)
             tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:)
+            IF ( ln_incpc_only ) THEN
+               tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_i(:,:,:)
+               tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_i(:,:,:)
+            ENDIF
          ENDIF
       ENDIF