Changeset 204 for LMDZ.3.3/trunk/libf/phylmd/phystokenc.F
- Timestamp:
- Apr 13, 2001, 12:44:53 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/phylmd/phystokenc.F
r199 r204 6 6 I pt,pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, 7 7 I pcoefh,yu1,yv1,ftsol,pctsrf, 8 I frac_impa,frac_nucl,8 I pfrac_impa,pfrac_nucl, 9 9 I pphis,paire,dtime,itap) 10 10 USE ioipsl … … 63 63 c ---------- 64 64 c 65 REAL frac_impa(klon,klev)66 REAL frac_nucl(klon,klev)65 REAL pfrac_impa(klon,klev) 66 REAL pfrac_nucl(klon,klev) 67 67 c 68 68 c Arguments necessaires pour les sources et puits de traceur … … 82 82 REAL coefh(klon,klev) ! flux detraine dans le panache descendant 83 83 REAL t(klon,klev) 84 REAL frac_impa(klon,klev) 85 REAL frac_nucl(klon,klev) 86 REAL rain(klon) 84 87 85 88 REAL pyu1(klon),pyv1(klon) … … 95 98 save t,mfu,mfd,en_u,de_u,en_d,de_d,coefh,dtcum 96 99 save iadvtr,irec 100 save frac_impa,frac_nucl,rain 97 101 save pyu1,pyv1,pftsol,ppsrf 98 102 … … 133 137 do k=1,klev 134 138 do i=1,klon 139 frac_impa(i,k)=1. 140 frac_nucl(i,k)=1. 135 141 mfu(i,k)=0. 136 142 mfd(i,k)=0. … … 144 150 enddo 145 151 do i=1,klon 152 rain(i)=0. 146 153 pyv1(i)=0. 147 154 pyu1(i)=0. … … 159 166 do k=1,klev 160 167 do i=1,klon 168 frac_impa(i,k)=frac_impa(i,k)*pfrac_impa(i,k) 169 frac_nucl(i,k)=frac_nucl(i,k)*pfrac_nucl(i,k) 161 170 mfu(i,k)=mfu(i,k)+pmfu(i,k)*pdtphys 162 171 mfd(i,k)=mfd(i,k)+pmfd(i,k)*pdtphys … … 187 196 do k=1,klev 188 197 do i=1,klon 198 c frac_impa=frac_impa : c'est la fraction cumulee qu'on stoke 199 c frac_nucl=frac_nucl : c'est la fraction cumulee qu'on stoke 189 200 mfu(i,k)=mfu(i,k)/dtcum 190 201 mfd(i,k)=mfd(i,k)/dtcum … … 198 209 enddo 199 210 do i=1,klon 211 rain(i)=rain(i)/dtcum 200 212 pyv1(i)=pyv1(i)/dtcum 201 213 pyu1(i)=pyu1(i)/dtcum … … 226 238 227 239 ccccc 240 print*,'AVANT ECRITURE' 228 241 CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d) 229 242 CALL histwrite(physid,"t",itap,zx_tmp_3d, 230 243 . iim*(jjm+1)*klev,ndex) 244 print*,'APRES ECRITURE' 231 245 232 246 CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d) … … 279 293 CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d) 280 294 CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d, 295 . iim*(jjm+1),ndex) 296 297 CALL gr_fi_ecrit(1,klon,iim,jjm+1, rain, zx_tmp_2d) 298 CALL histwrite(physid,"rain",itap,zx_tmp_2d, 281 299 . iim*(jjm+1),ndex) 282 300
Note: See TracChangeset
for help on using the changeset viewer.