Changeset 7662 for branches/UKMO
- Timestamp:
- 2017-02-09T16:07:37+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90
r7661 r7662 814 814 sbias_i(:,:,jk) = sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 815 815 ENDDO 816 817 816 tbias_i(:,:,:) =tbias_i(:,:,:)*tmask(:,:,:) 818 817 sbias_i(:,:,:) =sbias_i(:,:,:)*tmask(:,:,:) 819 IF (lwp) WRITE(numout,*) 'SUM(tbias_i)' SUM(tbias_i)820 818 ENDIF 821 822 819 ENDIF 823 820 … … 971 968 tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:) 972 969 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' 970 IF(lwp) WRITE(numout,*) 'SUM(tbias_i) is', SUM(tbias_i), 'before to be added to tsw' 975 971 IF(lwp) WRITE(numout,*) 'ln_incpc_only =', ln_incpc_only, 'tsw updated with IPC only and ln_dynhpg_imp = ',ln_dynhpg_imp 976 972 tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_i(:,:,:) … … 979 975 ENDIF 980 976 ENDIF 977 978 981 979 982 980 IF(lwp) WRITE(numout,*) 'dyn_bias( kt ) calculating rhd_pc, kt =', kt
Note: See TracChangeset
for help on using the changeset viewer.