Changeset 13214 for NEMO/trunk/src/OCE
- Timestamp:
- 2020-07-02T11:09:01+02:00 (4 years ago)
- Location:
- NEMO/trunk/src/OCE
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/IOM/iom.F90
r12649 r13214 1193 1193 & 'we accept this case, even if there is a possible mix-up between z and time dimension' ) 1194 1194 idmspc = idmspc - 1 1195 ELSE 1196 CALL ctl_stop( TRIM(clinfo), 'To keep iom lisibility, when reading a '//clrankpv//'D array,' , & 1197 & 'we do not accept data with '//cldmspc//' spatial dimensions', & 1198 & 'Use ncwa -a to suppress the unnecessary dimensions' ) 1195 !!GS: possibility to read 3D ABL atmopsheric forcing and use 1st level to force BULK simulation 1196 !ELSE 1197 ! CALL ctl_stop( TRIM(clinfo), 'To keep iom lisibility, when reading a '//clrankpv//'D array,', & 1198 ! & 'we do not accept data with '//cldmspc//' spatial dimensions' , & 1199 ! & 'Use ncwa -a to suppress the unnecessary dimensions' ) 1199 1200 ENDIF 1200 1201 ENDIF -
NEMO/trunk/src/OCE/SBC/sbcblk.F90
r13208 r13214 112 112 REAL(wp) :: rn_zu ! z(u) : height of wind measurements 113 113 ! 114 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: Cd_ice , Ch_ice , Ce_ice ! transfert coefficients over ice115 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: Cdn_oce, Chn_oce, Cen_oce ! neutral coeffs over ocean (L15 bulk scheme)116 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: t_zu, q_zu ! air temp. and spec. hum. at wind speed height (L15 bulk scheme)114 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: Cdn_oce, Chn_oce, Cen_oce ! neutral coeffs over ocean (L15 bulk scheme and ABL) 115 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: Cd_ice , Ch_ice , Ce_ice ! transfert coefficients over ice 116 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: t_zu, q_zu ! air temp. and spec. hum. at wind speed height (L15 bulk scheme) 117 117 118 118 LOGICAL :: ln_skin_cs ! use the cool-skin (only available in ECMWF and COARE algorithms) !LB … … 121 121 LOGICAL :: ln_humi_dpt ! humidity read in files ("sn_humi") is dew-point temperature [K] if .true. !LB 122 122 LOGICAL :: ln_humi_rlh ! humidity read in files ("sn_humi") is relative humidity [%] if .true. !LB 123 LOGICAL :: ln_tpot !!GS: flag to compute or not potential temperature 123 124 ! 124 125 INTEGER :: nhumi ! choice of the bulk algorithm … … 181 182 & ln_NCAR, ln_COARE_3p0, ln_COARE_3p6, ln_ECMWF, & ! bulk algorithm 182 183 & cn_dir , rn_zqt, rn_zu, & 183 & rn_pfac, rn_efac, ln_Cd_L12, ln_Cd_L15, 184 & rn_pfac, rn_efac, ln_Cd_L12, ln_Cd_L15, ln_tpot, & 184 185 & ln_crt_fbk, rn_stau_a, rn_stau_b, & ! current feedback 185 186 & ln_skin_cs, ln_skin_wl, ln_humi_sph, ln_humi_dpt, ln_humi_rlh ! cool-skin / warm-layer !LB … … 626 627 !#LB: because AGRIF hates functions that return something else than a scalar, need to 627 628 ! use scalar version of gamma_moist() ... 628 DO_2D_11_11 629 ztpot(ji,jj) = ptair(ji,jj) + gamma_moist( ptair(ji,jj), zqair(ji,jj) ) * rn_zqt 630 END_2D 629 IF( ln_tpot ) THEN 630 DO_2D_11_11 631 ztpot(ji,jj) = ptair(ji,jj) + gamma_moist( ptair(ji,jj), zqair(ji,jj) ) * rn_zqt 632 END_2D 633 ELSE 634 ztpot = ptair(:,:) 635 ENDIF 631 636 ENDIF 632 637 … … 970 975 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=putaui , clinfo1=' blk_ice: putaui : ' & 971 976 & , tab2d_2=pvtaui , clinfo2=' pvtaui : ' ) 972 ELSE 977 ELSE ! ln_abl 973 978 zztmp1 = 11637800.0_wp 974 979 zztmp2 = -5897.8_wp
Note: See TracChangeset
for help on using the changeset viewer.