Changeset 1497 for trunk/MESOSCALE

Timestamp:

Nov 12, 2015, 3:37:25 PM (9 years ago)

Author:

aslmd

Message:

added commented lines for hill LES case

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars

Files:

• module_initialize_les.F.hill (deleted)
• module_initialize_les.F.mars (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F.mars

@@ -94 +94 @@
    REAL    :: z_scale, xrad, yrad, zrad, rad, delt, cof1, cof2
 !   REAL, EXTERNAL :: interp_0
-   REAL    :: hm
+   REAL    :: hm, xa
    REAL    :: pi
 
@@ -331 +331 @@
   grid%p_top = interp_0( p_in, zk, config_flags%ztop + alt_input, nl_in )
 
+
+!!! add hill
+!             icm = ide/2
+!             jcm = jde/2
+!             !!3D hill 
+!             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!             !!MARS : mountain
+!             !!MARS : mountain ex. hm = 2000. xa = 6.0
+!               open(unit=22,file='ze_hill',form='formatted',status='old')
+!               rewind(22)
+!               read(22,*) hm, xa
+!               write(6,*) 'height, width ', hm, xa
+!               close(22)
+!             !!MARS
+!             !!MARS
+!             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+
+
   DO j=jts,jte
   DO i=its,ite
 !!MARS
     grid%ht(i,j) = alt_input
+
+!!! add hill
+!             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!             grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
+!             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
     grid%m_tsurf(i,j) = tsurf_input
 !!MARS