Changeset 7186
- Timestamp:
- 2016-11-03T17:15:25+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90
r7172 r7186 151 151 & sbias_i_out, & !: Output salinity bias field for incremental P correction 152 152 & tbias_asm_stscale, & !: Short time scale temperature bias field 153 & sbias_asm_stscale !: Short time scale salinity bias field 153 & sbias_asm_stscale, & !: Short time scale salinity bias field 154 & tbias_asm_stscale_out & !: Short time scale temperature bias output field 155 & sbias_asm_stscale_out !: Short time scale salinity bias output field 154 156 155 157 INTEGER, PRIVATE :: nn_lat_ramp ! choice of latitude dependent ramp … … 336 338 337 339 IF( ln_incpc ) ALLOCATE( tbias_asm_stscale(jpi,jpj,jpk), & 338 & sbias_asm_stscale(jpi,jpj,jpk)) 340 & sbias_asm_stscale(jpi,jpj,jpk), & 341 tbias_asm_stscale_out(jpi,jpj,jpk), & 342 sbias_asm_stscale_out(jpi,jpj,jpk)) 339 343 340 344 lalloc = .TRUE. … … 648 652 & - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 649 653 ENDDO 650 654 655 IF ( kt == nn_bias_itwrt ) THEN 656 DO jk = 1, jpk 657 tbias_asm_stscale_out(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * ztfrac ) + ( tbias_asm_stscale(:,:,jk) * zfrac1 ) 658 sbias_asm_stscale_out(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * ztfrac ) + ( sbias_asm_stscale(:,:,jk) * zfrac1 ) 659 ENDDO 660 ENDIF 651 661 ENDIF 652 662 … … 671 681 END DO 672 682 673 IF (ln_incpc) THEN 674 DO jk = 1, jpk 675 tbias_asm_stscale(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * ztfrac ) + ( tbias_asm_stscale(:,:,jk) * zfrac1 ) 676 sbias_asm_stscale(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * ztfrac ) + ( sbias_asm_stscale(:,:,jk) * zfrac1 ) 677 ENDDO 678 ENDIF 683 679 684 ENDIF 685 680 686 681 687 ELSE ! decay pressure correction from combined historical component and increments after IAU … … 706 712 & - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 707 713 ENDDO 714 715 IF ( kt == nn_bias_itwrt ) THEN 716 DO jk = 1, jpk 717 tbias_asm_stscale_out(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * zfrac ) + ( tbias_asm_stscale(:,:,jk) * zfrac1 ) 718 sbias_asm_stscale_out(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * zfrac ) + ( sbias_asm_stscale(:,:,jk) * zfrac1 ) 719 ENDDO 720 ENDIF 708 721 709 722 ENDIF … … 973 986 974 987 IF ( ln_incpc ) THEN 975 CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'tbias_asm_stscale' , tbias_asm_stscale )976 CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_asm_stscale' , sbias_asm_stscale )988 CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'tbias_asm_stscale' , tbias_asm_stscale_out ) 989 CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_asm_stscale' , sbias_asm_stscale_out ) 977 990 ENDIF 978 991
Note: See TracChangeset
for help on using the changeset viewer.