Changeset 7661


Ignore:
Timestamp:
2017-02-09T13:37:38+01:00 (4 years ago)
Author:
isabella
Message:

adding more messages

File:
1 edited

Legend:

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

    r7660 r7661  
    811811                   IF (lwp) WRITE(numout,*) 'obsoper or forecast mode - IPC - computing tbias_i and sbias_i' 
    812812                   DO jk = 1, jpkm1 
    813                       tbias_i(:,:,jk) = (tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ))*tmask(:,:,:) 
    814                       sbias_i(:,:,jk) = (sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ))*tmask(:,:,:) 
     813                      tbias_i(:,:,jk) = tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
     814                      sbias_i(:,:,jk) = sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    815815                   ENDDO 
     816                    
     817                   tbias_i(:,:,:) =tbias_i(:,:,:)*tmask(:,:,:) 
     818                   sbias_i(:,:,:) =sbias_i(:,:,:)*tmask(:,:,:) 
     819                   IF (lwp) WRITE(numout,*)  'SUM(tbias_i)' SUM(tbias_i) 
    816820               ENDIF 
    817821 
     
    967971            tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:) 
    968972            IF ( ln_incpc_only ) THEN 
     973               IF (lwp) WRITE(numout,*) 'SUM(tbias_i) is' SUM(tbias_i) 'before to be added to tsw' 
     974               IF (lwp) WRITE(numout,*) 'SUM(sbias_i) is' SUM(sbias_i) 'before to be added to tsw' 
    969975               IF(lwp) WRITE(numout,*) 'ln_incpc_only =', ln_incpc_only, 'tsw updated with IPC only and ln_dynhpg_imp = ',ln_dynhpg_imp 
    970976               tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_i(:,:,:) 
     
    978984       
    979985      ! is this needed? 
     986      IF(lwp) WRITE(numout,*) 'CALL lbc_lnk( rhd_pc, T, 1.0_wp )' 
    980987      CALL lbc_lnk( rhd_pc, 'T', 1.0_wp ) 
    981988 
    982989      ! Partial steps: now horizontal gradient of t,s,rd  
    983990      ! at the bottom ocean level 
     991      IF(lwp) WRITE(numout,*) 'horizontal gradient of t,s,rd ' 
    984992      IF( ln_zps    )  THEN 
    985993         CALL zps_hde( kt, jpts, tsw, gtsu, gtsv,  & 
Note: See TracChangeset for help on using the changeset viewer.