Changeset 143 for trunk/MESOSCALE/LMD_MM_MARS

Timestamp:

Jun 1, 2011, 12:25:21 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: modifications dans module_initialize_quarter_ss.F pour que la presence/absence de ze_hill et son contenu controle la topographie dans les simulations idealisees (plate, pente +/-, montagne, cratere) sans qu'il ne soit donc desormais necessaire de recompiler

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F (modified) (4 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F

@@ -127 +127 @@
 !!MARS
  REAL :: lon_input, lat_input, alt_input, tsurf_input
+ INTEGER :: ierr
 !!MARS
 
@@ -176 +177 @@
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!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 : NOFILE      === no mountain
+!!MARS : FILE xa=0.  === a linear slope with elevation hm
+!!MARS : FILE        === mountain height hm, width xa
+  open(unit=22,file='ze_hill',form='formatted',status='old',iostat=ierr)
+  IF (ierr .eq. 0) THEN
+    rewind(22) 
+    read(22,*) hm, xa
+    write(6,*) 'hm, xa ', hm, xa
+    close(22)
+  ENDIF
 !!MARS
 !!MARS
@@ -390 +394 @@
   DO i=its,ite
 !!MARS
-    grid%ht(i,j) = alt_input
-      ! flat surface
-      !!    grid%ht(i,j) = 0.
+IF (ierr .eq. 0) THEN
+  write(6,*) ' IDEALIZED TOPOGRAPHY '
+  IF (xa .ne. 0.) THEN
+      !!!2D hill
+      !grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2)
       !    grid%ht(i,j) = hm1*exp(-(( float(i-icm)/xa1)**2))   &
       !               *( (cos(pii*float(i-icm)/xal1))**2 )
-      !****Mars
-      !!3D hill
-      !        grid%ht(i,j) = hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
-
-!!!2D hill
-!grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2)
-
-!!3D hill
-grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
-
-!!!simple slope
-!grid%ht(i,j) = alt_input - hm * float(i)  
-
-
+    write(6,*) '3D hill or crater '
+    write(6,*) 'height, width ',hm,xa
+    grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
+  ELSE
+    write(6,*) 'linear slope '
+    write(6,*) 'height ',hm
+    grid%ht(i,j) = alt_input + hm * float(i)  
+  ENDIF
       !!!3D crater    
       !!        grid%ht(i,j) = hm - hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
@@ -425 +425 @@
       !           !grid%ht(i,j) = grid%ht(i,j) + alt_input  
       !   !if (rho_param .GE. dilat*10) ht(i,j) = hm
-
+ELSE
+    write(6,*) ' FLAT SURFACE '
+    grid%ht(i,j) = alt_input
+    !grid%ht(i,j) = 0.
+ENDIF
     grid%tsk(i,j) = tsurf_input
 !!MARS