Changeset 317

Timestamp:

Oct 11, 2011, 4:17:58 PM (14 years ago)

Author:

aslmd

Message:

GRAPHICS: - winds.py works for diagfi.nc

• still a problem with ortho projection
• better handling for time should be done

example:
winds.py -f diagfi1.nc -x -v fluxsurf_lw -n 1 -F 450

LMD_MM_MARS: cosmetic changes. checked that z0 changes are OK with old input files.

Location:

trunk

Files:

• LMDZ.MARS/libf/phymars/surfini.F (modified) (1 diff)
• MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py (modified) (6 diffs)
• MESOSCALE/LMD_MM_MARS/SRC/PYTHON/winds.py (modified) (1 diff)
• MESOSCALE_DEV/MANUAL/SRC/advance.tex (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/surfini.F

@@ -50 +50 @@
         
 #ifdef MESOSCALE
+         !write(*,*) "all qsurf to zero. dirty."
          do iq=1,nqmx
-         write(*,*) "all qsurf to zero. dirty."
          qsurf(ig,iq)=0.  !! on jette les inputs GCM
                           !! on regle juste watercaptag

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py

@@ -21 +21 @@
 def whatkindfile (nc):
     if 'controle' in nc.variables:   typefile = 'gcm'
-    elif 'phisinit' in nc.variables: typefile = 'gcmex'
+    elif 'phisinit' in nc.variables: typefile = 'gcm' 
     elif 'vert' in nc.variables:     typefile = 'mesoapi'
     elif 'U' in nc.variables:        typefile = 'meso'
@@ -152 +152 @@
     if zetime is not None \
        and 'vert' not in nc.variables:
-       #### strangely enough this does not work for api or ncrcat results!
+        #### strangely enough this does not work for api or ncrcat results!
         zetimestart = getattr(nc, 'START_DATE')
         zeday = int(zetime[8]+zetime[9]) - int(zetimestart[8]+zetimestart[9])
@@ -254 +254 @@
 
 def getcoorddef ( nc ):   
+    import numpy as np
     ## getcoord2d for predefined types
     typefile = whatkindfile(nc)
@@ -260 +261 @@
         lon2d = dumpbdy(lon2d,6)
         lat2d = dumpbdy(lat2d,6)
-    elif typefile in ['gcm','gcmex']:
+    elif typefile in ['gcm']: 
         [lon2d,lat2d] = getcoord2d(nc,nlat="latitude",nlon="longitude",is1d=True)
     elif typefile in ['geo']:
@@ -495 +496 @@
              "Um":           "%.0f",\
              "ALBBARE":      "%.2f",\
+             "TAU":          "%.1f",\
                     }
     if whichvar not in fmtvar:
@@ -520 +522 @@
              "QSURFICE":     "hot_r",\
              "ALBBARE":      "spectral",\
+             "TAU":          "YlOrBr_r",\
                      }
 #W --> spectral ou jet

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/winds.py

@@ -196 +196 @@
     ### Save the figure in a file in the data folder or an user-defined folder
     if typefile in ['meso','mesoapi']:   prefix = getprefix(nc)   
-    elif typefile in ['gcm','gcmex']:    prefix = 'LMD_GCM_'
+    elif typefile in ['gcm']:            prefix = 'LMD_GCM_'
     else:                                prefix = ''
     ###

trunk/MESOSCALE_DEV/MANUAL/SRC/advance.tex

@@ -130 +130 @@
 %%[set to 18 for newphys]
 
-%% attention init_TI n'est plus sÃparable... il faut run real.exe a cause du modele sous surface
+%% attention init_TI n'est plus separable... il faut run real.exe a cause du modele sous surface
 %% il faut ajouter init_Z0 ! ou il n'y a pas besoin de re executer real.exe
 
+%% si z0 n est pas dans le wrfinput on prend la valeur par defaut
+
 \clearemptydoublepage