Changeset 826 for trunk/LMDZ.TITAN/Tools

Timestamp:

Oct 30, 2012, 12:18:36 PM (12 years ago)

Author:

slebonnois

Message:

SL: amelioration angmom.F90 (Venus, Titan Tools) + correction bug GW non-oro (Venus)

File:

• trunk/LMDZ.TITAN/Tools/angmom.F90 (modified) (11 diffs)

trunk/LMDZ.TITAN/Tools/angmom.F90

@@ -16 +16 @@
 !                                       or if simple friction
 ! tdyn  -- 1D -- Dynamics torque IF dudyn is present 
+! tajs  -- 1D -- Torque from convective adjustment IF dudyn is present 
 ! tgwo  -- 1D -- Orographic GW torque IF dugwo is present 
 ! tgwno -- 1D -- Non-Orographic GW torque IF dugwno is present 
@@ -23 +24 @@
 ! - surface pressure and surface geopotential
 ! - zonal wind
-! Optional: dudyn, duvdf, dugwo, dugwno (acceleration terms from physiq param)
+! Optional: dudyn, duvdf, duajs, dugwo, dugwno (acceleration terms from physiq param)
 
 implicit none
@@ -67 +68 @@
 real,dimension(:,:,:,:),allocatable :: duvdf  ! Friction in BL
 real,dimension(:,:,:,:),allocatable :: dudyn  ! Dynamics 
+real,dimension(:,:,:,:),allocatable :: duajs  ! Convective adjustment
 real,dimension(:,:,:,:),allocatable :: dugwo  ! Orographic Gravity Waves
 real,dimension(:,:,:,:),allocatable :: dugwno ! Non-Orographic Gravity Waves
@@ -84 +86 @@
 real,dimension(:),allocatable :: tmou ! mountain torque (kg m2/s2)
 real,dimension(:),allocatable :: tdyn ! dynamics torque (kg m2/s2)
+real,dimension(:),allocatable :: tajs ! convective adjustment torque (kg m2/s2)
 real,dimension(:),allocatable :: tbls ! friction torque (kg m2/s2)
 real,dimension(:),allocatable :: tgwo ! oro GW torque (kg m2/s2)
@@ -99 +102 @@
 integer i,j,ilon,ilat,ilev,itim ! for loops
 integer idlsurf ! for option ideal surface
-logical :: flag_duvdf,flag_dudyn,flag_dugwo,flag_dugwno,lmdflag,dzflag
+logical :: flag_duvdf,flag_dudyn,flag_duajs,flag_dugwo,flag_dugwno,lmdflag,dzflag
 
 real :: deltalat,deltalon ! lat and lon intervals in radians
@@ -236 +239 @@
 allocate(phis(lonlength,latlength))
 
-if(lmdflag) then
-text="psol"
-else
 text="PS"
-endif
 call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
 if (ierr1.ne.NF_NOERR) then
-  write(*,*) "  looking for ps instead... "
-  text="ps"
+  write(*,*) "  looking for psol instead... "
+  text="psol"
   call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
-  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get ps ID"
+  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get psol ID"
 endif
 if (ierr2.ne.NF_NOERR) stop "Error: Failed reading surface pressure"
@@ -374 +373 @@
 else   ! lmdflag
   print*,"dugwo and dugwno not in CAM simulations"
+  flag_dugwo = .false.
+  flag_dugwno = .false.
+endif  ! lmdflag
+
+!===============================================================================
+! 2.1.65 Accelerations from convective adjustment
+!===============================================================================
+allocate(duajs(lonlength,latlength,altlength,timelength))
+
+if(lmdflag) then
+
+text="duajs"
+call get_var4d(infid,lonlength,latlength,altlength,timelength,text,duajs,ierr1,ierr2)
+if (ierr1.ne.NF_NOERR) then
+  write(*,*) "Failed to get duajs ID"
+  flag_duajs = .false.
+else
+  if (ierr2.ne.NF_NOERR) stop "Error: Failed reading duajs"
+  flag_duajs = .true.
+endif
+
+else   ! lmdflag
+  print*,"duajs not in CAM simulations"
+  flag_duajs = .false.
 endif  ! lmdflag
 
@@ -518 +541 @@
 
 !===============================================================================
-! 2.2.5 Mountain, friction, dynamics and GW torques
+! 2.2.5 Mountain, friction, convective adjustment, dynamics and GW torques
 !===============================================================================
 allocate(tmou(timelength))
 if (flag_dudyn)  allocate(tdyn(timelength))
+if (flag_duajs)  allocate(tajs(timelength))
 if (flag_duvdf)  allocate(tbls(timelength))
 if (flag_dugwo)  allocate(tgwo(timelength))
@@ -559 +583 @@
 enddo
   if (tdyn(itim).eq.0.) tdyn(itim) = miss_val
+endif
+
+if (flag_duajs) then
+tajs(itim) = 0.
+do ilon=1,lonlength
+ do ilat=1,latlength
+  do ilev=1,altlength
+    if (rayon(ilon,ilat,ilev,itim).ne.miss_val) then
+      tajs(itim) = tajs(itim) + duajs(ilon,ilat,ilev,itim) &
+       * rayon(ilon,ilat,ilev,itim)*cos(latrad(ilat))      &
+       * dmass(ilon,ilat,ilev,itim) &
+       / hadley
+    endif
+  enddo
+ enddo
+enddo
+  if (tajs(itim).eq.0.) tajs(itim) = miss_val
 endif
 
@@ -682 +723 @@
 !===============================================================================
 
+! Check variables to output
+
+do itim=1,timelength
+  if (flag_dudyn .and.( tdyn(itim).eq.miss_val)) flag_dudyn =.false.
+  if (flag_duajs .and.( tajs(itim).eq.miss_val)) flag_duajs =.false.
+  if (flag_duvdf .and.( tvdf(itim).eq.miss_val)) flag_duvdf =.false.
+  if (flag_dugwo .and.( tgwo(itim).eq.miss_val)) flag_dugwo =.false.
+  if (flag_dugwno.and.(tgwno(itim).eq.miss_val)) flag_dugwno=.false.
+enddo ! timelength
+
+if(dzflag) then
+do itim=1,timelength
+  if (tdyndz(itim).eq.miss_val) dzflag=.false.
+  if (tdisdz(itim).eq.miss_val) dzflag=.false.
+  if (tspgdz(itim).eq.miss_val) dzflag=.false.
+  if (tphydz(itim).eq.miss_val) dzflag=.false.
+enddo ! timelength
+endif ! dzflag
+
+
 ! 1D Variables
 
@@ -702 +763 @@
                 "tdyn      ", "Dynamics torque     ","E18kgm2s-2",miss_val,&
                  tdyn )
+endif
+
+if (flag_duajs) then
+call write_var1d(outfid,datashape1d,timelength,&
+                "tajs      ", "Dynamics torque     ","E18kgm2s-2",miss_val,&
+                 tajs )
 endif