Changeset 1532 for trunk/LMDZ.MARS/libf/phymars/iniwritesoil.F90

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/iniwritesoil.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/iniwritesoil.F90

@@ -1 @@
-subroutine iniwritesoil(nid,ngrid,inertia,area)
+subroutine iniwritesoil(nid,ngrid,inertia,area,nbplon,nbplat)
 
 ! initialization routine for 'writediagoil'. Here we create/define
@@ -17 +17 @@
 integer,intent(in) :: ngrid
 integer,intent(in) :: nid ! NetCDF output file ID
-real,intent(in) :: inertia(nbp_lon+1,nbp_lat,nsoilmx)
-real,intent(in) :: area(nbp_lon+1,nbp_lat) ! mesh area (m2)
+real,intent(in) :: inertia(nbplon,nbplat,nsoilmx)
+real,intent(in) :: area(nbplon,nbp_lat) ! mesh area (m2)
+integer,intent(in) :: nbplon,nbplat ! sizes of area
 
 ! Local variables:
@@ -33 +34 @@
 real,dimension(nbp_lon+1,nbp_lat,nsoilmx) :: data3 ! to store 3D data
 integer :: i,j,l,ig0
-real :: lon_reg_ext(nbp_lon+1) ! extended longitudes
+real,allocatable :: lon_reg_ext(:) ! extended longitudes
+
+
+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
 
 ! 1. Define the dimensions
@@ -40 +50 @@
 
 ! Define the dimensions
-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
 if (ierr.ne.NF_NOERR) then
   write(*,*)"iniwritesoil: Error, could not define longitude dimension"
@@ -82 +96 @@
 
 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 only:
+  ! add extra redundant point (180 degrees, since lon_reg starts at -180
+  lon_reg_ext(nbp_lon+1)=-lon_reg_ext(1)
+endif
 
 ! Write longitude to file