Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/phys_cosp.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/phys_cosp.F90
r5133 r5158 44 44 ! mr_ozone, !Concentration ozone (Kg/Kg) 45 45 ! dem_s !Cloud optical emissivity 46 ! dtau_s 47 ! emsfc_lw = 1. 46 ! dtau_s !Cloud optical thickness 47 ! emsfc_lw = 1. !Surface emissivity dans radlwsw.F90 48 48 49 49 !!! Outputs : … … 132 132 ! Declaration necessaires pour les sorties IOIPSL 133 133 integer :: ii 134 real:: ecrit_day,ecrit_hf,ecrit_mth, missing_val134 REAL :: ecrit_day,ecrit_hf,ecrit_mth, missing_val 135 135 logical :: ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, ok_all_xml 136 136 … … 150 150 real,dimension(Nlevlmdz) :: presnivs 151 151 integer :: itap,k,ip 152 real:: dtime,freq_cosp152 REAL :: dtime,freq_cosp 153 153 real,dimension(2) :: time_bnds 154 154 … … 249 249 250 250 zlev_half(:,1) = phis(:)/9.81 251 dok = 2, Nlevels252 doip = 1, Npoints251 DO k = 2, Nlevels 252 DO ip = 1, Npoints 253 253 zlev_half(ip,k) = phi(ip,k)/9.81 + & 254 254 (phi(ip,k)-phi(ip,k-1))/9.81 * (ph(ip,k)-p(ip,k)) / (p(ip,k)-p(ip,k-1)) … … 267 267 gbx%skt = skt !Skin temperature (K) 268 268 269 doip = 1, Npoints269 DO ip = 1, Npoints 270 270 if (fracTerLic(ip).ge.0.5) then 271 271 gbx%land(ip) = 1.
Note: See TracChangeset
for help on using the changeset viewer.