Changeset 2567

Timestamp:

Oct 5, 2021, 8:34:09 AM (3 years ago)

Author:

emillour

Message:

Mars GCM utilities:
Minor fixes to run with picky gfortran 10.3.0 which requires one element arrays
(rather than scalars) when calling NetCDF routines, andf that stop statements
should not be followed by strings. While at it replaced tabs with spaces.
EM

Location:

trunk/LMDZ.MARS/util

Files:

• aeroptical.F90 (modified) (21 diffs)
• concatnc.F90 (modified) (45 diffs)
• hrecast.F90 (modified) (41 diffs)
• localtime.F90 (modified) (41 diffs)
• lslin.F90 (modified) (34 diffs)
• solzenangle.F90 (modified) (41 diffs)
• streamfunction.F90 (modified) (37 diffs)
• xvik/fit_Iceinertia_MONSicedepth.F (modified) (14 diffs)
• zrecast.F90 (modified) (53 diffs)

trunk/LMDZ.MARS/util/aeroptical.F90

@@ -1304 +1304 @@
   hybrid=.true.
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed reading aps"
+    write(*,*) "init2 Error: Failed reading aps"
+    stop
   endif
 
@@ -1316 +1317 @@
   ierr=nf90_inq_varid(gcmfid,"bps",tmpvarid)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to get bps ID."
+    write(*,*) "init2 Error: Failed to get bps ID."
+    stop
   endif
   ierr=NF90_GET_VAR(gcmfid,tmpvarid,bps)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed reading bps"
+    write(*,*) "init2 Error: Failed reading bps"
+    stop
   endif
 endif
@@ -1338 +1341 @@
     apbp=.true.
     if (ierr.ne.nf90_noerr) then
-      stop "Error: Failed reading ap"
+      write(*,*) "Error: Failed reading ap"
+      stop
     endif
   endif
@@ -1357 +1361 @@
     apbp=.true.
     if (ierr.ne.nf90_noerr) then
-      stop "Error: Failed reading bp"
+      write(*,*) "Error: Failed reading bp"
+      stop
     endif
   endif
@@ -1372 +1377 @@
   ierr=NF90_GET_VAR(gcmfid,tmpvarid,sigma)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed reading sigma"
+    write(*,*) "init2 Error: Failed reading sigma"
+    stop
   endif
 endif ! of if (.not.hybrid)
@@ -1389 +1395 @@
   ierr=NF90_GET_VAR(gcmfid,tmpvarid,aire)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed reading aire"
+    write(*,*) "init2 Error: Failed reading aire"
+    stop
   endif
   area = .true.
@@ -1407 +1414 @@
   ierr=NF90_GET_VAR(gcmfid,tmpvarid,phisinit)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -1426 +1434 @@
   ierr=NF90_DEF_VAR(outfid,"aps",nf90_float,(/layerdimout/),tmpvarid)
   if (ierr.ne.nf90_noerr) then
-     stop "init2 Error: Failed to define the variable aps"
+     write(*,*) "init2 Error: Failed to define the variable aps"
+     stop
   endif
   ! Write the attributes
@@ -1437 +1446 @@
   ierr=NF90_PUT_VAR(outfid,tmpvarid,aps)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to write aps"
+    write(*,*) "init2 Error: Failed to write aps"
+    stop
   endif
 
@@ -1446 +1456 @@
   ierr=NF90_DEF_VAR(outfid,"bps",nf90_float,(/layerdimout/),tmpvarid)
   if (ierr.ne.nf90_noerr) then
-     stop "init2 Error: Failed to define the variable bps"
+     write(*,*) "init2 Error: Failed to define the variable bps"
+     stop
   endif
   ! Write the attributes
@@ -1457 +1468 @@
   ierr=NF90_PUT_VAR(outfid,tmpvarid,bps)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to write bps"
+    write(*,*) "init2 Error: Failed to write bps"
+    stop
   endif
 
@@ -1468 +1480 @@
     ierr=NF90_DEF_VAR(outfid,"ap",nf90_float,(/interlayerdimout/),tmpvarid)
     if (ierr.ne.nf90_noerr) then
-       stop "init2 Error: Failed to define the variable ap"
+       write(*,*) "init2 Error: Failed to define the variable ap"
+       stop
     endif
     ! Write the attributes
@@ -1479 +1492 @@
     ierr=NF90_PUT_VAR(outfid,tmpvarid,ap)
     if (ierr.ne.nf90_noerr) then
-      stop "Error: Failed to write ap"
+      write(*,*) "Error: Failed to write ap"
+      stop
     endif
 
@@ -1489 +1503 @@
     ierr=NF90_DEF_VAR(outfid,"bp",nf90_float,(/interlayerdimout/),tmpvarid)
     if (ierr.ne.nf90_noerr) then
-       stop "init2 Error: Failed to define the variable bp"
+       write(*,*) "init2 Error: Failed to define the variable bp"
+       stop
     endif
     ! Write the attributes
@@ -1500 +1515 @@
     ierr=NF90_PUT_VAR(outfid,tmpvarid,bp)
     if (ierr.ne.nf90_noerr) then
-      stop "Error: Failed to write bp"
+      write(*,*) "Error: Failed to write bp"
+      stop
     endif
   endif ! of if (apbp)
@@ -1511 +1527 @@
   ierr=NF90_DEF_VAR(outfid,"sigma",nf90_float,(/layerdimout/),tmpvarid)
   if (ierr.ne.nf90_noerr) then
-     stop "init2 Error: Failed to define the variable sigma"
+     write(*,*) "init2 Error: Failed to define the variable sigma"
+     stop
   endif
   ! Write the attributes
@@ -1522 +1539 @@
   ierr=NF90_PUT_VAR(outfid,tmpvarid,sigma)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to write sigma"
+    write(*,*) "init2 Error: Failed to write sigma"
+    stop
   endif
 endif ! of if (hybrid)
@@ -1537 +1555 @@
   ierr=NF90_DEF_VAR(outfid,"aire",nf90_float,(/londimout,latdimout/),tmpvarid)
   if (ierr.ne.nf90_noerr) then
-     stop "init2 Error: Failed to define the variable aire"
+     write(*,*) "init2 Error: Failed to define the variable aire"
+     stop
   endif
   ! Write the attributes
@@ -1548 +1567 @@
   ierr=NF90_PUT_VAR(outfid,tmpvarid,aire)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to write aire"
+    write(*,*) "init2 Error: Failed to write aire"
+    stop
   endif
 endif ! of if (area)
@@ -1559 +1579 @@
   ierr=NF90_DEF_VAR(outfid,"phisinit",nf90_float,(/londimout,latdimout/),tmpvarid)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to define the variable phisinit"
+    write(*,*) "init2 Error: Failed to define the variable phisinit"
+    stop
   endif
   ! Write the attributes
@@ -1570 +1591 @@
   ierr=NF90_PUT_VAR(outfid,tmpvarid,phisinit)
   if (ierr.ne.nf90_noerr) then
-    stop "init2 Error: Failed to write phisinit"
+    write(*,*) "init2 Error: Failed to write phisinit"
+    stop
   endif
 

trunk/LMDZ.MARS/util/concatnc.F90

@@ -28 +28 @@
 include "netcdf.inc" ! NetCDF definitions
 
-character (len=80), dimension(1000) :: file
+character (len=100), dimension(1000) :: file
 ! file(): input file(s) names(s)
 character (len=30), dimension(16) :: notconcat
 ! notconcat(): names of the (16) variables that won't be concatenated
-character (len=50), dimension(:), allocatable :: var
+character (len=100), dimension(:), allocatable :: var
 ! var(): name(s) of variable(s) that will be concatenated
 character (len=100) :: tmpvar,tmpfile,long_name,units
@@ -153 +153 @@
 if(nbfile==0) then
    write(*,*) "no file... game over"
-   stop ""
+   stop
 endif
 
@@ -183 +183 @@
 if (ierr.NE.NF_NOERR) then
    write(*,*) 'ERROR: Pb opening file '//file(1)
-   stop ""
+   stop
 endif
 
@@ -250 +250 @@
 write(*,*) "(an empty line , i.e: just <Enter>, implies end of list)"
 nbvar=0
-read(*,'(a50)') tmpvar
+read(*,'(a100)') tmpvar
 do while ((tmpvar/=' ').AND.(trim(tmpvar)/='all'))
    nbvar=nbvar+1
    var(nbvar)=tmpvar
-   read(*,'(a50)') tmpvar
+   read(*,'(a100)') tmpvar
 enddo
 
@@ -287 +287 @@
 else if(nbvar==0) then
    write(*,*) "no variable... game over"
-   stop ""
+   stop
 endif ! of if (tmpvar=="all")
 
@@ -316 +316 @@
          write(*,*) 'ERROR: Pb opening file '//file(i)
          write(*,*) NF_STRERROR(ierr)
-         stop ""
+         stop
       endif
    endif
@@ -338 +338 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <latitude> or <lat> is missing in file '//file(i)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_DIMLEN(nid,latdim,latlen)
@@ -347 +347 @@
        ierr=NF_INQ_VARID(nid,"longitude",lonvar)
        if (ierr.NE.NF_NOERR) then
-          write(*,*) 'ERROR: Field <longitude> is missing in file '//file(i)
-          stop "" 
+          write(*,*) 'ERROR: Field <longitude> is missing in file '//file(i)
+          stop 
        endif
        ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
@@ -357 +357 @@
        if (ierr.NE.NF_NOERR) then
           write(*,*) 'ERROR: Field <lon> is missing in file '//file(i)
-          stop "" 
+          stop 
        endif
        ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
@@ -367 +367 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <altitude> is missing in file '//file(i)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,altdim,altlen)
@@ -389 +389 @@
 !==============================================================================
    if (XIOS.EQ.0) then
-        ierr=NF_INQ_DIMID(nid,"index",ctldim)
-        if (ierr.NE.NF_NOERR) then
-          write(*,*) 'Dimension <index> is missing in file '//file(i)
-          write(*,*) "The program continues..."
-          ctllen=0
-          !stop ""
+        ierr=NF_INQ_DIMID(nid,"index",ctldim)
+        if (ierr.NE.NF_NOERR) then
+          write(*,*) 'Dimension <index> is missing in file '//file(i)
+          write(*,*) "The program continues..."
+          ctllen=0
+          !stop
        else
