Changeset 2147

Timestamp:

Jul 16, 2019, 12:36:42 PM (5 years ago)

Author:

emillour

Message:

Mars GCM utilities:
Add the information in concat when time axis is set to "Solar longitude", and try to check in localtime that the input time axis in indeed in sols.
EM

Location:

trunk/LMDZ.MARS/util

Files:

• concatnc.F90 (modified) (4 diffs)
• concatnc.def (modified) (1 diff)
• localtime.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/util/concatnc.F90

@@ -404 +404 @@
       call initiate (filename,lat,lon,alt,ctl,GCM_layers,nout,&
        latdimout,londimout,altdimout,timedimout,&
-       layerdimout,interlayerdimout,timevarout)
+       layerdimout,interlayerdimout,timevarout,axis)
    ! Initialize output file's aps,bps,ap,bp and phisinit variables
      call init2(nid,lonlen,latlen,altlen,GCM_layers,&
@@ -619 +619 @@
 subroutine initiate (filename,lat,lon,alt,ctl,GCM_layers,nout,&
          latdimout,londimout,altdimout,timedimout,&
-         layerdimout,interlayerdimout,timevarout)
+         layerdimout,interlayerdimout,timevarout,axis)
 !==============================================================================
 ! Purpose:
@@ -662 +662 @@
 integer, intent(out):: timevarout
 ! timevarout: [netcdf] "Time" (considered as a variable) ID
+character (len=4),intent(in) :: axis
+! axis: "ls" or "sol"
 
 !==============================================================================
@@ -701 +703 @@
 ! 3. Write "Time" and its attributes
 !==============================================================================
-
-call def_var(nout,"Time","Time","years since 0000-00-0 00:00:00",1,&
+if (axis=="sol") then
+  call def_var(nout,"Time","Time","days since 0000-00-0 00:00:00",1,&
              (/timedimout/),timevarout,ierr)
-
+else ! Ls
+  call def_var(nout,"Time","Solar longitude","days since 0000-00-0 00:00:00",1,&
+             (/timedimout/),timevarout,ierr)
+endif
 !==============================================================================
 ! 4. Write "latitude" (data and attributes)

trunk/LMDZ.MARS/util/concatnc.def

@@ -4 +4 @@
 
 515              ! Starting sol of the run stored in the first input file?
-sol              ! output timescale ("sol" or "Ls")
+sol              ! output timescale ("sol" or "ls")
 tsurf            
 ps

trunk/LMDZ.MARS/util/localtime.F90

@@ -302 +302 @@
 !  write(*,*) "timelen: ",timelen
 
+   ! Sanity Check: Does "Time" has a "title" attribute
+   ! and is it "Solar longitude" ?
+   ierr=nf_get_att_text(nid,timevar,"title",title)
+   if ((ierr.EQ.NF_NOERR).and.(title.eq."Solar longitude")) then
+     write(*,*) "ERROR: Time axis in input file is in Solar Longitude!"
+     write(*,*) "       localtime requires sols as time axis!"
+     write(*,*) " Might as well stop here."
+     stop
+   endif
+
    ! allocate time() array and fill it with values from input file
    allocate(time(timelen))
@@ -310 +320 @@
    ierr = NF_GET_VAR_REAL(nid,timevar,time)
 #endif
+   if (ierr.NE.NF_NOERR) then
+     write(*,*) "Error , failed to load Time"
+     stop
+   endif
 
 
@@ -396 +410 @@
         ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,lt_out(it))
 #endif
+        if (ierr.NE.NF_NOERR) then
+          write(*,*) "Error , failed to write Time"
+          stop
+        endif
      enddo
   end if