Changeset 322 for trunk/LMDZ.MARS/util

Timestamp:

Oct 21, 2011, 10:59:42 AM (13 years ago)

Author:

emillour

Message:

Mars Utility:
Bug fix in concatnc by TN: outputs in Ls would sometimes be wrong by
+ or - 360 at transitions between years.
EM

File:

• trunk/LMDZ.MARS/util/concatnc.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/util/concatnc.F90

@@ -1085 +1085 @@
 ierr=NF_INQ_VARID(nid,"Time",nvarid)
 if (ierr.NE.NF_NOERR) then
-   write(*,*) 'ERROR: Field <Time> not found'
+   write(*,*) 'ERROR in change_time_axis: Field <Time> not found'
    print*, NF_STRERROR(ierr)
    stop ""
@@ -1097 +1097 @@
 ierr = NF_GET_VAR_REAL(nid,nvarid,time)
 #endif
+if (ierr.ne.NF_NOERR) then
+  write(*,*) "ERROR in change_time_axis: Failed to load Time"
+  stop
+endif
 
 !==============================================================================
@@ -1107 +1111 @@
 
 !==============================================================================
-! 2. Write
+! 2.1 Check if there are not jumps in Ls (rounding problems in sol2ls)
+!==============================================================================
+
+do i=1,timelen-1
+   if ((ls(i+1)-ls(i)) > 350) then
+       write(*,*) "+ 360° Ls jump solved:", ls(i), ls(i+1), "at timestep", i
+      ls(i+1) = ls(i+1) - 360
+       write(*,*) " corrected to now be:", ls(i), ls(i+1)
+   else if ((ls(i)-ls(i+1)) > 350) then
+       write(*,*) "- 360° Ls jump solved:", ls(i), ls(i+1), "at timestep", i
+      ls(i+1) = ls(i+1) + 360
+       write(*,*) " corrected to now be:", ls(i), ls(i+1)
+   endif
+enddo
+
+!==============================================================================
+! 3. Write
 !==============================================================================
 
@@ -1115 +1135 @@
 ierr = NF_PUT_VAR_REAL(nid,nvarid,ls)
 #endif
+if (ierr.ne.NF_NOERR) then
+  write(*,*) "ERROR in change_time_axis: Failed to write Ls"
+  stop
+endif
 
 end subroutine change_time_axis