Changeset 999 for trunk/LMDZ.MARS/libf/phymars/soil_settings.F

Timestamp:

Jul 2, 2013, 9:40:28 AM (11 years ago)

Author:

tnavarro

Message:

Possibility to store multiple initial states in one start/startfi. This is RETROCOMPATIBLE. New option ecrithist in run.def to write data in start/startfi every ecrithist dynamical timestep. New option timestart in run.def to initialize the GCM with the time timestart stored in start

File:

• trunk/LMDZ.MARS/libf/phymars/soil_settings.F (modified) (5 diffs)

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

@@ -1 @@
-      subroutine soil_settings(nid,ngrid,nsoil,tsurf,tsoil)
+      subroutine soil_settings(nid,ngrid,nsoil,tsurf,tsoil,indextime)
 
       use netcdf
@@ -12 +12 @@
 !                      r4 or r8 restarts independently of having compiled
 !                      the GCM in r4 or r8)
+!                June 2013 TN : Possibility to read files with a time axis
+!
 !
 !  This subroutine reads from a NetCDF file (opened by the caller)
@@ -42 +44 @@
       integer nsoil     ! # of soil layers
       real tsurf(ngrid) ! surface temperature
+      integer indextime ! position on time axis
 !  output:
       real tsoil(ngridmx,nsoilmx)       ! soil temperature
@@ -53 +56 @@
       integer ndims     ! # of dimensions of read <inertiedat> data
       integer ig,iloop  ! loop counters
+      
+      integer edges(3),corner(3) ! to read a specific time
 
       logical :: olddepthdef=.false. ! flag
@@ -298 +303 @@
         endif
        else ! put values in tsoil
+        corner(1)=1
+        corner(2)=1
+        corner(3)=indextime
+        edges(1)=ngridmx
+        edges(2)=nsoilmx
+        edges(3)=1
+        !ierr=nf90_get_var(nid,nvarid,tsoil,corner,edges)
         ierr=nf90_get_var(nid,nvarid,tsoil)
         if (ierr.ne.nf90_noerr) then