Changeset 2748
- Timestamp:
- 2011-04-19T09:04:38+02:00 (13 years ago)
- File:
-
- 1 edited
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r2715 r2748 755 755 !!---------------------------------------------------------------------- 756 756 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released, iwrk_in_use, iwrk_not_released 757 USE wrk_nemo, ONLY: dU10 => wrk_2d_1 ! dU [m/s]758 USE wrk_nemo, ONLY: dT => wrk_2d_ 2! air/sea temperature difference [K]759 USE wrk_nemo, ONLY: dq => wrk_2d_ 3! air/sea humidity difference [K]760 USE wrk_nemo, ONLY: Cd_n10 => wrk_2d_ 4! 10m neutral drag coefficient761 USE wrk_nemo, ONLY: Ce_n10 => wrk_2d_ 5! 10m neutral latent coefficient762 USE wrk_nemo, ONLY: Ch_n10 => wrk_2d_ 6! 10m neutral sensible coefficient763 USE wrk_nemo, ONLY: sqrt_Cd_n10 => wrk_2d_ 7! root square of Cd_n10764 USE wrk_nemo, ONLY: sqrt_Cd => wrk_2d_ 8! root square of Cd765 USE wrk_nemo, ONLY: T_vpot => wrk_2d_ 9! virtual potential temperature [K]766 USE wrk_nemo, ONLY: T_star => wrk_2d_ 10! turbulent scale of tem. fluct.767 USE wrk_nemo, ONLY: q_star => wrk_2d_ 11! turbulent humidity of temp. fluct.768 USE wrk_nemo, ONLY: U_star => wrk_2d_ 12! turb. scale of velocity fluct.769 USE wrk_nemo, ONLY: L => wrk_2d_ 13! Monin-Obukov length [m]770 USE wrk_nemo, ONLY: zeta_u => wrk_2d_ 14! stability parameter at height zu771 USE wrk_nemo, ONLY: zeta_t => wrk_2d_ 15! stability parameter at height zt772 USE wrk_nemo, ONLY: U_n10 => wrk_2d_ 16! neutral wind velocity at 10m [m]773 USE wrk_nemo, ONLY: xlogt => wrk_2d_ 17, xct => wrk_2d_18, zpsi_hu => wrk_2d_19, zpsi_ht => wrk_2d_20, zpsi_m => wrk_2d_21757 USE wrk_nemo, ONLY: dU10 => wrk_2d_14 ! dU [m/s] 758 USE wrk_nemo, ONLY: dT => wrk_2d_15 ! air/sea temperature difference [K] 759 USE wrk_nemo, ONLY: dq => wrk_2d_16 ! air/sea humidity difference [K] 760 USE wrk_nemo, ONLY: Cd_n10 => wrk_2d_17 ! 10m neutral drag coefficient 761 USE wrk_nemo, ONLY: Ce_n10 => wrk_2d_18 ! 10m neutral latent coefficient 762 USE wrk_nemo, ONLY: Ch_n10 => wrk_2d_19 ! 10m neutral sensible coefficient 763 USE wrk_nemo, ONLY: sqrt_Cd_n10 => wrk_2d_20 ! root square of Cd_n10 764 USE wrk_nemo, ONLY: sqrt_Cd => wrk_2d_21 ! root square of Cd 765 USE wrk_nemo, ONLY: T_vpot => wrk_2d_22 ! virtual potential temperature [K] 766 USE wrk_nemo, ONLY: T_star => wrk_2d_23 ! turbulent scale of tem. fluct. 767 USE wrk_nemo, ONLY: q_star => wrk_2d_24 ! turbulent humidity of temp. fluct. 768 USE wrk_nemo, ONLY: U_star => wrk_2d_25 ! turb. scale of velocity fluct. 769 USE wrk_nemo, ONLY: L => wrk_2d_26 ! Monin-Obukov length [m] 770 USE wrk_nemo, ONLY: zeta_u => wrk_2d_27 ! stability parameter at height zu 771 USE wrk_nemo, ONLY: zeta_t => wrk_2d_28 ! stability parameter at height zt 772 USE wrk_nemo, ONLY: U_n10 => wrk_2d_29 ! neutral wind velocity at 10m [m] 773 USE wrk_nemo, ONLY: xlogt => wrk_2d_30, xct => wrk_2d_31, zpsi_hu => wrk_2d_32, zpsi_ht => wrk_2d_33, zpsi_m => wrk_2d_34 774 774 USE wrk_nemo, ONLY: stab => iwrk_2d_1 ! 1st guess stability test integer 775 775 !! … … 798 798 !! * Start 799 799 800 IF( wrk_in_use(2, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21) .OR. & 801 iwrk_in_use(2, 1) ) THEN 800 IF( wrk_in_use(2, 14,15,16,17,18,19, & 801 20,21,22,23,24,25,26,27,28,29, & 802 30,31,32,33,34) .OR. & 803 iwrk_in_use(2, 1) ) THEN 802 804 CALL ctl_stop('TURB_CORE_2Z: requested workspace arrays unavailable') ; RETURN 803 END 805 ENDIF 804 806 805 807 !! Initial air/sea differences … … 876 878 END DO 877 879 !! 878 IF( wrk_not_released(2, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21) .OR. & 879 iwrk_not_released(2, 1) ) CALL ctl_stop('TURB_CORE_2Z: requested workspace arrays unavailable') 880 IF( wrk_not_released(2, 14,15,16,17,18,19, & 881 & 20,21,22,23,24,25,26,27,28,29, & 882 & 30,31,32,33,34 ) .OR. & 883 iwrk_not_released(2, 1) ) & 884 CALL ctl_stop('TURB_CORE_1Z: failed to release workspace arrays') 880 885 ! 881 886 END SUBROUTINE TURB_CORE_2Z … … 885 890 !------------------------------------------------------------------------------- 886 891 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 887 USE wrk_nemo, ONLY: X2 => wrk_2d_3 3888 USE wrk_nemo, ONLY: X => wrk_2d_3 4889 USE wrk_nemo, ONLY: stabit => wrk_2d_3 5892 USE wrk_nemo, ONLY: X2 => wrk_2d_35 893 USE wrk_nemo, ONLY: X => wrk_2d_36 894 USE wrk_nemo, ONLY: stabit => wrk_2d_37 890 895 !! 891 896 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: zta … … 895 900 !------------------------------------------------------------------------------- 896 901 897 IF( wrk_in_use(2, 3 3,34,35) ) THEN902 IF( wrk_in_use(2, 35,36,37) ) THEN 898 903 CALL ctl_stop('psi_m: requested workspace arrays unavailable') ; RETURN 899 904 ENDIF … … 904 909 & + (1. - stabit)*(2*log((1. + X)/2) + log((1. + X2)/2) - 2*atan(X) + pi/2) ! Unstable 905 910 906 IF( wrk_not_released(2, 3 3,34,35) ) CALL ctl_stop('psi_m: failed to release workspace arrays')911 IF( wrk_not_released(2, 35,36,37) ) CALL ctl_stop('psi_m: failed to release workspace arrays') 907 912 ! 908 913 END FUNCTION psi_m … … 912 917 !------------------------------------------------------------------------------- 913 918 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 914 USE wrk_nemo, ONLY: X2 => wrk_2d_3 3915 USE wrk_nemo, ONLY: X => wrk_2d_3 4916 USE wrk_nemo, ONLY: stabit => wrk_2d_3 5919 USE wrk_nemo, ONLY: X2 => wrk_2d_35 920 USE wrk_nemo, ONLY: X => wrk_2d_36 921 USE wrk_nemo, ONLY: stabit => wrk_2d_37 917 922 ! 918 923 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: zta … … 921 926 !------------------------------------------------------------------------------- 922 927 923 IF( wrk_in_use(2, 3 3,34,35) ) THEN928 IF( wrk_in_use(2, 35,36,37) ) THEN 924 929 CALL ctl_stop('psi_h: requested workspace arrays unavailable') ; RETURN 925 930 ENDIF … … 930 935 & + (1. - stabit)*(2.*log( (1. + X2)/2. )) ! Unstable 931 936 932 IF( wrk_not_released(2, 3 3,34,35) ) CALL ctl_stop('psi_h: failed to release workspace arrays')937 IF( wrk_not_released(2, 35,36,37) ) CALL ctl_stop('psi_h: failed to release workspace arrays') 933 938 ! 934 939 END FUNCTION psi_h
Note: See TracChangeset
for help on using the changeset viewer.
