 Timestamp:
 20170208T14:55:40+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
r7653 r7656 771 771 772 772 IF ( ln_incpc ) THEN 773 IF (lwp) WRITE(numout,*) 'after end of IAU  IPC  saving s/tbias_asm_stscale' 773 774 DO jk = 1, jpkm1 774 775 tbias_asm_stscale_out(:,:,jk) = tbias_asm_stscale(:,:,jk) * zfrac … … 785 786 786 787 ELSE ! no assimilation increments, simply decay pressure correction (e.g for forecasts) 788 ! this occurs at obsoper step as well ( ln_asmiau is false ) 787 789 788 790 DO jk = 1, jpkm1 … … 806 808 807 809 IF (ln_incpc) THEN 808 810 IF (lwp) WRITE(numout,*) 'obsoper or forecast mode  IPC  computing tbias_i and sbias_i' 809 811 DO jk = 1, jpkm1 810 tbias_i(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0  fbcoef_stscale(:,:)) )811 sbias_i(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0  fbcoef_stscale(:,:)) )812 tbias_i(:,:,jk) = tbias_i(:,:,jk) + ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0  fbcoef_stscale(:,:)) ) 813 sbias_i(:,:,jk) = sbias_i(:,:,jk) + ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0  fbcoef_stscale(:,:)) ) 812 814 ENDDO 813 814 IF ( kt == nn_bias_itwrt ) THEN815 DO jk = 1, jpkm1816 tbias_asm_stscale_out(:,:,jk) = ( tbias_asm_stscale(:,:,jk) * zfrac1 )817 sbias_asm_stscale_out(:,:,jk) = ( sbias_asm_stscale(:,:,jk) * zfrac1 )818 ENDDO819 ENDIF820 821 815 ENDIF 822 816 … … 912 906 913 907 IF ( kt == nn_bias_itwrt ) THEN 914 915 916 IF ( ln_bias_pc_app ) THEN917 908 tbias_p_out(:,:,:) = tbias_p(:,:,:) 918 909 sbias_p_out(:,:,:) = sbias_p(:,:,:) 919 910 IF (ln_incpc) THEN 920 911 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 912 sbias_i_out(:,:,:) = sbias_i(:,:,:) 929 913 ENDIF 930 931 ELSE932 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.0935 sbias_p_out(:,:,:) = 0.0936 tbias_asm_out(:,:,:) = 0.0937 sbias_asm_out(:,:,:) = 0.0938 IF (ln_incpc) THEN939 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.0941 sbias_i_out(:,:,:) = 0.0942 tbias_asm_stscale_out(:,:,:) = 0.0943 sbias_asm_stscale_out(:,:,:) = 0.0944 ENDIF945 ENDIF946 914 ENDIF 947 915
Note: See TracChangeset
for help on using the changeset viewer.