Changeset 1992 for LMDZ5/trunk/libf/phylmd/inistats.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (11 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/inistats.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/inistats.F) (1 diff)

LMDZ5/trunk/libf/phylmd/inistats.F90

@@ -1 @@
-      subroutine inistats(ierr)
+SUBROUTINE inistats(ierr)
 
-      implicit none
+  IMPLICIT NONE
 
-#include "dimensions.h"
-#include "paramet.h"
-#include "comgeom.h"
-#include "comvert.h"
-#include "comconst.h"
-#include "statto.h"
-#include "netcdf.inc"
+  include "dimensions.h"
+  include "paramet.h"
+  include "comgeom.h"
+  include "comvert.h"
+  include "comconst.h"
+  include "statto.h"
+  include "netcdf.inc"
 
-      integer,intent(out) :: ierr
-      integer :: nid
-      integer :: l,nsteppd
-      real, dimension(llm) ::  sig_s
-      integer :: idim_lat,idim_lon,idim_llm,idim_llmp1,idim_time
-      real, dimension(istime) :: lt
-      integer :: nvarid
+  INTEGER, INTENT (OUT) :: ierr
+  INTEGER :: nid
+  INTEGER :: l, nsteppd
+  REAL, DIMENSION (llm) :: sig_s
+  INTEGER :: idim_lat, idim_lon, idim_llm, idim_llmp1, idim_time
+  REAL, DIMENSION (istime) :: lt
+  INTEGER :: nvarid
 
-      write (*,*) 
-      write (*,*) '                        || STATS ||'
-      write (*,*) 
-      write (*,*) 'daysec',daysec
-      write (*,*) 'dtphys',dtphys
-      nsteppd=nint(daysec/dtphys)
-      write (*,*) 'nsteppd=',nsteppd
-      if (abs(float(nsteppd)-daysec/dtphys).gt.1.e-8*daysec)
-     &   stop'Dans Instat:  1jour .ne. n pas physiques'
+  WRITE (*, *)
+  WRITE (*, *) '                        || STATS ||'
+  WRITE (*, *)
+  WRITE (*, *) 'daysec', daysec
+  WRITE (*, *) 'dtphys', dtphys
+  nsteppd = nint(daysec/dtphys)
+  WRITE (*, *) 'nsteppd=', nsteppd
+  IF (abs(float(nsteppd)-daysec/dtphys)>1.E-8*daysec) &
+    STOP 'Dans Instat:  1jour .ne. n pas physiques'
 
-      if(mod(nsteppd,istime).ne.0)
-     &   stop'Dans Instat:  1jour .ne. n*istime pas physiques'
+  IF (mod(nsteppd,istime)/=0) STOP &
+    'Dans Instat:  1jour .ne. n*istime pas physiques'
 
-      istats=nsteppd/istime
-      write (*,*) 'istats=',istats
-      write (*,*) 'Storing ',istime,'times per day'
-      write (*,*) 'thus every ',istats,'physical timestep '
-      write (*,*) 
+  istats = nsteppd/istime
+  WRITE (*, *) 'istats=', istats
+  WRITE (*, *) 'Storing ', istime, 'times per day'
+  WRITE (*, *) 'thus every ', istats, 'physical timestep '
+  WRITE (*, *)
 
-      do l= 1, llm
-         sig_s(l)=((ap(l)+ap(l+1))/preff+bp(l)+bp(l+1))/2.
-         pseudoalt(l)=log(preff/presnivs(l))*8.  
-      enddo
+  DO l = 1, llm
+    sig_s(l) = ((ap(l)+ap(l+1))/preff+bp(l)+bp(l+1))/2.
+    pseudoalt(l) = log(preff/presnivs(l))*8.
+  END DO
 
-      ierr = NF_CREATE("stats.nc",NF_CLOBBER,nid)
-      if (ierr.ne.NF_NOERR) then
-         write (*,*) NF_STRERROR(ierr)
-         stop ""
-      endif
+  ierr = nf_create('stats.nc', nf_clobber, nid)
+  IF (ierr/=nf_noerr) THEN
+    WRITE (*, *) nf_strerror(ierr)
+    STOP ''
+  END IF
 
-      ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_lat)
-      ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_lon)
-      ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
-      ierr = NF_DEF_DIM (nid, "llmp1", llm+1, idim_llmp1)
-      ierr = NF_DEF_DIM (nid, "Time", NF_UNLIMITED, idim_time)
+  ierr = nf_def_dim(nid, 'latitude', jjp1, idim_lat)
+  ierr = nf_def_dim(nid, 'longitude', iip1, idim_lon)
+  ierr = nf_def_dim(nid, 'altitude', llm, idim_llm)
+  ierr = nf_def_dim(nid, 'llmp1', llm+1, idim_llmp1)
+  ierr = nf_def_dim(nid, 'Time', nf_unlimited, idim_time)
 
-      ierr = NF_ENDDEF(nid)
-      call def_var_stats(nid,"Time","Time",
-     &            "hours since 0000-00-0 00:00:00",1,
-     &            idim_time,nvarid,ierr)
-! Time is initialised later by mkstats subroutine
+  ierr = nf_enddef(nid)
+  CALL def_var_stats(nid, 'Time', 'Time', 'hours since 0000-00-0 00:00:00', &
+    1, idim_time, nvarid, ierr)
+  ! Time is initialised later by mkstats subroutine
 
