Changeset 7653
 Timestamp:
 20170207T16:37:29+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
r7652 r7653 404 404 tbias_asm(:,:,:) = 0.0_wp 405 405 sbias_asm(:,:,:) = 0.0_wp 406 tbias_asm_out(:,:,:) = 0.0_wp 407 sbias_asm_out(:,:,:) = 0.0_wp 406 408 ENDIF 407 409 … … 409 411 tbias_asm_stscale(:,:,:) = 0.0_wp 410 412 sbias_asm_stscale(:,:,:) = 0.0_wp 413 tbias_asm_stscale_out(:,:,:) = 0.0_wp 414 sbias_asm_stscale_out(:,:,:) = 0.0_wp 411 415 ENDIF 412 416 … … 645 649 & ( t_asm_mem**(real(kt)/ztsday) * sbias_asm(:,:,jk) + & 646 650 & ztfrac * t_asm_upd * s_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof1(:,:) 647 tbias_p(:,:,jk) = tbias_p(:,:,jk) + & 651 652 tbias_p(:,:,jk) = tbias_p(:,:,jk) + & 648 653 & ( t_asm_mem**(real(kt)/ztsday) * tbias_asm(:,:,jk) + & 649 654 & ztfrac * t_asm_upd * t_bkginc(:,:,jk) * tmask(:,:,jk) ) * zcof2(:,:) … … 907 912 908 913 IF ( kt == nn_bias_itwrt ) THEN 909 tbias_p_out(:,:,:) = tbias_p(:,:,:) 910 sbias_p_out(:,:,:) = sbias_p(:,:,:) 911 IF (ln_incpc) THEN 912 tbias_i_out(:,:,:) = tbias_i(:,:,:) 913 sbias_i_out(:,:,:) = sbias_i(:,:,:) 914 915 916 IF ( ln_bias_pc_app ) THEN 917 tbias_p_out(:,:,:) = tbias_p(:,:,:) 918 sbias_p_out(:,:,:) = sbias_p(:,:,:) 919 IF (ln_incpc) THEN 920 tbias_i_out(:,:,:) = tbias_i(:,:,:) 921 sbias_i_out(:,:,:) = sbias_i(:,:,:) 922 IF (ln_incpc_only) THEN 923 IF(lwp) WRITE(numout,*) ' t/p bias_p_out and tbias_asm_out are set to zero if incpc only' 924 tbias_p_out(:,:,:) = 0.0 925 sbias_p_out(:,:,:) = 0.0 926 tbias_asm_out(:,:,:) = 0.0 927 sbias_asm_out(:,:,:) = 0.0 928 ENDIF 929 ENDIF 930 931 ELSE 932 933 IF(lwp) WRITE(numout,*) ' ln_bias_pc_app is =',ln_bias_pc_app, 'so t/p bias_p_out and tbias_asm_out are set to zero' 934 tbias_p_out(:,:,:) = 0.0 935 sbias_p_out(:,:,:) = 0.0 936 tbias_asm_out(:,:,:) = 0.0 937 sbias_asm_out(:,:,:) = 0.0 938 IF (ln_incpc) THEN 939 IF(lwp) WRITE(numout,*) ' ln_bias_pc_app is =',ln_bias_pc_app,' ln_incpc =',ln_incpc,'so t/p bias_i_out and tbias_asm_stscale_out are set to zero' 940 tbias_i_out(:,:,:) = 0.0 941 sbias_i_out(:,:,:) = 0.0 942 tbias_asm_stscale_out(:,:,:) = 0.0 943 sbias_asm_stscale_out(:,:,:) = 0.0 944 ENDIF 914 945 ENDIF 915 946 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.