Changeset 3158 for trunk/LMDZ.MARS/libf/aeronomars
- Timestamp:
- Dec 13, 2023, 4:09:00 PM (13 months ago)
- Location:
- trunk/LMDZ.MARS/libf/aeronomars
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/conduction.F ¶
r2615 r3158 82 82 83 83 zt(1)=pt(ig,1)+pdt(ig,1)*ptimestep 84 c zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g85 c zlev(1)=0.086 84 zlay(1)=zzlay(ig,1) 87 85 zlev(1)=zzlev(ig,1) … … 90 88 91 89 zt(i)=pt(ig,i)+pdt(ig,i)*ptimestep 92 c tmean=zt(i) 93 c if(zt(i).ne.zt(i-1)) 94 c & tmean=(zt(i)-zt(i-1))/log(zt(i)/zt(i-1)) 95 c zlay(i)= zlay(i-1) 96 c & -log(pplay(ig,i)/pplay(ig,i-1))*Rnew(ig,i-1)*tmean/g 97 c zlev(i)= zlev(i-1) 98 c & -log(pplev(ig,i)/pplev(ig,i-1))*Rnew(ig,i-1)*tmean/g 99 zlay(i)=zzlay(ig,i) 100 zlev(i)=zzlev(ig,i) 90 zlay(i)=zzlay(ig,i) 91 zlev(i)=zzlev(ig,i) 101 92 enddo 102 103 c zlev(nlayer+1)= zlev(nlayer) 104 c & -log(max(pplev(ig,nlayer+1),1.e-30)/pplev(ig,nlayer)) 105 c & *Rnew(ig,nlayer)*tmean/g 106 c if(pplev(ig,nlayer+1).eq.0.) 107 c & zlev(nlayer+1)=zlev(nlayer)+(zlay(nlayer)-zlay(nlayer-1)) 108 93 109 94 zlev(nlayer+1)= zlev(nlayer)+10000. 110 95 -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90 ¶
r2615 r3158 411 411 enddo 412 412 413 ! zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))414 ! & *Rnew(ig,1)*zt(ig,1)/g415 413 zlocal(1)=zzlay(ig,1) 416 414 zlocal(1)=zlocal(1)/1000. -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/moldiff.F ¶
r2615 r3158 308 308 c 309 309 310 c zlocal(1)=-log(pplay(ig,1)/pplev(ig,1))* Rnew(ig,1)*tt(1)/g311 310 zlocal(1)=zzlay(ig,1) 312 311 … … 325 324 enddo 326 325 327 c tmean=tt(l)328 c if(tt(l).ne.tt(l-1))329 c & tmean=(tt(l)-tt(l-1))/log(tt(l)/tt(l-1))330 c zlocal(l)= zlocal(l-1)331 c & -log(pplay(ig,l)/pplay(ig,l-1))*rnew(ig,l)*tmean/g332 326 zlocal(l)=zzlay(ig,l) 333 327 enddo 334 328 335 c zlocal(nz)= zlocal(nz-1)336 c & -log(pplay(ig,nz)/pplay(ig,nz-1))*rnew(ig,nz)*tmean/g337 329 zlocal(nz)=zzlay(ig,nz) 338 330 -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/molvis.F ¶
r2615 r3158 95 95 zt(1)=pt(ig,1)+(pdteuv(ig,1)+pdtconduc(ig,1))*ptimestep 96 96 zvel(1)=pvel(ig,1) 97 c zlay(1)=-log(pplay(ig,1)/pplev(ig,1))*Rnew(ig,1)*zt(1)/g98 c zlev(1)=0.099 100 97 zlay(1)=zzlay(ig,1) 101 98 zlev(1)=zzlev(ig,1) … … 104 101 zt(l)=pt(ig,l)+(pdteuv(ig,l)+pdtconduc(ig,l))*ptimestep 105 102 zvel(l)=pvel(ig,l) 106 c tmean=zt(l) 107 c if(zt(l).ne.zt(l-1)) tmean=(zt(l)-zt(l-1))/log(zt(l)/zt(l-1)) 108 c zlay(l)= zlay(l-1) 109 c & -log(pplay(ig,l)/pplay(ig,l-1))*Rnew(ig,l-1)*tmean/g 110 c zlev(l)= zlev(l-1) 111 c & -log(pplev(ig,l)/pplev(ig,l-1))*Rnew(ig,l-1)*tmean/g 112 zlay(l)=zzlay(ig,l) 113 zlev(l)=zzlev(ig,l) 103 zlay(l)=zzlay(ig,l) 104 zlev(l)=zzlev(ig,l) 114 105 enddo 115 116 c zlev(nz+1)= zlev(nz)117 c & -log(max(pplev(ig,nz+1),1.e-30)/pplev(ig,nz))118 c & *Rnew(ig,nz)*tmean/g119 c if(pplev(ig,nz+1).eq.0.)120 c & zlev(nz+1)=zlev(nz)+(zlay(nz)-zlay(nz-1))121 106 122 107 zlev(nz+1)= zlev(nz)+10000.
Note: See TracChangeset
for help on using the changeset viewer.