Changeset 7661 for branches/UKMO

Timestamp:

2017-02-09T13:37:38+01:00 (7 years ago)

Author:

isabella

Message:

adding more messages

File:

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

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

@@ -811 +811 @@
                    IF (lwp) WRITE(numout,*) 'obsoper or forecast mode - IPC - computing tbias_i and sbias_i'
                    DO jk = 1, jpkm1
-                      tbias_i(:,:,jk) = (tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ))*tmask(:,:,:)
-                      sbias_i(:,:,jk) = (sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ))*tmask(:,:,:)
+                      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
+                   
+                   tbias_i(:,:,:) =tbias_i(:,:,:)*tmask(:,:,:)
+                   sbias_i(:,:,:) =sbias_i(:,:,:)*tmask(:,:,:)
+                   IF (lwp) WRITE(numout,*)  'SUM(tbias_i)' SUM(tbias_i)
                ENDIF
 
@@ -967 +971 @@
             tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:)
             IF ( ln_incpc_only ) THEN
+               IF (lwp) WRITE(numout,*) 'SUM(tbias_i) is' SUM(tbias_i) 'before to be added to tsw'
+               IF (lwp) WRITE(numout,*) 'SUM(sbias_i) is' SUM(sbias_i) 'before to be added to tsw'
                IF(lwp) WRITE(numout,*) 'ln_incpc_only =', ln_incpc_only, 'tsw updated with IPC only and ln_dynhpg_imp = ',ln_dynhpg_imp
                tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_i(:,:,:)
@@ -978 +984 @@
       
       ! is this needed?
+      IF(lwp) WRITE(numout,*) 'CALL lbc_lnk( rhd_pc, T, 1.0_wp )'
       CALL lbc_lnk( rhd_pc, 'T', 1.0_wp )
 
       ! Partial steps: now horizontal gradient of t,s,rd 
       ! at the bottom ocean level
+      IF(lwp) WRITE(numout,*) 'horizontal gradient of t,s,rd '
       IF( ln_zps    )  THEN
          CALL zps_hde( kt, jpts, tsw, gtsu, gtsv,  &