Changeset 1126 for trunk/LMDZ.COMMON/libf/dyn3dpar
- Timestamp:
- Dec 17, 2013, 1:02:44 PM (11 years ago)
- Location:
- trunk/LMDZ.COMMON/libf/dyn3dpar
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/dyn3dpar/calfis_p.F
r1107 r1126 260 260 REAL,allocatable,save :: tmpvarbar(:) 261 261 REAL,allocatable,save :: tmpvarbarp1(:) 262 real :: zz1,zz2 262 263 263 264 c----------------------------------------------------------------------- … … 314 315 c------------------------------------------------------------------ 315 316 c moyennes globales pour les profils de pression et de temperature 316 if(planet_type.eq."titan" ) then317 if(planet_type.eq."titan".or.planet_type.eq."venus") then 317 318 call AllGather_Field(pp,iip1*jjp1,llmp1) 318 319 call AllGather_Field(pteta,iip1*jjp1,llm) 319 320 call AllGather_Field(ppk,iip1*jjp1,llm) 321 call AllGather_Field(pphi,iip1*jjp1,llm) 322 call AllGather_Field(pphis,iip1*jjp1,1) 320 323 ALLOCATE(plevmoy(llm+1)) 321 324 ALLOCATE(playmoy(llm)) … … 323 326 ALLOCATE(tetamoy(llm)) 324 327 ALLOCATE(pkmoy(llm)) 328 ALLOCATE(phimoy(0:llm)) 329 ALLOCATE(zlevmoy(llm+1)) 330 ALLOCATE(zlaymoy(llm)) 325 331 plevmoy=0. 326 332 do l=1,llmp1 … … 333 339 tetamoy=0. 334 340 pkmoy=0. 341 phimoy=0. 342 do i=1,iip1 343 do j=1,jjp1 344 phimoy(0)=phimoy(0)+pphis(i,j)/(iip1*jjp1) 345 enddo 346 enddo 335 347 do l=1,llm 336 348 do i=1,iip1 … … 338 350 tetamoy(l)=tetamoy(l)+pteta(i,j,l)/(iip1*jjp1) 339 351 pkmoy(l)=pkmoy(l)+ppk(i,j,l)/(iip1*jjp1) 352 phimoy(l)=phimoy(l)+pphi(i,j,l)/(iip1*jjp1) 340 353 enddo 341 354 enddo 342 355 enddo 343 playmoy = preff * (pkmoy/cpp) ** unskap356 playmoy(:) = preff * (pkmoy(:)/cpp) ** unskap 344 357 call tpot2t_p(1,llm,tetamoy,tmoy,pkmoy) 358 c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE: 359 zlaymoy(:) = g*rad*rad/(g*rad-phimoy(:))-rad 360 zlevmoy(1) = phimoy(0)/g 361 DO l=2,llm 362 zz1=(playmoy(l-1)+plevmoy(l))/(playmoy(l-1)-plevmoy(l)) 363 zz2=(plevmoy(l) +playmoy(l))/(plevmoy(l) -playmoy(l)) 364 zlevmoy(l)=(zz1*zlaymoy(l-1)+zz2*zlaymoy(l))/(zz1+zz2) 365 ENDDO 366 zlevmoy(llmp1)=zlaymoy(llm)+(zlaymoy(llm)-zlevmoy(llm)) 345 367 c------------------- 346 368 c + lat index 347 369 allocate(klat(klon)) 348 370 do ig0=1,klon 349 350 371 j=index_j(ig0) 372 klat(ig0)=j 351 373 enddo 352 374 endif ! planet_type=titan -
trunk/LMDZ.COMMON/libf/dyn3dpar/moyzon_mod.F90
r1107 r1126 24 24 ! to be changed... 25 25 REAL,ALLOCATABLE,SAVE :: plevmoy(:),playmoy(:) 26 REAL,ALLOCATABLE,SAVE :: zlevmoy(:),zlaymoy(:),phimoy(:) 26 27 REAL,ALLOCATABLE,SAVE :: tmoy(:),tetamoy(:),pkmoy(:) 27 28 INTEGER,ALLOCATABLE,SAVE :: klat(:)
Note: See TracChangeset
for help on using the changeset viewer.