-          ierr=NF_INQ_VARID(nid,"controle",ctlvar)
-          if (ierr.NE.NF_NOERR) then
-            write(*,*) 'Field <controle> is missing in file '//file(i)
-            write(*,*) "The program continues..."
-            ctllen=0
-            !stop ""
-          else
+          ierr=NF_INQ_VARID(nid,"controle",ctlvar)
+          if (ierr.NE.NF_NOERR) then
+            write(*,*) 'Field <controle> is missing in file '//file(i)
+            write(*,*) "The program continues..."
+            ctllen=0
+            !stop
+          else
             ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
-          endif
+          endif
        endif
    else
@@ -412 +412 @@
        write(*,*) "The program continues..."
        ctllen=0
-       !stop ""
+       !stop
      else
        ierr=NF_INQ_VARID(nid,"controle",ctlvar)
        if (ierr.NE.NF_NOERR) then
-         write(*,*) 'Field <controle> is missing in file '//file(i)
-         write(*,*) "The program continues..."
-         ctllen=0
-         !stop ""
+          write(*,*) 'Field <controle> is missing in file '//file(i)
+           write(*,*) "The program continues..."
+          ctllen=0
+          !stop
         else
-         ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
+          ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
        endif
      endif
@@ -439 +439 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <Time> is missing in file '//file(i)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_VARID(nid,"Time",timevar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <Time> is missing in file '//file(i)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,timedim,timelen)
@@ -452 +452 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <time_counter> is missing in file '//file(i)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_VARID(nid,"time_counter",timevar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <time_counter> is missing in file '//file(i)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,timedim,timelen)
@@ -542 +542 @@
       if (memolatlen/=latlen) then
            write(*,*) "ERROR: Not the same latitude axis"
-           stop ""
+           stop
        else if (memolonlen/=lonlen) then
            write(*,*) "ERROR: Not the same longitude axis"
-           stop ""
+           stop
        else if (memoaltlen/=altlen) then
            write(*,*) "ERROR: Not the same altitude axis"
-           stop ""
+           stop
        else if (memoctllen/=ctllen) then
            write(*,*) "ERROR: Not the same controle axis"
-           stop ""
+           stop
        endif
    endif ! of if (i==1)
@@ -582 +582 @@
        if (rep.eq.1) write(*,*) "Sol", int(output_time)  
 
-       ierr= NF_PUT_VARA_REAL(nout,timevarout,k,1,output_time)
+       ierr= NF_PUT_VARA_REAL(nout,timevarout,(/k/),(/1/),output_time)
    end do
 !  use the last output_time value to update memotime   
@@ -602 +602 @@
       if (ierr.NE.NF_NOERR) then
          write(*,*) 'ERROR: Field <',var(j),'> not found in file '//file(i)
-         stop ""
+         stop
       endif
       ierr=nf_inq_varndims(nid,varid,ndim)
@@ -630 +630 @@
          corner(4)=1
 
-         ! length (along dimensions) of block of data to be written
+         ! length (along dimensions) of block of data to be written
          edges(1)=timelen 
          edges(2)=1 
@@ -648 +648 @@
          corner(4)=1
 
-         ! length (along dimensions) of block of data to be written
+         ! length (along dimensions) of block of data to be written
          edges(1)=lonlen 
          edges(2)=latlen 
@@ -667 +667 @@
          corner(4)=reptime+1
 
-         ! length (along dimensions) of block of data to be written
+         ! length (along dimensions) of block of data to be written
          edges(1)=lonlen
          edges(2)=latlen
@@ -699 +699 @@
       if (ierr.ne.NF_NOERR) then
          write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-         stop ""
+         stop
       endif
 
@@ -812 +812 @@
 if (ierr.NE.NF_NOERR) then
    WRITE(*,*)'ERROR: Impossible to create the file.'
-   stop ""
+   stop
 endif
 
@@ -963 +963 @@
   hybrid=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading aps"
+    write(*,*) "init2 Error: Failed reading aps"
+    stop
   endif
 
@@ -975 +976 @@
   ierr=NF_INQ_VARID(infid,"bps",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to get bps ID."
+    write(*,*) "init2 Error: Failed to get bps ID."
+    stop
   endif
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading bps"
+    write(*,*) "init2 Error: Failed reading bps"
+    stop
   endif
 endif
@@ -997 +1000 @@
   apbp=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading ap"
+    write(*,*) "Error: Failed reading ap"
+    stop
   endif
 endif
@@ -1015 +1019 @@
   apbp=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading bp"
+    write(*,*) "Error: Failed reading bp"
+    stop
   endif
 endif
@@ -1034 +1039 @@
     sig=.true.
     if (ierr.ne.NF_NOERR) then
-      stop "init2 Error: Failed reading sigma"
+      write(*,*) "init2 Error: Failed reading sigma"
+      stop
     endif
   endif
@@ -1052 +1058 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,aire)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading aire"
+    write(*,*) "init2 Error: Failed reading aire"
+    stop
   endif
   area = .true.
@@ -1070 +1077 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -1087 +1095 @@
              (/layerdimout/),apsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to def_var aps"
+    write(*,*) "init2 Error: Failed to def_var aps"
+    stop
   endif
 
@@ -1093 +1102 @@
   ierr=NF_PUT_VAR_REAL(outfid,apsid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to write aps"
+    write(*,*) "init2 Error: Failed to write aps"
+    stop
   endif
 
@@ -1100 +1110 @@
              (/layerdimout/),bpsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to def_var bps"
+    write(*,*) "init2 Error: Failed to def_var bps"
+    stop
   endif
 
@@ -1106 +1117 @@
   ierr=NF_PUT_VAR_REAL(outfid,bpsid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to write bps"
+    write(*,*) "init2 Error: Failed to write bps"
+    stop
   endif
 
@@ -1114 +1126 @@
              (/interlayerdimout/),apid,ierr)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to def_var ap"
+      write(*,*) "Error: Failed to def_var ap"
+      stop
     endif
 
@@ -1120 +1133 @@
     ierr=NF_PUT_VAR_REAL(outfid,apid,ap)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to write ap"
+      write(*,*) "Error: Failed to write ap"
+      stop
     endif
 
@@ -1127 +1141 @@
              (/interlayerdimout/),bpid,ierr)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to def_var bp"
+      write(*,*) "Error: Failed to def_var bp"
+      stop
     endif
 
@@ -1133 +1148 @@
     ierr=NF_PUT_VAR_REAL(outfid,bpid,bp)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to write bp"
+      write(*,*) "Error: Failed to write bp"
+      stop
     endif
   endif ! of if (apbp)
@@ -1142 +1158 @@
              (/layerdimout/),sigmaid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to def_var sigma"
+    write(*,*) "init2 Error: Failed to def_var sigma"
+    stop
   endif
 ! write sigma
   ierr=NF_PUT_VAR_REAL(outfid,sigmaid,sigma)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to write sigma"
+    write(*,*) "init2 Error: Failed to write sigma"
+    stop
   endif
 endif ! of if (hybrid)
@@ -1160 +1178 @@
            (/londimout,latdimout/),tmpvarid,ierr)
   if (ierr.ne.NF_NOERR) then
-     stop "init2 Error: Failed to def_var aire"
+     write(*,*) "init2 Error: Failed to def_var aire"
+     stop
   endif
   
@@ -1166 +1185 @@
   ierr=NF_PUT_VAR_REAL(outfid,tmpvarid,aire)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to write aire"
+    write(*,*) "init2 Error: Failed to write aire"
+    stop
   endif
 endif ! of if (area)
@@ -1176 +1196 @@
             (/londimout,latdimout/),phisinitid,ierr)
    if (ierr.ne.NF_NOERR) then
-     stop "init2 Error: Failed to def_var phisinit"
+     write(*,*) "init2 Error: Failed to def_var phisinit"
+     stop
    endif
 
@@ -1182 +1203 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinitid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed to write phisinit"
+    write(*,*) "init2 Error: Failed to write phisinit"
+    stop
   endif
 
@@ -1296 +1318 @@
    write(*,*) 'ERROR in change_time_axis: Field <Time> not found'
    print*, NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -1494 +1516 @@
    !write(*,*) 'NF_NOERR', NF_NOERR
    !CALL abort
-   stop ""
+   stop
 endif
 
@@ -1505 +1527 @@
 !    WRITE(*,*) 'NF_NOERR', NF_NOERR
 !    CALL abort
-   stop ""
+   stop
 endif
 

trunk/LMDZ.MARS/util/hrecast.F90

@@ -116 +116 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'ERROR: Pb opening file ',trim(infile)
-   stop ""
+   stop
 endif
 
@@ -189 +189 @@
   write(*,*) 'No variables to process !?'
   write(*,*) 'Might as well stop here'
-  stop ""
+  stop
 else
   write(*,*) ""
