Changeset 1768 for trunk

Timestamp:

Aug 25, 2017, 2:23:32 AM (7 years ago)

Author:

aslmd

Message:

LES. complement to previous commit. make read levels as an alternative for stretch_grid for non-mars cases (as it used to be before the previous commit)

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F (modified) (2 diffs)

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

@@ -253 +253 @@
                                 (1.-exp(-1./z_scale))
      ENDDO
+   ELSE
+     !!MARS
+     !!MARS
+     IF (planet .ne. "mars") THEN
+       open(unit=12,file='levels',form='formatted',status='old')
+       rewind(12)
+       DO k=1, kde
+        read(12,*) grid%znw(k)
+        write(6,*) 'read level ', k,grid%znw(k)
+       ENDDO
+       close(12)
+     ENDIF
+     !!MARS
+     !!MARS
+     !  !DO k=1, kde
+     !  ! grid%znw(k) = 1. - float(k-1)/float(kde-1)
+     !  !ENDDO
    ENDIF
-   !ELSE
-   !  !DO k=1, kde
-   !  ! grid%znw(k) = 1. - float(k-1)/float(kde-1)
-   !  !ENDDO
-
-
+
+   !! SPECIFIC FOR LES PBL MARS
    IF (planet == "mars") THEN
      !!!MARS
@@ -269 +282 @@
        grid%znw(k) = grid%znw(4) * ( 1. - float(k-4)/float(kde-4) )
      ENDDO
-!!!MARS
-!!!MARS
-!  open(unit=12,file='levels',form='formatted',status='old')
-!  rewind(12)
-!  DO k=1, kde
-!  read(12,*) grid%znw(k)
-!  write(6,*) 'read level ', k,grid%znw(k)
-!  ENDDO
-!  close(12)
-!!!MARS
-!!!MARS
    ENDIF