Changeset 102 for trunk/libf/phytitan/aaam_bud.F

Timestamp:

Mar 24, 2011, 11:46:27 AM (14 years ago)

Author:

slebonnois

Message:

SL : corrections et modifications dans phytitan correspondant a celles
faites apres compilation Venus. Titan pas encore compile.

File:

• trunk/libf/phytitan/aaam_bud.F (modified) (10 diffs)

trunk/libf/phytitan/aaam_bud.F

@@ -2 +2 @@
      i                   rea,rg,ome,      
      i                   plat,plon,phis,
-     i                   dragu,liftu,phyu,
-     i                   dragv,liftv,phyv,
+     i                   dragu,liftu,clu,
+     i                   dragv,liftv,clv,
      i                   p, u, v)
 c
+      use dimphy
       implicit none
 c======================================================================
@@ -17 +18 @@
 C    C               ra,rg,romega,
 C    C               rlat,rlon,pphis,
-C    C               zustrdr,zustrli,zustrph,
-C    C               zvstrdr,zvstrli,zvstrph,
+C    C               zustrdr,zustrli,zustrcl,
+C    C               zvstrdr,zvstrli,zvstrcl,
 C    C               paprs,u,v)
 C
@@ -39 +40 @@
 c dragu---input-R-orodrag stress (zonal)
 c liftu---input-R-orolift stress (zonal)
-c phyu----input-R-Stress total de la physique (zonal)
+c clu-----input-R-Boundary layer stress (zonal)
 c dragv---input-R-orodrag stress (Meridional)
 c liftv---input-R-orolift stress (Meridional)
-c phyv----input-R-Stress total de la physique (Meridional)
+c clv-----input-R-Boundary layer stress (Meridional)
 c p-------input-R-Pressure (Pa) at model half levels
 c u-------input-R-Horizontal wind (m/s)
@@ -78 +79 @@
 
 #include "dimensions.h"
-#include "dimphy.h"
 c
 c ARGUMENTS
@@ -85 +85 @@
       REAL rjour,rsec,rea,rg,ome
       REAL plat(nlon),plon(nlon),phis(nlon)
-      REAL dragu(nlon),liftu(nlon),phyu(nlon)             
-      REAL dragv(nlon),liftv(nlon),phyv(nlon)             
+      REAL dragu(nlon),liftu(nlon),clu(nlon)             
+      REAL dragv(nlon),liftv(nlon),clv(nlon)             
       REAL p(nlon,nlev+1), u(nlon,nlev), v(nlon,nlev)
 c
@@ -148 +148 @@
           ssou(i,1)=dragu(l)+liftu(l)
           ssov(i,1)=dragv(l)+liftv(l)
-          blsu(i,1)=phyu(l)-dragu(l)-liftu(l)
-          blsv(i,1)=phyv(l)-dragv(l)-liftv(l)
+          blsu(i,1)=clu(l)
+          blsv(i,1)=clv(l)
 
         enddo
@@ -162 +162 @@
           ssou(iim+1,j)=dragu(l+1)+liftu(l+1)
           ssov(iim+1,j)=dragv(l+1)+liftv(l+1)
-          blsu(iim+1,j)=phyu(l+1)-dragu(l+1)-liftu(l+1)
-          blsv(iim+1,j)=phyv(l+1)-dragv(l+1)-liftv(l+1)
+          blsu(iim+1,j)=clu(l+1)
+          blsv(iim+1,j)=clv(l+1)
       zlon(iim+1)=-plon(l+1)*xpi/180.
       zlat(j)=plat(l+1)*xpi/180.
@@ -182 +182 @@
           ssou(i,j)=dragu(l)+liftu(l)
           ssov(i,j)=dragv(l)+liftv(l)
-          blsu(i,j)=phyu(l)-dragu(l)-liftu(l)
-          blsv(i,j)=phyv(l)-dragv(l)-liftv(l)
+          blsu(i,j)=clu(l)
+          blsv(i,j)=clv(l)
       zlon(i)=plon(l)*xpi/180.
 
@@ -214 +214 @@
           ssou(i,jjm+1)=dragu(l)+liftu(l)
           ssov(i,jjm+1)=dragv(l)+liftv(l)
-          blsu(i,jjm+1)=phyu(l)-dragu(l)-liftu(l)
-          blsv(i,jjm+1)=phyv(l)-dragv(l)-liftv(l)
+          blsu(i,jjm+1)=clu(l)
+          blsv(i,jjm+1)=clv(l)
       ub(i,jjm+1)=ub(1,jjm+1)                               
       vb(i,jjm+1)=vb(1,jjm+1)                                
@@ -333 +333 @@
      c      tmou(3)/hadley,tsso(3)/hadley,tbls(3)/hadley 
 100   format(F12.5,15(1x,F12.5))
-c100   format(F12.5,5(1x,F12.5))
+c00   format(F12.5,5(1x,F12.5))
 
       write(iam+1,*)((zs(i,j),i=1,iim),j=1,jjm+1)