Changeset 1532 for trunk/LMDZ.MARS/libf/phymars/iniwrite.F

Timestamp:

Apr 7, 2016, 3:53:15 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• Some fixes for buggy outputs in 1D introduced by previous code modifications.
• made statto.h a module.
• ecritphy in dyn3d/control_mod.F90 should be an integer, not a real.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/iniwrite.F (modified) (6 diffs)

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

@@ -1 @@
-      SUBROUTINE iniwrite(nid,idayref,phis,area)
+      SUBROUTINE iniwrite(nid,idayref,phis,area,nbplon,nbplat)
 
       use comsoil_h, only: mlayer, nsoilmx
@@ -35 +35 @@
       integer,intent(in) :: nid        ! NetCDF file ID
       INTEGER*4,intent(in) :: idayref  ! date (initial date for this run)
-      real,intent(in) :: phis(nbp_lon+1,nbp_lat) ! surface geopotential
-      real,intent(in) :: area(nbp_lon+1,nbp_lat) ! mesh area (m2)
+      real,intent(in) :: phis(nbplon,nbp_lat) ! surface geopotential
+      real,intent(in) :: area(nbplon,nbp_lat) ! mesh area (m2)
+      integer,intent(in) :: nbplon,nbplat ! sizes of area and phis arrays
 
 c   Local:
@@ -44 +45 @@
       REAL tab_cntrl(length) ! run parameters are stored in this array
       INTEGER ierr
-      REAl :: lon_reg_ext(nbp_lon+1) ! extended longitudes
+      REAl,ALLOCATABLE :: lon_reg_ext(:) ! extended longitudes
 
       integer :: nvarid,idim_index,idim_rlonv
@@ -51 +52 @@
       integer, dimension(2) :: id  
 c-----------------------------------------------------------------------
+
+      IF (nbp_lon*nbp_lat==1) THEN
+        ! 1D model
+        ALLOCATE(lon_reg_ext(1))
+      ELSE
+        ! 3D model
+        ALLOCATE(lon_reg_ext(nbp_lon+1))
+      ENDIF
 
       DO l=1,length
@@ -104 +113 @@
 !      ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
       ierr = NF_DEF_DIM (nid, "latitude", nbp_lat, idim_rlatu)
-      ierr = NF_DEF_DIM (nid, "longitude", nbp_lon+1, idim_rlonv)
+      IF (nbp_lon*nbp_lat==1) THEN
+        ierr = NF_DEF_DIM (nid, "longitude", 1, idim_rlonv)
+      ELSE
+        ierr = NF_DEF_DIM (nid, "longitude", nbp_lon+1, idim_rlonv)
+      ENDIF
 !      ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
       ierr = NF_DEF_DIM (nid, "interlayer", (nbp_lev+1), idim_llmp1)
@@ -166 +179 @@
 c
 c --------------------------
+      
       lon_reg_ext(1:nbp_lon)=lon_reg(1:nbp_lon)
-      !add extra redundant point (180 degrees, since lon_reg starts at -180
-      lon_reg_ext(nbp_lon+1)=-lon_reg_ext(1)
+      IF (nbp_lon*nbp_lat/=1) THEN
+        ! In 3D, add extra redundant point (180 degrees,
+        ! since lon_reg starts at -180)
+        lon_reg_ext(nbp_lon+1)=-lon_reg_ext(1)
+      ENDIF
 
       ierr = NF_REDEF (nid)