-      call def_var_stats(nid,"latitude","latitude",
-     &            "degrees_north",1,idim_lat,nvarid,ierr)
+  CALL def_var_stats(nid, 'latitude', 'latitude', 'degrees_north', 1, &
+    idim_lat, nvarid, ierr)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu/pi*180)
+  ierr = nf_put_var_double(nid, nvarid, rlatu/pi*180)
 #else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu/pi*180)
+  ierr = nf_put_var_real(nid, nvarid, rlatu/pi*180)
 #endif
-      call def_var_stats(nid,"longitude","East longitude",
-     &            "degrees_east",1,idim_lon,nvarid,ierr)
+  CALL def_var_stats(nid, 'longitude', 'East longitude', 'degrees_east', 1, &
+    idim_lon, nvarid, ierr)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv/pi*180)
+  ierr = nf_put_var_double(nid, nvarid, rlonv/pi*180)
 #else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv/pi*180)
+  ierr = nf_put_var_real(nid, nvarid, rlonv/pi*180)
 #endif
 
-! Niveaux verticaux, aps et bps
-      ierr = NF_REDEF (nid)
-! presnivs
+  ! Niveaux verticaux, aps et bps
+  ierr = nf_redef(nid)
+  ! presnivs
 #ifdef NC_DOUBLE
-      ierr = NF_DEF_VAR (nid,"presnivs", NF_DOUBLE, 1,idim_llm,nvarid)
+  ierr = nf_def_var(nid, 'presnivs', nf_double, 1, idim_llm, nvarid)
 #else
-      ierr = NF_DEF_VAR (nid,"presnivs", NF_FLOAT, 1,idim_llm,nvarid)
+  ierr = nf_def_var(nid, 'presnivs', nf_float, 1, idim_llm, nvarid)
 #endif
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",15,
-     &                        "Vertical levels")
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"Pa")
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',4,"down")
-      ierr = NF_ENDDEF(nid)
+  ierr = nf_put_att_text(nid, nvarid, 'long_name', 15, 'Vertical levels')
+  ierr = nf_put_att_text(nid, nvarid, 'units', 2, 'Pa')
+  ierr = nf_put_att_text(nid, nvarid, 'positive', 4, 'down')
+  ierr = nf_enddef(nid)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid, 
-     &                          presnivs(1:llm))
+  ierr = nf_put_var_double(nid, nvarid, presnivs(1:llm))
 #else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid, 
-     &                        presnivs(1:llm))
-#endif 
-! Pseudo alts
+  ierr = nf_put_var_real(nid, nvarid, presnivs(1:llm))
+#endif
+  ! Pseudo alts
 #ifdef NC_DOUBLE
-      ierr = NF_DEF_VAR (nid,"altitude", NF_DOUBLE, 1,idim_llm,nvarid)
+  ierr = nf_def_var(nid, 'altitude', nf_double, 1, idim_llm, nvarid)
 #else
-      ierr = NF_DEF_VAR (nid,"altitude", NF_FLOAT, 1,idim_llm,nvarid)
+  ierr = nf_def_var(nid, 'altitude', nf_float, 1, idim_llm, nvarid)
 #endif
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",8,"altitude")
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
-      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
-      ierr = NF_ENDDEF(nid)
+  ierr = nf_put_att_text(nid, nvarid, 'long_name', 8, 'altitude')
+  ierr = nf_put_att_text(nid, nvarid, 'units', 2, 'km')
+  ierr = nf_put_att_text(nid, nvarid, 'positive', 2, 'up')
+  ierr = nf_enddef(nid)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
+  ierr = nf_put_var_double(nid, nvarid, pseudoalt)
 #else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
-#endif 
-!      call def_var_stats(nid,"aps","hybrid pressure at midlayers"," ", 
-!     &            1,idim_llm,nvarid,ierr)
-!#ifdef NC_DOUBLE
-!      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
-!#else
-!      ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
-!#endif
+  ierr = nf_put_var_real(nid, nvarid, pseudoalt)
+#endif
+  ! call def_var_stats(nid,"aps","hybrid pressure at midlayers"," ",
+  ! &            1,idim_llm,nvarid,ierr)
+  ! #ifdef NC_DOUBLE
+  ! ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
+  ! #else
+  ! ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
+  ! #endif
 
-!      call def_var_stats(nid,"bps","hybrid sigma at midlayers"," ", 
-!     &            1,idim_llm,nvarid,ierr)
-!#ifdef NC_DOUBLE
-!      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
-!#else 
-!      ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
-!#endif
+  ! call def_var_stats(nid,"bps","hybrid sigma at midlayers"," ",
+  ! &            1,idim_llm,nvarid,ierr)
+  ! #ifdef NC_DOUBLE
+  ! ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
+  ! #else
+  ! ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
+  ! #endif
 
-      ierr=NF_CLOSE(nid)
+  ierr = nf_close(nid)
 
-      end subroutine inistats
+END SUBROUTINE inistats