@@ -209 +209 @@
 ierr=NF_INQ_DIMID(infid,"latitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get latitude dimension ID"
+  write(*,*) "Error: Failed to get latitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"latitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get latitude ID"
+    write(*,*) "Error: Failed to get latitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,inlatlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get latitude length"
+      write(*,*) "Error: Failed to get latitude length"
+      stop
     else
       allocate(inlat(inlatlength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,inlat)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading latitude"
+        write(*,*) "Error: Failed reading latitude"
+        stop
       endif
     endif
@@ -244 +248 @@
 ierr=NF_INQ_DIMID(infid,"longitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get longitude dimension ID"
+  write(*,*) "Error: Failed to get longitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"longitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get longitude ID"
+    write(*,*) "Error: Failed to get longitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,inlonlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get longitude length"
+      write(*,*) "Error: Failed to get longitude length"
+      stop
     else
       allocate(inlon(inlonlength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,inlon)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading longitude"
+        write(*,*) "Error: Failed reading longitude"
+        stop
       endif
     endif
@@ -279 +287 @@
 ierr=NF_INQ_DIMID(infid,"altitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get altitude dimension ID"
+  write(*,*) "Error: Failed to get altitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"altitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get altitude ID"
+    write(*,*) "Error: Failed to get altitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,altlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get altitude length"
+      write(*,*) "Error: Failed to get altitude length"
+      stop
     else
       allocate(alt(altlength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,alt)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading altitude"
+        write(*,*) "Error: Failed reading altitude"
+        stop
       endif
     endif
@@ -301 +313 @@
 ierr=NF_INQ_DIMID(infid,"Time",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get Time dimension ID"
+  write(*,*) "Error: Failed to get Time dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"Time",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get Time ID"
+    write(*,*) "Error: Failed to get Time ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get Time length"
+      write(*,*) "Error: Failed to get Time length"
+      stop
     else
       allocate(time(timelength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading Time"
+        write(*,*) "Error: Failed reading Time"
+        stop
       endif
     endif
@@ -332 +348 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,sigma)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading sigma"
+    write(*,*) "Error: Failed reading sigma"
+    stop
   endif
 endif
@@ -341 +358 @@
   ierr=NF_INQ_VARID(infid,"aps",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get aps ID"
+    write(*,*) "Error: Failed to get aps ID"
+    stop
   else
     allocate(aps(altlength))
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,aps)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading aps"
+      write(*,*) "Error: Failed reading aps"
+      stop
     endif
   endif
@@ -353 +372 @@
   ierr=NF_INQ_VARID(infid,"bps",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get bps ID"
+    write(*,*) "Error: Failed to get bps ID"
+    stop
   else
     allocate(bps(altlength))
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading bps"
+      write(*,*) "Error: Failed reading bps"
+      stop
     endif
   endif
@@ -377 +398 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,inphisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading phisinit"
+    write(*,*) "Error: Failed reading phisinit"
+    stop
   endif
 endif
@@ -525 +547 @@
 ierr=NF_DEF_DIM(outfid,"longitude",lonlength,lon_dimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define longitude dimension"
+  write(*,*) "Error: Could not define longitude dimension"
+  stop
 endif
 
@@ -531 +554 @@
 ierr=NF_DEF_DIM(outfid,"latitude",latlength,lat_dimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define latitude dimension"
+  write(*,*) "Error: Could not define latitude dimension"
+  stop
 endif
 
@@ -537 +561 @@
 ierr=NF_DEF_DIM(outfid,"altitude",altlength,alt_dimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define altitude dimension"
+  write(*,*) "Error: Could not define altitude dimension"
+  stop
 endif
 
@@ -543 +568 @@
 ierr=NF_DEF_DIM(outfid,"Time",NF_UNLIMITED,time_dimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define latitude dimension"
+  write(*,*) "Error: Could not define latitude dimension"
+  stop
 endif
 
@@ -556 +582 @@
 ierr=NF_DEF_VAR(outfid,"longitude",NF_REAL,1,datashape(1),lon_varid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define longitude variable"
+  write(*,*) "Error: Could not define longitude variable"
+  stop
 endif
 
@@ -563 +590 @@
 ierr=NF_PUT_ATT_TEXT(outfid,lon_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for longitude"
+  write(*,*) "Error: Problem writing long_name for longitude"
+  stop
 endif
 text='degrees_east'
 ierr=NF_PUT_ATT_TEXT(outfid,lon_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for longitude"
+  write(*,*) "Error: Problem writing units for longitude"
+  stop
 endif
 
@@ -575 +604 @@
 ierr=NF_DEF_VAR(outfid,"latitude",NF_REAL,1,datashape(2),lat_varid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define latitude variable"
+  write(*,*) "Error: Could not define latitude variable"
+  stop
 endif
 
@@ -582 +612 @@
 ierr=NF_PUT_ATT_TEXT(outfid,lat_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for latitude"
+  write(*,*) "Error: Problem writing long_name for latitude"
+  stop
 endif
 text='degrees_north'
 ierr=NF_PUT_ATT_TEXT(outfid,lat_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for latitude"
+  write(*,*) "Error: Problem writing units for latitude"
+  stop
 endif
 
@@ -594 +626 @@
 ierr=NF_DEF_VAR(outfid,"altitude",NF_REAL,1,datashape(3),alt_varid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define altitude variable"
+  write(*,*) "Error: Could not define altitude variable"
+  stop
 endif
 
@@ -629 +662 @@
   ierr=NF_DEF_VAR(outfid,"sigma",NF_REAL,1,alt_dimid,sigma_varid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not define sigma variable"
+    write(*,*) "Error: Could not define sigma variable"
+    stop
   endif
 else ! hybrid coordinates
   ierr=NF_DEF_VAR(outfid,"aps",NF_REAL,1,alt_dimid,aps_varid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not define aps variable"
+    write(*,*) "Error: Could not define aps variable"
+    stop
   endif
   ierr=NF_DEF_VAR(outfid,"bps",NF_REAL,1,alt_dimid,bps_varid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not define bps variable"
+    write(*,*) "Error: Could not define bps variable"
+    stop
   endif
 endif
@@ -647 +683 @@
   ierr=NF_PUT_ATT_TEXT(outfid,sigma_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for sigma"
+    write(*,*) "Error: Problem writing long_name for sigma"
+    stop
   endif
 else ! hybrid coordinates
@@ -653 +690 @@
   ierr=NF_PUT_ATT_TEXT(outfid,aps_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for aps"
+    write(*,*) "Error: Problem writing long_name for aps"
+    stop
   endif
   text="hybrid sigma at midlayers"
   ierr=NF_PUT_ATT_TEXT(outfid,bps_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for bps"
+    write(*,*) "Error: Problem writing long_name for bps"
+    stop
   endif
 endif ! of if (have_sigma)
@@ -666 +705 @@
 ierr=NF_DEF_VAR(outfid,"Time",NF_REAL,1,datashape(4),time_varid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not define Time variable"
+  write(*,*) "Error: Could not define Time variable"
+  stop
 endif
 
@@ -683 +723 @@
   ierr=NF_PUT_ATT_TEXT(outfid,time_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for Time"
+    write(*,*) "Error: Problem writing units for Time"
+    stop
   endif
 endif
@@ -694 +735 @@
     ierr=NF_DEF_VAR(outfid,"phisinit",NF_REAL,2,datashape,phisinit_varid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Could not define phisinit variable"
+      write(*,*) "Error: Could not define phisinit variable"
+      stop
     endif
   endif
@@ -701 +743 @@
   ierr=NF_PUT_ATT_TEXT(outfid,phisinit_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for phisinit"
+    write(*,*) "Error: Problem writing long_name for phisinit"
+    stop
   endif
 endif ! of if (have_geopot)
@@ -714 +757 @@
   if (ierr.ne.NF_NOERR) then
     write(*,*) 'Error, failed to get ID for input variable ',trim(var(i))
-    stop ""
+    stop
   endif
   
@@ -738 +781 @@
   if (ierr.ne.NF_NOERR) then
     write(*,*) 'Error, could not define variable ',trim(var(i))
-    stop ""
+    stop
   endif
   
@@ -747 +790 @@
     write(*,*) 'Error, could not get number of attributes for variable ',&
                trim(var(i))
-    stop ""
+    stop
   endif
   ! inititialize j == number of attributes written to output
@@ -766 +809 @@
       if (ierr.ne.NF_NOERR) then
         write(*,*) "Error failed to copy title attribute:",trim(text)
-      stop ""
+        stop
       endif
     endif
@@ -775 +818 @@
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error failed to copy long_name attribute:",trim(text)
-      stop ""
+      stop
     endif
   endif
@@ -789 +832 @@
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error failed to copy units attribute:",trim(text)
-      stop ""
+      stop
     endif
   endif
@@ -805 +848 @@
   ierr=NF_PUT_ATT_REAL(outfid,var_id(i),'missing_value',NF_REAL,1,miss_val)
   if (ierr.ne.NF_NOERR) then
-    stop "Error, failed to write missing_value attribute"
+    write(*,*) "Error, failed to write missing_value attribute"
+    stop
   endif
   
@@ -823 +867 @@
 ierr=NF_ENDDEF(outfid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not switch out of define mode"
+  write(*,*) "Error: Could not switch out of define mode"
+  stop
 endif
 
@@ -829 +874 @@
 ierr=NF_PUT_VAR_REAL(outfid,lon_varid,lon)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not write longitude data to output file"
+  write(*,*) "Error: Could not write longitude data to output file"
+  stop
 endif
 
@@ -835 +881 @@
 ierr=NF_PUT_VAR_REAL(outfid,lat_varid,lat)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not write latitude data to output file"
+  write(*,*) "Error: Could not write latitude data to output file"
+  stop
 endif
 
@@ -841 +888 @@
 ierr=NF_PUT_VAR_REAL(outfid,alt_varid,alt)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not write altitude data to output file"
+  write(*,*) "Error: Could not write altitude data to output file"
+  stop
 endif
 
@@ -848 +896 @@
   ierr=NF_PUT_VAR_REAL(outfid,sigma_varid,sigma)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not write sigma data to output file"
+    write(*,*) "Error: Could not write sigma data to output file"
+    stop
   endif
 else ! hybrid coordinates
   ierr=NF_PUT_VAR_REAL(outfid,aps_varid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not write aps data to output file"
+    write(*,*) "Error: Could not write aps data to output file"
+    stop
   endif
   ierr=NF_PUT_VAR_REAL(outfid,bps_varid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not write bps data to output file"
+    write(*,*) "Error: Could not write bps data to output file"
+    stop
   endif
 endif
 
 ! write time
-ierr=NF_PUT_VARA_REAL(outfid,time_varid,1,timelength,time)
-if (ierr.ne.NF_NOERR) then
-  stop "Error: Could not write Time data to output file"
+ierr=NF_PUT_VARA_REAL(outfid,time_varid,(/1/),(/timelength/),time)
+if (ierr.ne.NF_NOERR) then
+  write(*,*) "Error: Could not write Time data to output file"
+  stop
 endif
 
@@ -880 +932 @@
   ierr=NF_INQ_VARID(infid,"phisinit",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get phisinit ID"
+    write(*,*) "Error: Failed to get phisinit ID"
+    stop
   endif
   ! Get physinit
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,in_2d_data)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading input phisinit"
+    write(*,*) "Error: Failed reading input phisinit"
+    stop
   endif
   
@@ -903 +957 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinit_varid,out_2d_data)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Could not write phisinit data to output file"
+    write(*,*) "Error: Could not write phisinit data to output file"
+    stop
   endif
 endif ! of if (have_geopot)
@@ -1073 +1128 @@
 
        REAL airen ((imnmx2+1)*(jmnmx2+1)) ! aire dans la nouvelle grille
-       REAL airentotn   ! aire totale pole nord dans la nouvelle grille
-       REAL airentots   ! aire totale pole sud dans la nouvelle grille
+       REAL airentotn   ! aire totale pole nord dans la nouvelle grille
+       REAL airentots   ! aire totale pole sud dans la nouvelle grille
 !    Info sur les ktotal intersection entre les cases new/old grille
 

trunk/LMDZ.MARS/util/localtime.F90

@@ -12 +12 @@
 include "netcdf.inc" ! NetCDF definitions
 
-character (len=50)  file
+character (len=100)  file
 ! file(): input file(s) names(s)
 character (len=30), dimension(16) :: notprocessed
 ! notprocessed(): names of the (16) variables that won't be processed
-character (len=50), dimension(:), allocatable :: var
+character (len=100), dimension(:), allocatable :: var
 ! var(): name(s) of variable(s) that will be processed
 character (len=100) :: tmpvar,long_name,units
@@ -121 +121 @@
 if (len_trim(file).eq.0) then
    write(*,*) "no file... game over"
-   stop ""
+   stop
 endif
 
@@ -132 +132 @@
    write(*,*) 'ERROR: Pb opening file '//trim(file)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -140 +140 @@
    write(*,*) 'ERROR: Pb with NF_INQ_NVARS'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -199 +199 @@
 write(*,*) "(an empty line , i.e: just <Enter>, implies end of list)"
 nbvar=0
-read(*,'(a50)') tmpvar
+read(*,'(a100)') tmpvar
 do while ((tmpvar/=' ').AND.(trim(tmpvar)/='all'))
    nbvar=nbvar+1
    var(nbvar)=tmpvar
-   read(*,'(a50)') tmpvar
+   read(*,'(a100)') tmpvar
 enddo
 
@@ -219 +219 @@
 else if(nbvar==0) then
    write(*,*) "no variable... game over"
-   stop ""
+   stop
 endif ! of if (tmpvar=="all")
 
@@ -238 +238 @@
       write(*,*) 'ERROR: Pb opening file '//trim(file)
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 
@@ -251 +251 @@
          if (ierr.NE.NF_NOERR) then
             write(*,*) 'ERROR: Dimension <latitude> and <lat> is missing in file '//trim(file)
-            stop ""  
+            stop  
          endif
    endif
@@ -260 +260 @@
       if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <latitude> and <lat> is missing in file '//trim(file)
-        stop ""  
+        stop  
       endif
    endif
@@ -272 +272 @@
          if (ierr.NE.NF_NOERR) then
             write(*,*) 'ERROR: Dimension <longitude> and <lon> is missing in file '//trim(file)
-            stop ""  
+            stop  
          endif 
    endif
@@ -281 +281 @@
       if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <longitude> and <lon> is missing in file '//trim(file)
-        stop ""  
+        stop  
       endif
    endif
@@ -290 +290 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <altitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_VARID(nid,"altitude",altvar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <altitude> is missing in file '//trim(file)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,altdim,altlen)
@@ -321 +321 @@
          ctllen=0
       endif
-      !stop ""  
+      !stop  
    endif
 !   ierr=NF_INQ_DIMID(nid,"controle_axe",ctldim)
 !   if (ierr.NE.NF_NOERR) then
-!      !stop ""  
+!      !stop  
 !   endif
    ierr=NF_INQ_VARID(nid,"controle",ctlvar)
@@ -331 +331 @@
       write(*,*) 'Field <controle> is missing in file '//trim(file)
       ctllen=0
-      !stop ""
+      !stop
    else
       ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
@@ -414 +414 @@
      if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Dimension <Time> or is missing in file'//trim(file)
-        stop ""
+        stop
      endif
    endif
@@ -423 +423 @@
      if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <Time> is missing in file '//trim(file)
-        stop ""
+        stop
      endif
    endif
@@ -556 +556 @@
 
      do it=1,timelen_lt
-        ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,lt_hour(it))
+        ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),lt_hour(it))
      enddo
    else
      do it=1,timelen_lt
-        ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,lt_out(it))
+        ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),lt_out(it))
         if (ierr.NE.NF_NOERR) then
           write(*,*) "Error , failed to write Time"
@@ -582 +582 @@
       if (ierr.NE.NF_NOERR) then
          write(*,*) 'ERROR: Field <',var(j),'> not found in file'//file
-         stop ""
+         stop
       endif
       ierr=nf_inq_varndims(nid,varid,ndim)
@@ -718 +718 @@
       if (ierr.ne.NF_NOERR) then
          write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-         stop ""
+         stop
       endif
 
@@ -820 +820 @@
    WRITE(*,*)'ERROR: Impossible to create the file ',trim(filename)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -831 +831 @@
    WRITE(*,*)'initiate: error failed to define dimension <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "longitude", lonlen, londimout)
@@ -837 +837 @@
    WRITE(*,*)'initiate: error failed to define dimension <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "altitude", altlen, altdimout)
@@ -843 +843 @@
    WRITE(*,*)'initiate: error failed to define dimension <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 if (size(ctl).ne.0) then
@@ -850 +850 @@
     WRITE(*,*)'initiate: error failed to define dimension <index>'
     write(*,*) NF_STRERROR(ierr)
-    stop ""
+    stop
   endif
 endif
@@ -857 +857 @@
    WRITE(*,*)'initiate: error failed to define dimension <Time>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -867 +867 @@
    WRITE(*,*)'initiate: error failed to switch out of define mode'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -888 +888 @@
    WRITE(*,*)'initiate: error failed writing variable <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -903 +903 @@
    WRITE(*,*)'initiate: error failed writing variable <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -926 +926 @@
    WRITE(*,*)'initiate: error failed writing variable <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -948 +948 @@
       WRITE(*,*)'initiate: error failed writing variable <controle>'
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 endif
@@ -1014 +1014 @@
   hybrid=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading aps"
+    write(*,*) "init2 Error: Failed reading aps"
+    stop
   endif
 
@@ -1031 +1032 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "init2 Error: Failed reading bps"
+      write(*,*) "init2 Error: Failed reading bps"
+      stop
     endif
   endif
@@ -1050 +1052 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -1069 +1072 @@
                (/layerdimout/),apsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var aps"
+    write(*,*) "Error: Failed to def_var aps"
+    stop
   endif
 
@@ -1075 +1079 @@
   ierr=NF_PUT_VAR_REAL(outfid,apsid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write aps"
+    write(*,*) "Error: Failed to write aps"
+    stop
   endif
   
@@ -1082 +1087 @@
                (/layerdimout/),bpsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var bps"
+    write(*,*) "Error: Failed to def_var bps"
+    stop
   endif
 
@@ -1088 +1094 @@
   ierr=NF_PUT_VAR_REAL(outfid,bpsid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write bps"
+    write(*,*) "Error: Failed to write bps"
+    stop
   endif
   
@@ -1105 +1112 @@
               (/londimout,latdimout/),phisinitid,ierr)
   if (ierr.ne.NF_NOERR) then
-     stop "Error: Failed to def_var phisinit"
+     write(*,*) "Error: Failed to def_var phisinit"
+     stop
   endif
 
@@ -1111 +1119 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinitid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write phisinit"
+    write(*,*) "Error: Failed to write phisinit"
+    stop
   endif
 
@@ -1226 +1235 @@
      !write(*,*) 'NF_NOERR', NF_NOERR
      !CALL abort
-     stop ""
+     stop
   endif
 endif
@@ -1238 +1247 @@
   !    WRITE(*,*) 'NF_NOERR', NF_NOERR
   !    CALL abort
-     stop ""
+     stop
   endif
 endif

trunk/LMDZ.MARS/util/lslin.F90

@@ -131 +131 @@
    write(*,*) 'Failed to open file '//infile
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -149 +149 @@
     write(*,*) "Failed to get number of dims for variable number:",start_var
     write(*,*) NF_STRERROR(ierr)
-    stop ""
+    stop
   endif
 enddo
@@ -188 +188 @@
       write(*,*) 'ERROR: Field <latitude> is missing'
       write(*,*) NF_STRERROR(ierr)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_DIMLEN(nid,latdim,latlen)
@@ -198 +198 @@
       write(*,*) 'ERROR: Field <longitude> is missing'
       write(*,*) NF_STRERROR(ierr)
-      stop "" 
+      stop 
    endif
    ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
@@ -208 +208 @@
       write(*,*) 'ERROR: Field <altitude> is missing'
       write(*,*) NF_STRERROR(ierr)
-!     stop ""
+!     stop
        altlen = 1 
    else
@@ -239 +239 @@
       write(*,*) ' Dimension <index> is missing in file '//trim(infile)
       ctllen=0
-      !stop ""  
+      !stop  
    endif
    ierr=NF_INQ_VARID(nid,"controle",ctlvar)
@@ -245 +245 @@
       write(*,*) 'Field <controle> is missing in file '//trim(infile)
       ctllen=0
-      !stop ""
+      !stop
    else
       ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
@@ -289 +289 @@
    write(*,*) 'ERROR: Dimension <Time> is missing'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr=NF_INQ_VARID(nid,"Time",timevar)
@@ -295 +295 @@
    write(*,*) 'ERROR: Field <Time> is missing'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -407 +407 @@
 
 do k=1,Nls
-  ierr= NF_PUT_VARA_REAL(nout,timevarout,k,1,timels(k))
+  ierr= NF_PUT_VARA_REAL(nout,timevarout,(/k/),(/1/),timels(k))
 enddo
 
@@ -426 +426 @@
       write(*,*) 'ERROR: Field <',var(j),'> not found'
       write(*,*) NF_STRERROR(ierr)
-      stop "Better stop now..."
+      write(*,*) "Better stop now..."
+      stop
    endif
    
@@ -589 +590 @@
    if (ierr.ne.NF_NOERR) then
      write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-     stop ""
+     stop
    endif
 
@@ -716 +717 @@
    WRITE(*,*)'ERROR: Impossible to create the file ',trim(filename)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -727 +728 @@
    WRITE(*,*)'initiate: error failed to define dimension <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "longitude", lonlen, londimout)
@@ -733 +734 @@
    WRITE(*,*)'initiate: error failed to define dimension <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "altitude", altlen, altdimout)
@@ -739 +740 @@
    WRITE(*,*)'initiate: error failed to define dimension <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 if (size(ctl).ne.0) then
@@ -746 +747 @@
     WRITE(*,*)'initiate: error failed to define dimension <index>'
     write(*,*) NF_STRERROR(ierr)
-    stop ""
+    stop
   endif
 endif
@@ -753 +754 @@
    WRITE(*,*)'initiate: error failed to define dimension <Time>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -761 +762 @@
    WRITE(*,*)'initiate: error failed to switch out of define mode'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -785 +786 @@
    WRITE(*,*)'initiate: error failed writing variable <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -799 +800 @@
    WRITE(*,*)'initiate: error failed writing variable <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -822 +823 @@
    WRITE(*,*)'initiate: error failed writing variable <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -844 +845 @@
       WRITE(*,*)'initiate: error failed writing variable <controle>'
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 endif
@@ -917 +918 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,aps)
     if (ierr.ne.NF_NOERR) then
-     stop "init2 error: Failed reading aps"
+     write(*,*) "init2 error: Failed reading aps"
+     stop
      aps_ok=.false.
     endif
@@ -945 +947 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "init2 Error: Failed reading bps"
+      write(*,*) "init2 Error: Failed reading bps"
+      stop
       bps_ok=.false.
     endif
@@ -966 +969 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -985 +989 @@
                (/altdimout/),apsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var aps"
+    write(*,*) "Error: Failed to def_var aps"
+    stop
   endif
 
@@ -991 +996 @@
   ierr=NF_PUT_VAR_REAL(outfid,apsid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write aps"
+    write(*,*) "Error: Failed to write aps"
+    stop
   endif
 ENDIF ! of IF (aps_ok)
@@ -1000 +1006 @@
                (/altdimout/),bpsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var bps"
+    write(*,*) "Error: Failed to def_var bps"
+    stop
   endif
 
@@ -1006 +1013 @@
   ierr=NF_PUT_VAR_REAL(outfid,bpsid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write bps"
+    write(*,*) "Error: Failed to write bps"
+    stop
   endif
 ENDIF ! of IF (bps_ok)
@@ -1020 +1028 @@
               (/londimout,latdimout/),phisinitid,ierr)
   if (ierr.ne.NF_NOERR) then
-     stop "Error: Failed to def_var phisinit"
+     write(*,*) "Error: Failed to def_var phisinit"
+     stop
   endif
 
@@ -1026 +1035 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinitid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write phisinit"
+    write(*,*) "Error: Failed to write phisinit"
+    stop
   endif
 
@@ -1141 +1151 @@
 !   print*,'missing_value: valid_range attribution failed'
 !   print*, NF_STRERROR(ierr)
-!   stop ""
+!   stop
 !endif
 !*********************************************************
@@ -1150 +1160 @@
    print*, 'missing_value: missing value attribution failed'
    print*, NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 

trunk/LMDZ.MARS/util/solzenangle.F90

@@ -17 +17 @@
 character (len=30), dimension(16) :: notprocessed
 ! notprocessed(): names of the (16) variables that won't be processed
-character (len=50), dimension(:), allocatable :: var
+character (len=100), dimension(:), allocatable :: var
 ! var(): name(s) of variable(s) that will be processed
 character (len=100) :: tmpvar,long_name,units
@@ -141 +141 @@
 if (len_trim(file).eq.0) then
    write(*,*) "no file... game over"
-   stop ""
+   stop
 endif
 
@@ -152 +152 @@
    write(*,*) 'ERROR: Pb opening file '//trim(file)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -160 +160 @@
    write(*,*) 'ERROR: Pb with NF_INQ_NVARS'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -218 +218 @@
 write(*,*) "(an empty line , i.e: just <Enter>, implies end of list)"
 nbvar=0
-read(*,'(a50)') tmpvar
+read(*,'(a100)') tmpvar
 do while ((tmpvar/=' ').AND.(trim(tmpvar)/='all'))
    nbvar=nbvar+1
    var(nbvar)=tmpvar
-   read(*,'(a50)') tmpvar
+   read(*,'(a100)') tmpvar
 enddo
 
@@ -238 +238 @@
 else if(nbvar==0) then
    write(*,*) "no variable... game over"
-   stop ""
+   stop
 endif ! of if (tmpvar=="all")
 
@@ -252 +252 @@
    write(*,*) 'ERROR: Pb opening file '//trim(file)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -262 +262 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <latitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_VARID(nid,"latitude",latvar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <latitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_DIMLEN(nid,latdim,latlen)
@@ -275 +275 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <longitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_VARID(nid,"longitude",lonvar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <longitude> is missing in file'//trim(file)
-      stop "" 
+      stop 
    endif
    ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
@@ -288 +288 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <altitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_VARID(nid,"altitude",altvar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <altitude> is missing in file'//trim(file)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,altdim,altlen)
@@ -315 +315 @@
       write(*,*) ' Dimension <index> is missing in file '//trim(file)
       ctllen=0
-      !stop ""  
+      !stop  
    endif
    ierr=NF_INQ_VARID(nid,"controle",ctlvar)
@@ -321 +321 @@
       write(*,*) 'Field <controle> is missing in file '//trim(file)
       ctllen=0
-      !stop ""
+      !stop
    else
       ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
@@ -401 +401 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <Time> is missing in file'//trim(file)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_VARID(nid,"Time",timevar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <Time> is missing in file'//trim(file)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,timedim,timelen)
@@ -583 +583 @@
 
      do it=1,nsol
-       ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,intsol(it)*24.)
+       ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),intsol(it)*24.)
      enddo
    else
      do it=1,nsol
-       ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,intsol(it))
+       ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),intsol(it))
        if (ierr.NE.NF_NOERR) then
          write(*,*) "Error , failed to write Time"
@@ -609 +609 @@
       if (ierr.NE.NF_NOERR) then
          write(*,*) 'ERROR: Field <',var(j),'> not found in file'//file
-         stop ""
+         stop
       endif
       ierr=nf_inq_varndims(nid,varid,ndim)
@@ -768 +768 @@
        if (ierr.ne.NF_NOERR) then
          write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-         stop ""
+         stop
        endif
 
@@ -814 +814 @@
    if (ierr.ne.NF_NOERR) then
       write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 
@@ -929 +929 @@
    WRITE(*,*)'ERROR: Impossible to create the file ',trim(filename)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -940 +940 @@
    WRITE(*,*)'initiate: error failed to define dimension <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "longitude", lonlen, londimout)
@@ -946 +946 @@
    WRITE(*,*)'initiate: error failed to define dimension <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "altitude", altlen, altdimout)
@@ -952 +952 @@
    WRITE(*,*)'initiate: error failed to define dimension <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 if (size(ctl).ne.0) then
@@ -959 +959 @@
     WRITE(*,*)'initiate: error failed to define dimension <index>'
     write(*,*) NF_STRERROR(ierr)
-    stop ""
+    stop
   endif
 endif
@@ -966 +966 @@
    WRITE(*,*)'initiate: error failed to define dimension <Time>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -973 +973 @@
    WRITE(*,*)'initiate: error failed to define dimension <GCM_layers>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -981 +981 @@
    WRITE(*,*)'initiate: error failed to switch out of define mode'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 !==============================================================================
@@ -992 +992 @@
 ierr = NF_REDEF (nout)
 ierr = NF_PUT_ATT_TEXT(nout,timevarout,'comment',135,&
-"The true unit of the variable Time is the integer value &
- of sol at lon=0deg. A false unit is put to enable reading from Grads or Ferret")
+"The true unit of the variable Time is the integer value "//&
+"of sol at lon=0deg. A false unit is put to enable reading from Grads or Ferret")
 ! End netcdf define mode
 ierr = NF_ENDDEF(nout)
@@ -1008 +1008 @@
    WRITE(*,*)'initiate: error failed writing variable <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -1022 +1022 @@
    WRITE(*,*)'initiate: error failed writing variable <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -1045 +1045 @@
    WRITE(*,*)'initiate: error failed writing variable <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -1067 +1067 @@
       WRITE(*,*)'initiate: error failed writing variable <controle>'
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 endif
@@ -1134 +1134 @@
   hybrid=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading aps"
+    write(*,*) "init2 Error: Failed reading aps"
+    stop
   endif
 
@@ -1151 +1152 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "init2 Error: Failed reading bps"
+      write(*,*) "init2 Error: Failed reading bps"
+      stop
     endif
   endif
@@ -1170 +1172 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -1189 +1192 @@
                (/layerdimout/),apsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var aps"
+    write(*,*) "Error: Failed to def_var aps"
+    stop
   endif
 
@@ -1195 +1199 @@
   ierr=NF_PUT_VAR_REAL(outfid,apsid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write aps"
+    write(*,*) "Error: Failed to write aps"
+    stop
   endif
   
@@ -1202 +1207 @@
                (/layerdimout/),bpsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var bps"
+    write(*,*) "Error: Failed to def_var bps"
+    stop
   endif
 
@@ -1208 +1214 @@
   ierr=NF_PUT_VAR_REAL(outfid,bpsid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write bps"
+    write(*,*) "Error: Failed to write bps"
+    stop
   endif
   
@@ -1225 +1232 @@
               (/londimout,latdimout/),phisinitid,ierr)
   if (ierr.ne.NF_NOERR) then
-     stop "Error: Failed to def_var phisinit"
+     write(*,*) "Error: Failed to def_var phisinit"
+     stop
   endif
 
@@ -1231 +1239 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinitid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write phisinit"
+    write(*,*) "Error: Failed to write phisinit"
+    stop
   endif
 
@@ -1347 +1356 @@
      !write(*,*) 'NF_NOERR', NF_NOERR
      !CALL abort
-     stop ""
+     stop
   endif
 endif
@@ -1360 +1369 @@
   !    WRITE(*,*) 'NF_NOERR', NF_NOERR
   !    CALL abort
-     stop ""
+     stop
   endif
 endif

trunk/LMDZ.MARS/util/streamfunction.F90

@@ -129 +129 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'ERROR: Pb opening input file'
-   stop ""
+   stop
 endif
 
@@ -143 +143 @@
 ierr=NF_INQ_DIMID(infid,"latitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get latitude dimension ID"
+  write(*,*) "Error: Failed to get latitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"latitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get latitude ID"
+    write(*,*) "Error: Failed to get latitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,latlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get latitude length"
+      write(*,*) "Error: Failed to get latitude length"
+      stop
     else
       allocate(lat(latlength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,lat)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading latitude"
+        write(*,*) "Error: Failed reading latitude"
+        stop
       endif
     endif
@@ -165 +169 @@
 ierr=NF_INQ_DIMID(infid,"longitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get longitude dimension ID"
+  write(*,*) "Error: Failed to get longitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"longitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get longitude ID"
+    write(*,*) "Error: Failed to get longitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,lonlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get longitude length"
+      write(*,*) "Error: Failed to get longitude length"
+      stop
     else
       allocate(lon(lonlength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,lon)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading longitude"
+        write(*,*) "Error: Failed reading longitude"
+        stop
       endif
     endif
@@ -187 +195 @@
 ierr=NF_INQ_DIMID(infid,"Time",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get Time dimension ID"
+  write(*,*) "Error: Failed to get Time dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"Time",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get Time ID"
+    write(*,*) "Error: Failed to get Time ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get Time length"
+      write(*,*) "Error: Failed to get Time length"
+      stop
     else
       allocate(time(timelength))
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading Time"
+        write(*,*) "Error: Failed reading Time"
+        stop
       endif
     endif
@@ -209 +221 @@
 ierr=NF_INQ_DIMID(infid,"altitude",tmpdimid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get altitude dimension ID"
+  write(*,*) "Error: Failed to get altitude dimension ID"
+  stop
 else
   ierr=NF_INQ_VARID(infid,"altitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-     stop "Error: Failed to get altitude length"
+     write(*,*) "Error: Failed to get altitude length"
+     stop
   else
      ierr=NF_INQ_DIMLEN(infid,tmpdimid,altlength)
      if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get altitude length"
+        write(*,*) "Error: Failed to get altitude length"
+        stop
      else   
         allocate(alt(altlength))
         ierr=NF_GET_VAR_REAL(infid,tmpvarid,alt)
         if (ierr.ne.NF_NOERR) then
-           stop "Error: Failed reading Altitude"
+           write(*,*) "Error: Failed reading Altitude"
+           stop
         endif
      endif
@@ -235 +251 @@
 ierr=NF_INQ_VARID(infid,"aps",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-   stop "Error: Failed to get aps ID"
+   write(*,*) "Error: Failed to get aps ID"
+   stop
 else
    allocate(aps(altlength))
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,aps)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading aps"
+      write(*,*) "Error: Failed reading aps"
+      stop
    endif
 endif
@@ -247 +265 @@
 ierr=NF_INQ_VARID(infid,"bps",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-   stop "Error: Failed to get bps ID"
+   write(*,*) "Error: Failed to get bps ID"
+   stop
 else
    allocate(bps(altlength))
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading bps"
+      write(*,*) "Error: Failed reading bps"
+      stop
    endif
 endif
@@ -259 +279 @@
 ierr=NF_INQ_VARID(infid,"ps",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-   stop "Error: Failed to get ps ID"
+   write(*,*) "Error: Failed to get ps ID"
+   stop
 else
    allocate(ps(lonlength,latlength,timelength))
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,ps)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading ps"
+      write(*,*) "Error: Failed reading ps"
+      stop
    endif
 endif
@@ -284 +306 @@
 ierr=NF_INQ_VARID(infid,"u",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get u ID"
+  write(*,*) "Error: Failed to get u ID"
+  stop
 else
   allocate(u(lonlength,latlength,altlength,timelength))
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,u)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading zonal wind"
+    write(*,*) "Error: Failed reading zonal wind"
+    stop
   endif
 endif
@@ -296 +320 @@
 ierr=NF_INQ_VARID(infid,"v",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get v ID"
+  write(*,*) "Error: Failed to get v ID"
+  stop
 else
   allocate(v(lonlength,latlength,altlength,timelength))
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,v)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading zonal wind"
+    write(*,*) "Error: Failed reading zonal wind"
+    stop
   endif
 endif
@@ -322 +348 @@
      if (ierr.ne.NF_NOERR) then
         write(*,*) "Problem: Could not find/open these files"
-        stop "Might as well stop here"
+        write(*,*) "Might as well stop here"
+        stop
      endif
   endif
@@ -329 +356 @@
    ierr=NF_INQ_VARID(infid2,"aire",tmpvarid)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get aire ID"
+      write(*,*) "Error: Failed to get aire ID"
+      stop
    endif
    ! Get aire
    ierr=NF_GET_VAR_REAL(infid2,tmpvarid,aire)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading aire"
+      write(*,*) "Error: Failed reading aire"
+      stop
    endif
    ! Close file
@@ -342 +371 @@
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,aire)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading aire"
+      write(*,*) "Error: Failed reading aire"
+      stop
    endif
 endif
@@ -363 +393 @@
      if (ierr.ne.NF_NOERR) then
         write(*,*) "Problem: Could not find/open these files"
-        stop "Might as well stop here"
+        write(*,*) "Might as well stop here"
+        stop
      endif
   endif
@@ -371 +402 @@
    ierr=NF_INQ_VARID(infid2,"phisinit",tmpvarid)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get phisinit ID"
+      write(*,*) "Error: Failed to get phisinit ID"
+      stop
    endif
    ! Get phisinit
    ierr=NF_GET_VAR_REAL(infid2,tmpvarid,phisinit)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading phisinit"
+      write(*,*) "Error: Failed reading phisinit"
+      stop
    endif
    ! Close file
@@ -384 +417 @@
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading phisinit"
+      write(*,*) "Error: Failed reading phisinit"
+      stop
    endif
 endif
@@ -407 +441 @@
      if (ierr.ne.NF_NOERR) then
         write(*,*) "Problem: Could not find/open these files"
-        stop "Might as well stop here"
+        write(*,*) "Might as well stop here"
+        stop
      endif
   endif
@@ -415 +450 @@
    ierr=NF_INQ_VARID(infid2,"cv",tmpvarid)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get cv ID"
+      write(*,*) "Error: Failed to get cv ID"
+      stop
    endif
    ! Get cv
    ierr=NF_GET_VAR_REAL(infid2,tmpvarid,cv)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading cv"
+      write(*,*) "Error: Failed reading cv"
+      stop
    endif
    ! Close file
@@ -428 +465 @@
    ierr=NF_GET_VAR_REAL(infid,tmpvarid,cv)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading cv"
+      write(*,*) "Error: Failed reading cv"
+      stop
    endif
 endif
@@ -436 +474 @@
 ierr=NF_INQ_VARID(infid,"temp",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get temp ID"
+  write(*,*) "Error: Failed to get temp ID"
+  stop
 else
   allocate(temp(lonlength,latlength,altlength,timelength))
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,temp)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading temperature"
+    write(*,*) "Error: Failed reading temperature"
+    stop
   endif
 endif
@@ -464 +504 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,rho)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading mass density"
+    write(*,*) "Error: Failed reading mass density"
+    stop
   endif
 endif
@@ -666 +707 @@
 lon_fake(1)=0.
 #ifdef NC_DOUBLE
-        ierr= NF_PUT_VARA_DOUBLE(nout,lonvarout,1,1,lon_fake)
+        ierr= NF_PUT_VARA_DOUBLE(nout,lonvarout,(/1/),(/1/),lon_fake)
 #else
-        ierr= NF_PUT_VARA_REAL(nout,lonvarout,1,1,lon_fake)
+        ierr= NF_PUT_VARA_REAL(nout,lonvarout,(/1/),(/1/),lon_fake)
 #endif
 
 ! 3.1.3 New time dimension/value in output file
 #ifdef NC_DOUBLE
-        ierr= NF_PUT_VARA_DOUBLE(nout,timevarout,1,1,lon_fake)
+        ierr= NF_PUT_VARA_DOUBLE(nout,timevarout,(/1/),(/1/),lon_fake)
 #else
-        ierr= NF_PUT_VARA_REAL(nout,timevarout,1,1,lon_fake)
+        ierr= NF_PUT_VARA_REAL(nout,timevarout,(/1/),(/1/),lon_fake)
 #endif
 
@@ -690 +731 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable psi'
-   stop ""
+   stop
 endif
 
@@ -705 +746 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable momave'
-   stop ""
+   stop
 endif
 
@@ -721 +762 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable u'
-   stop ""
+   stop
 endif
 
@@ -737 +778 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable v'
-   stop ""
+   stop
 endif
 
@@ -752 +793 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable rho'
-   stop ""
+   stop
 endif
 
@@ -767 +808 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable temp'
-   stop ""
+   stop
 endif
 
@@ -782 +823 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable ps'
-   stop ""
+   stop
 endif
 
@@ -798 +839 @@
 if (ierr.ne.NF_NOERR) then
    write(*,*) 'Error, could not define variable phisinit'
-   stop ""
+   stop
 endif
 
@@ -874 +915 @@
 if (ierr.NE.NF_NOERR) then
    WRITE(*,*)'ERROR: Impossible to create the file.'
-   stop ""
+   stop
 endif
 
@@ -983 +1024 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,aps)
   if (ierr.ne.NF_NOERR) then
-   stop "error: Failed reading aps"
+   write(*,*) "error: Failed reading aps"
+   stop
   endif
   aps_ok=.true.
@@ -997 +1039 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading bps"
+    write(*,*) "Error: Failed reading bps"
+    stop
   endif
   bps_ok=.true.
@@ -1015 +1058 @@
         (/altdimout/),apsid,ierr)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to def_var aps"
+      write(*,*) "Error: Failed to def_var aps"
+      stop
    endif
 
@@ -1025 +1069 @@
 #endif
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to write aps"
+      write(*,*) "Error: Failed to write aps"
+      stop
    endif
 ENDIF 
@@ -1034 +1079 @@
         (/altdimout/),bpsid,ierr)
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to def_var bps"
+      write(*,*) "Error: Failed to def_var bps"
+      stop
    endif
 
@@ -1044 +1090 @@
 #endif
    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to write bps"
+      write(*,*) "Error: Failed to write bps"
+      stop
    endif
 ENDIF

trunk/LMDZ.MARS/util/xvik/fit_Iceinertia_MONSicedepth.F

@@ -132 +132 @@
     
        maxiceidiff = 3000  ! maximum allowed difference between N and S best fit ice thermal inertia
-       deltaicei= 20     ! step in ice thermal inertia
+       deltaicei= 20         ! step in ice thermal inertia
        maxiceddiff = 30.e-4  ! maximum allowed difference between N and S best fit ice depth
        deltaiced= 0.5e-4 ! step in ice depth coefficient
@@ -332 +332 @@
         read(13,*) solgcm(n), patm(n,3), pcapn(n,3),pcaps(n,3)
         read(14,*) solgcm(n), patm(n,4), pcapn(n,4),pcaps(n,4)
-        
+        
         elseif (time_unit == 2) then
         read(21,*) lsgcm(n),  pvl_gcm(n,1)
@@ -358 +358 @@
         read(13,*) solgcm(n), lsgcm(n), patm(n,3), pcapn(n,3),pcaps(n,3)
         read(14,*) solgcm(n), lsgcm(n), patm(n,4), pcapn(n,4),pcaps(n,4)
-        
-        else 
-        write(*,*) 'Wrong integer for xvik files format :',
+        
+        else 
+        write(*,*) 'Wrong integer for xvik files format :',
      &' must be 1, 2 or 3'
         stop
 
         endif
-        
+        
 c       Checking total CO2 inventory for all runs :
         do i=1,4
@@ -381 +381 @@
       end do
 
-        close(11)
-        close(12)
-        close(13)
-        close(14)
-        close(21)
-        close(22)
-        close(23)
-        close(24)
-        
-        
+        close(11)
+        close(12)
+        close(13)
+        close(14)
+        close(21)
+        close(22)
+        close(23)
+        close(24)
+        
+        
 c      Smoothing simulated GCM Viking 1 pressure curves
 c      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -493 +493 @@
 c     to plot COST(DN,DS) (plot_test=1), loops must be in this order : DN->DS->IN->IS            
         write(*,*)'   DN        DS      cost       IN       IS       Ps'
-        icedn=icedmin ! initialization          
+        icedn=icedmin ! initialization                
         do while (icedn.le.icedmax) ! loop on northern ice depth coefficient
 !         albn=albmin ! initialization
@@ -510 +510 @@
               cost =0. !initialization
               do n=1,nsol  
-                !write(*,*) n
+                !write(*,*) n
                ! Pressure corresponding to Northern cap
                pcapn_new=pcapn_sm(n,refrun)+
@@ -527 +527 @@
                ! and observed pressures at VL1 site
                cost= cost+ ( pvl_obs(n) - pvl1)**2 
-                
+                
               end do ! of do n=1,nsol
               ! store parameters which lead to minimum cost
@@ -539 +539 @@
                  icedsfit=iceds
               end if
-           
+           
               if(cost.lt.cost4plot) then
 c                RMS, best pressure and best albedos for these icedsivities value 
@@ -575 +575 @@
 c     to plot COST(IN,IS) (plot_test=2), loops must be in this order : IN->IS->DN->DS            
         write(*,*)'   DN        DS      cost       IN       IS       Ps'
-        write(33,*)'   DN        DS      cost       IN       IS     
-     &Ps'       
-        icein=iceimin ! initialization          
+        write(33,*)'   DN        DS      cost       IN       IS     Ps'
+        icein=iceimin ! initialization                
         do while (icein.le.iceimax) ! loop on northern ice depth coefficient
 !         albn=albmin ! initialization
@@ -594 +593 @@
               cost =0. !initialization
               do n=1,nsol  
-                !write(*,*) n
+                !write(*,*) n
                ! Pressure corresponding to Northern cap
                pcapn_new=pcapn_sm(n,refrun)+
@@ -611 +610 @@
                ! and observed pressures at VL1 site
                cost= cost+ ( pvl_obs(n) - pvl1)**2 
-                
+                
               end do ! of do n=1,nsol
               ! store parameters which lead to minimum cost
@@ -623 +622 @@
                  icedsfit=iceds
               end if
-           
+           
               if(cost.lt.cost4plot) then
 c                RMS, best pressure and best albedos for these icedsivities value 
@@ -683 +682 @@
      &      + (fonc2s(icedsfit) -fonc2s(icedgcm_ref)) * dpdes(n)
            pcaps_new= max(pcaps_new,0.)
-        
+        
       call sol2ls(sol(n),solconv)   
       
@@ -722 +721 @@
 
 
-      if (nbig.lt.3*nmax) stop 'Must increase nbig in runave'
+      if (nbig.lt.3*nmax) then
+        write(*,*) 'Must increase nbig in runave'
+        stop
+      endif
 
 c     Reindexation des donnees

trunk/LMDZ.MARS/util/zrecast.F90

@@ -176 +176 @@
    write(*,*) 'ERROR: Pb opening file ',trim(infile)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -270 +270 @@
   write(*,*) 'No variables to process !?'
   write(*,*) 'Might as well stop here'
-  stop ""
+  stop
 else
   write(*,*) ""
@@ -294 +294 @@
   ierr=NF_INQ_DIMID(infid,"lat",tmpdimid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get lat dimension ID"
+    write(*,*) "Error: Failed to get lat dimension ID"
+    stop
   else
     ierr=NF_INQ_VARID(infid,"lat",tmpvarid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get lat ID"
+      write(*,*) "Error: Failed to get lat ID"
+      stop
     else
       ierr=NF_INQ_DIMLEN(infid,tmpdimid,latlength)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get lat length"
+        write(*,*) "Error: Failed to get lat length"
+        stop
       else
         allocate(lat(latlength),stat=ierr)
@@ -312 +315 @@
         ierr=NF_GET_VAR_REAL(infid,tmpvarid,lat)
         if (ierr.ne.NF_NOERR) then
-          stop "Error: Failed reading lat"
+          write(*,*) "Error: Failed reading lat"
+          stop
         endif
       endif
@@ -320 +324 @@
   ierr=NF_INQ_VARID(infid,"latitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get latitude ID"
+    write(*,*) "Error: Failed to get latitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,latlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get latitude length"
+      write(*,*) "Error: Failed to get latitude length"
+      stop
     else
       allocate(lat(latlength),stat=ierr)
@@ -334 +340 @@
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,lat)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading latitude"
+        write(*,*) "Error: Failed reading latitude"
+        stop
       endif
     endif
@@ -347 +354 @@
   ierr=NF_INQ_DIMID(infid,"lon",tmpdimid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get lon dimension ID"
+    write(*,*) "Error: Failed to get lon dimension ID"
+    stop
   else
     ierr=NF_INQ_VARID(infid,"lon",tmpvarid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get lon ID"
+      write(*,*) "Error: Failed to get lon ID"
+      stop
     else
       ierr=NF_INQ_DIMLEN(infid,tmpdimid,lonlength)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get lon length"
+        write(*,*) "Error: Failed to get lon length"
+        stop
       else
         allocate(lon(lonlength),stat=ierr)
@@ -365 +375 @@
         ierr=NF_GET_VAR_REAL(infid,tmpvarid,lon)
         if (ierr.ne.NF_NOERR) then
-          stop "Error: Failed reading longitude"
+          write(*,*) "Error: Failed reading longitude"
+          stop
         endif
       endif
@@ -373 +384 @@
   ierr=NF_INQ_VARID(infid,"longitude",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get longitude ID"
+    write(*,*) "Error: Failed to get longitude ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,lonlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get longitude length"
+      write(*,*) "Error: Failed to get longitude length"
+      stop
     else
       allocate(lon(lonlength),stat=ierr)
@@ -387 +400 @@
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,lon)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading longitude"
+        write(*,*) "Error: Failed reading longitude"
+        stop
       endif
     endif
@@ -404 +418 @@
     ierr=NF_INQ_DIMID(infid,"time_counter",tmpdimid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get time_counter dimension ID"
+      write(*,*) "Error: Failed to get time_counter dimension ID"
+      stop
     else
       ierr=NF_INQ_VARID(infid,"time_counter",tmpvarid)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get time_counter ID"
+        write(*,*) "Error: Failed to get time_counter ID"
+        stop
       else
         ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
         if (ierr.ne.NF_NOERR) then
-          stop "Error: Failed to get time_counter length"
+          write(*,*) "Error: Failed to get time_counter length"
+          stop
         else
           allocate(time(timelength),stat=ierr)
@@ -422 +439 @@
           ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
           if (ierr.ne.NF_NOERR) then
-            stop "Error: Failed reading time"
+            write(*,*) "Error: Failed reading time"
+            stop
           endif
         endif
@@ -430 +448 @@
     ierr=NF_INQ_VARID(infid,"time",tmpvarid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get time ID"
+      write(*,*) "Error: Failed to get time ID"
+      stop
     else
       ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get time length"
+        write(*,*) "Error: Failed to get time length"
+        stop
       else
         allocate(time(timelength),stat=ierr)
@@ -444 +464 @@
         ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
         if (ierr.ne.NF_NOERR) then
-          stop "Error: Failed reading time"
+          write(*,*) "Error: Failed reading time"
+          stop
         endif
       endif
@@ -452 +473 @@
   ierr=NF_INQ_VARID(infid,"Time",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get Time ID"
+    write(*,*) "Error: Failed to get Time ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get Time length"
+      write(*,*) "Error: Failed to get Time length"
+      stop
     else
       allocate(time(timelength),stat=ierr)
@@ -466 +489 @@
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading Time"
+        write(*,*) "Error: Failed reading Time"
+        stop
       endif
     endif
@@ -479 +503 @@
   ierr=NF_INQ_DIMID(infid,"sigma",tmpdimid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get sigma dimension ID"
+    write(*,*) "Error: Failed to get sigma dimension ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,altlength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get altitude length"
+      write(*,*) "Error: Failed to get altitude length"
+      stop
     endif
   endif
@@ -489 +515 @@
   ierr=NF_INQ_DIMLEN(infid,tmpdimid,altlength)
   if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get altitude length"
+      write(*,*) "Error: Failed to get altitude length"
+      stop
   endif
 endif
@@ -504 +531 @@
     ierr=NF_INQ_VARID(infid,"controle",tmpvarid)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get controle ID"
+      write(*,*) "Error: Failed to get controle ID"
+      stop
     else
       ierr=NF_INQ_DIMLEN(infid,tmpdimid,ctllength)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed to get controle length"
+        write(*,*) "Error: Failed to get controle length"
+        stop
       else
         allocate(ctl(ctllength),stat=ierr)
@@ -518 +547 @@
         ierr=NF_GET_VAR_REAL(infid,tmpvarid,ctl)
         if (ierr.ne.NF_NOERR) then
-          stop "Error: Failed reading controle"
+          write(*,*) "Error: Failed reading controle"
+          stop
         endif
       endif
@@ -526 +556 @@
   ierr=NF_INQ_VARID(infid,"controle",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get controle ID"
+    write(*,*) "Error: Failed to get controle ID"
+    stop
   else
     ierr=NF_INQ_DIMLEN(infid,tmpdimid,ctllength)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed to get controle length"
+      write(*,*) "Error: Failed to get controle length"
+      stop
     else
       allocate(ctl(ctllength),stat=ierr)
@@ -540 +572 @@
       ierr=NF_GET_VAR_REAL(infid,tmpvarid,ctl)
       if (ierr.ne.NF_NOERR) then
-        stop "Error: Failed reading controle"
+        write(*,*) "Error: Failed reading controle"
+        stop
       endif
     endif
@@ -563 +596 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,sigma)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading sigma"
+    write(*,*) "Error: Failed reading sigma"
+    stop
   endif
 endif
@@ -572 +606 @@
   ierr=NF_INQ_VARID(infid,"aps",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get aps ID"
+    write(*,*) "Error: Failed to get aps ID"
+    stop
   else
     allocate(aps(altlength),stat=ierr)
@@ -582 +617 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,aps)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading aps"
+      write(*,*) "Error: Failed reading aps"
+      stop
     endif
   endif
@@ -589 +625 @@
   ierr=NF_INQ_VARID(infid,"bps",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get bps ID"
+    write(*,*) "Error: Failed to get bps ID"
+    stop
   else
     allocate(bps(altlength),stat=ierr)
@@ -599 +636 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading bps"
+      write(*,*) "Error: Failed reading bps"
+      stop
     endif
   endif
@@ -631 +669 @@
       if (ierr.ne.NF_NOERR) then
         write(*,*) "Error: Could not open that file either"
-        stop "Might as well stop here"
+        write(*,*) "Might as well stop here"
+        stop
       endif
     endif
@@ -639 +678 @@
   ierr=NF_INQ_VARID(infid2,"phisinit",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get phisinit ID"
+    write(*,*) "Error: Failed to get phisinit ID"
+    stop
   endif
   ! Get physinit
   ierr=NF_GET_VAR_REAL(infid2,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading phisinit"
+    write(*,*) "Error: Failed reading phisinit"
+    stop
   endif
   ! Close file
@@ -652 +693 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading phisinit"
+    write(*,*) "Error: Failed reading phisinit"
+    stop
   endif
 endif
@@ -918 +960 @@
 ierr=NF_INQ_VARID(infid,"ps",tmpvarid)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Failed to get ps ID"
+  write(*,*) "Error: Failed to get ps ID"
+  stop
 else
   allocate(ps(lonlength,latlength,timelength),stat=ierr)
@@ -929 +972 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,ps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading surface pressure"
+    write(*,*) "Error: Failed reading surface pressure"
+    stop
   endif
 endif
@@ -1022 +1066 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,temp)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading atmospheric temperature"
+    write(*,*) "Error: Failed reading atmospheric temperature"
+    stop
   endif
   have_temp=.true.
@@ -1046 +1091 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,rho)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed reading atmospheric density"
+    write(*,*) "Error: Failed reading atmospheric density"
+    stop
   endif
 endif
@@ -1234 +1280 @@
 ierr=NF_PUT_ATT_TEXT(outfid,lon_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for longitude"
+  write(*,*) "Error: Problem writing long_name for longitude"
+  stop
 endif
 text='degrees_east'
 ierr=NF_PUT_ATT_TEXT(outfid,lon_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for longitude"
+  write(*,*) "Error: Problem writing units for longitude"
+  stop
 endif
 
@@ -1255 +1303 @@
 ierr=NF_PUT_ATT_TEXT(outfid,lat_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for latitude"
+  write(*,*) "Error: Problem writing long_name for latitude"
+  stop
 endif
 text='degrees_north'
 ierr=NF_PUT_ATT_TEXT(outfid,lat_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for latitude"
+  write(*,*) "Error: Problem writing units for latitude"
+  stop
 endif
 
@@ -1277 +1327 @@
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for altitude"
+    write(*,*) "Error: Problem writing long_name for altitude"
+    stop
   endif
   text='Pa'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for altitude"
+    write(*,*) "Error: Problem writing units for altitude"
+    stop
   endif
   text='down'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'positive',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing positive for altitude"
+    write(*,*) "Error: Problem writing positive for altitude"
+    stop
   endif
 else if (ztype.eq.2) then ! above areoid vertical coordinate
@@ -1293 +1346 @@
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for altitude"
+    write(*,*) "Error: Problem writing long_name for altitude"
+    stop
   endif
   text='m'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for altitude"
+    write(*,*) "Error: Problem writing units for altitude"
+    stop
   endif
   text='up'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'positive',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing positive for altitude"
+    write(*,*) "Error: Problem writing positive for altitude"
+    stop
   endif
 else if (ztype.eq.3) then ! above surface vertical coordinate
@@ -1309 +1365 @@
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for altitude"
+    write(*,*) "Error: Problem writing long_name for altitude"
+    stop
   endif
   text='m'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for altitude"
+    write(*,*) "Error: Problem writing units for altitude"
+    stop
   endif
   text='up'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'positive',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing positive for altitude"
+    write(*,*) "Error: Problem writing positive for altitude"
+    stop
   endif
 else ! radial distance to center of planet
@@ -1325 +1384 @@
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for altitude"
+    write(*,*) "Error: Problem writing long_name for altitude"
+    stop
   endif
   text='m'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for altitude"
+    write(*,*) "Error: Problem writing units for altitude"
+    stop
   endif
   text='up'
   ierr=NF_PUT_ATT_TEXT(outfid,alt_varid,'positive',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing positive for altitude"
+    write(*,*) "Error: Problem writing positive for altitude"
+    stop
   endif
 endif ! of if (have_sigma)
@@ -1352 +1414 @@
   ierr=NF_PUT_ATT_TEXT(outfid,ctl_varid,'title',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing title for controle"
+    write(*,*) "Error: Problem writing title for controle"
+    stop
   endif
 endif
@@ -1387 +1450 @@
   ierr=NF_PUT_ATT_TEXT(outfid,sigma_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for sigma"
+    write(*,*) "Error: Problem writing long_name for sigma"
+    stop
   endif
 else ! hybrid coordinates
@@ -1393 +1457 @@
   ierr=NF_PUT_ATT_TEXT(outfid,aps_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for aps"
+    write(*,*) "Error: Problem writing long_name for aps"
+    stop
   endif
   text="hybrid sigma at midlayers"
   ierr=NF_PUT_ATT_TEXT(outfid,bps_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for bps"
+    write(*,*) "Error: Problem writing long_name for bps"
+    stop
   endif
 endif ! of if (have_sigma)
@@ -1415 +1481 @@
 ierr=NF_PUT_ATT_TEXT(outfid,time_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for Time"
+  write(*,*) "Error: Problem writing long_name for Time"
+  stop
 endif
 text='days since 0000-01-1 00:00:00'
 ierr=NF_PUT_ATT_TEXT(outfid,time_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for Time"
+  write(*,*) "Error: Problem writing units for Time"
+  stop
 endif
 
@@ -1430 +1498 @@
 ierr=NF_PUT_ATT_TEXT(outfid,phis_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for phisinit"
+  write(*,*) "Error: Problem writing long_name for phisinit"
+  stop
 endif
 text='m2.s2'
 ierr=NF_PUT_ATT_TEXT(outfid,phis_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for geopotential phisinit"
+  write(*,*) "Error: Problem writing units for geopotential phisinit"
+  stop
 endif
 
@@ -1457 +1527 @@
 ierr=NF_PUT_ATT_TEXT(outfid,ps_varid,'long_name',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing long_name for surface pressure"
+  write(*,*) "Error: Problem writing long_name for surface pressure"
+  stop
 endif
 text='Pa'
 ierr=NF_PUT_ATT_TEXT(outfid,ps_varid,'units',len_trim(text),text)
 if (ierr.ne.NF_NOERR) then
-  stop "Error: Problem writing units for surface pressure"
+  write(*,*) "Error: Problem writing units for surface pressure"
+  stop
 endif
 
@@ -1480 +1552 @@
   ierr=NF_PUT_ATT_TEXT(outfid,za_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for zareoid"
+    write(*,*) "Error: Problem writing long_name for zareoid"
+    stop
   endif
   text='m'
   ierr=NF_PUT_ATT_TEXT(outfid,za_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for zareoid"
+    write(*,*) "Error: Problem writing units for zareoid"
+    stop
   endif
   ! zareoid missing value
   ierr=NF_PUT_ATT_REAL(outfid,za_varid,'missing_value',NF_REAL,1,miss_val)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing missing_value for zareoid"
+    write(*,*) "Error: Problem writing missing_value for zareoid"
+    stop
   endif
 else ! above areoid or above local surface vertical coordinate
@@ -1504 +1579 @@
   ierr=NF_PUT_ATT_TEXT(outfid,p_varid,'long_name',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing long_name for pressure"
+    write(*,*) "Error: Problem writing long_name for pressure"
+    stop
   endif
   text='Pa'
   ierr=NF_PUT_ATT_TEXT(outfid,p_varid,'units',len_trim(text),text)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing units for pressure"
+    write(*,*) "Error: Problem writing units for pressure"
+    stop
   endif
   ! pressure missing value
   ierr=NF_PUT_ATT_REAL(outfid,p_varid,'missing_value',NF_REAL,1,miss_val)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Problem writing missing_value for pressure"
+    write(*,*) "Error: Problem writing missing_value for pressure"
+    stop
   endif
 endif
@@ -1588 +1666 @@
       if (ierr.ne.NF_NOERR) then
         write(*,*) "Error failed to copy title attribute:",trim(text)
-      stop ""
+      stop
       endif
     else ! no 'title' attribute
@@ -1611 +1689 @@
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error failed to copy long_name attribute:",trim(text)
-      stop ""
+      stop
     endif
   endif
@@ -1625 +1703 @@
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error failed to copy units attribute:",trim(text)
-      stop ""
+      stop
     endif
   endif
@@ -1641 +1719 @@
   ierr=NF_PUT_ATT_REAL(outfid,var_id(i),'missing_value',NF_REAL,1,miss_val)
   if (ierr.ne.NF_NOERR) then
-    stop "Error, failed to write missing_value attribute"
+    write(*,*) "Error, failed to write missing_value attribute"
+    stop
   endif
 
@@ -1746 +1825 @@
 
 ! Write time
-ierr=NF_PUT_VARA_REAL(outfid,time_varid,1,timelength,time)
+ierr=NF_PUT_VARA_REAL(outfid,time_varid,(/1/),(/timelength/),time)
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Error: Could not write Time data to output file"