 Timestamp:
 20170209T13:37:38+01:00 (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90
r7660 r7661 811 811 IF (lwp) WRITE(numout,*) 'obsoper or forecast mode  IPC  computing tbias_i and sbias_i' 812 812 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(:,:)) ) 815 815 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) 816 820 ENDIF 817 821 … … 967 971 tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal)  sbias_p(:,:,:)  sbias_i(:,:,:) 968 972 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' 969 975 IF(lwp) WRITE(numout,*) 'ln_incpc_only =', ln_incpc_only, 'tsw updated with IPC only and ln_dynhpg_imp = ',ln_dynhpg_imp 970 976 tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem)  tbias_i(:,:,:) … … 978 984 979 985 ! is this needed? 986 IF(lwp) WRITE(numout,*) 'CALL lbc_lnk( rhd_pc, T, 1.0_wp )' 980 987 CALL lbc_lnk( rhd_pc, 'T', 1.0_wp ) 981 988 982 989 ! Partial steps: now horizontal gradient of t,s,rd 983 990 ! at the bottom ocean level 991 IF(lwp) WRITE(numout,*) 'horizontal gradient of t,s,rd ' 984 992 IF( ln_zps ) THEN 985 993 CALL zps_hde( kt, jpts, tsw, gtsu, gtsv, &
Note: See TracChangeset
for help on using the changeset viewer.