Changeset 819

Timestamp:

Oct 26, 2012, 11:31:16 AM (12 years ago)

Author:

slebonnois

Message:

• newstart-VT.F and start2archive-VT.F added in dyn3d/, their name indicating that they are tools for Venus and Titan.
• doc for these tools added
• correction of a bug on angmom (Venus and Titan postprocessing tools)
• one tool (lat_torques) removed because not updated until further needed
• correction of a bug in phyetat0 (Venus)

Location:

trunk

Files:

• DOC/documentation/newstart.pdf (added)
• DOC/documentation/newstart.tex (added)
• LMDZ.COMMON/libf/dyn3d/newstart-VT.F (added)
• LMDZ.COMMON/libf/dyn3d/start2archive-VT.F (added)
• LMDZ.TITAN/Tools/README (modified) (1 diff)
• LMDZ.TITAN/Tools/angmom.F90 (modified) (1 diff)
• LMDZ.TITAN/Tools/lat_torques.F90 (deleted)
• LMDZ.VENUS/Tools/README (modified) (1 diff)
• LMDZ.VENUS/Tools/angmom.F90 (modified) (1 diff)
• LMDZ.VENUS/Tools/lat_torques.F90 (deleted)
• LMDZ.VENUS/libf/phyvenus/phyetat0.F (modified) (1 diff)

trunk/LMDZ.TITAN/Tools/README

@@ -2 +2 @@
 DOC/documentation/vt-tools.pdf
 
+List of tools:
+zrecast
+angmom
+energy
+fft 
+psi
+stability
+tem
+tmc
+
 Known bugs:
-- when using a file histmth_P.nc, missing surface pressure...
-I need to add it in zrecast.
 

trunk/LMDZ.TITAN/Tools/angmom.F90

@@ -236 +236 @@
 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 psol instead... "
-  text="psol"
+  write(*,*) "  looking for ps instead... "
+  text="ps"
   call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
-  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get psol ID"
+  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get ps ID"
 endif
 if (ierr2.ne.NF_NOERR) stop "Error: Failed reading surface pressure"

trunk/LMDZ.VENUS/Tools/README

@@ -2 +2 @@
 DOC/documentation/vt-tools.pdf
 
+List of tools:
+zrecast
+angmom
+energy
+fft 
+psi
+stability
+tem
+tmc
+
 Known bugs:
-- when using a file histmth_P.nc, missing surface pressure...
-I need to add it in zrecast.
 

trunk/LMDZ.VENUS/Tools/angmom.F90

@@ -236 +236 @@
 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 psol instead... "
-  text="psol"
+  write(*,*) "  looking for ps instead... "
+  text="ps"
   call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
-  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get psol ID"
+  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get ps ID"
 endif
 if (ierr2.ne.NF_NOERR) stop "Error: Failed reading surface pressure"

trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F

@@ -9 +9 @@
      .           fder,radsol,
      .    zmea, zstd, zsig, zgam, zthe, zpic, zval,  
-     .           tabcntr0,
      .           t_ancien,ancien_ok)