Changeset 1910 for LMDZ5/branches/testing/libf/dyn3d

Timestamp:

Nov 29, 2013, 9:40:25 AM (11 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r1860:1909 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (2 props)
• libf (modified) (1 prop)
• libf/dyn3d (modified) (1 prop)
• libf/dyn3d/PVtheta.F (modified) (1 prop)
• libf/dyn3d/abort_gcm.F (modified) (1 prop)
• libf/dyn3d/academic.h (modified) (1 prop)
• libf/dyn3d/adaptdt.F (modified) (1 prop)
• libf/dyn3d/addfi.F (modified) (1 prop)
• libf/dyn3d/advect.F (modified) (1 prop)
• libf/dyn3d/advn.F (modified) (1 prop)
• libf/dyn3d/advtrac.F90 (modified) (1 prop)
• libf/dyn3d/advx.F (modified) (1 prop)
• libf/dyn3d/advxp.F (modified) (1 prop)
• libf/dyn3d/advy.F (modified) (1 prop)
• libf/dyn3d/advyp.F (modified) (1 prop)
• libf/dyn3d/advz.F (modified) (1 prop)
• libf/dyn3d/advzp.F (modified) (1 prop)
• libf/dyn3d/bernoui.F (modified) (1 prop)
• libf/dyn3d/bilan_dyn.F (modified) (1 prop)
• libf/dyn3d/caladvtrac.F (modified) (1 prop)
• libf/dyn3d/caldyn.F (modified) (1 prop)
• libf/dyn3d/caldyn0.F (modified) (1 prop)
• libf/dyn3d/calfis.F (modified) (1 prop)
• libf/dyn3d/ce0l.F90 (modified) (1 prop)
• libf/dyn3d/clesph0.h (modified) (1 prop)
• libf/dyn3d/coefpoly.F (modified) (1 prop)
• libf/dyn3d/com_io_dyn_mod.F90 (modified) (1 prop)
• libf/dyn3d/comconst.h (modified) (1 prop)
• libf/dyn3d/comdissip.h (modified) (1 prop)
• libf/dyn3d/comdissipn.h (modified) (1 prop)
• libf/dyn3d/comdissnew.h (modified) (1 prop)
• libf/dyn3d/comgeom.h (modified) (1 prop)
• libf/dyn3d/comgeom2.h (modified) (1 prop)
• libf/dyn3d/comvert.h (modified) (1 prop)
• libf/dyn3d/conf_dat2d.F (modified) (1 prop)
• libf/dyn3d/conf_dat3d.F (modified) (1 prop)
• libf/dyn3d/conf_gcm.F (modified) (1 prop)
• libf/dyn3d/conf_planete.F90 (modified) (1 prop)
• libf/dyn3d/control_mod.F90 (modified) (1 prop)
• libf/dyn3d/convflu.F (modified) (1 prop)
• libf/dyn3d/convmas.F (modified) (1 prop)
• libf/dyn3d/coordij.F (modified) (1 prop)
• libf/dyn3d/covcont.F (modified) (1 prop)
• libf/dyn3d/covnat.F (modified) (1 prop)
• libf/dyn3d/defrun.F (modified) (1 prop)
• libf/dyn3d/description.h (modified) (1 prop)
• libf/dyn3d/diagedyn.F (modified) (1 prop)
• libf/dyn3d/dissip.F (modified) (1 prop)
• libf/dyn3d/disvert.F90 (modified) (1 prop)
• libf/dyn3d/disvert_noterre.F (modified) (1 prop)
• libf/dyn3d/diverg.F (modified) (1 prop)
• libf/dyn3d/diverg_gam.F (modified) (1 prop)
• libf/dyn3d/divergf.F (modified) (1 prop)
• libf/dyn3d/divergst.F (modified) (1 prop)
• libf/dyn3d/divgrad.F (modified) (1 prop)
• libf/dyn3d/divgrad2.F (modified) (1 prop)
• libf/dyn3d/dteta1.F (modified) (1 prop)
• libf/dyn3d/dudv1.F (modified) (1 prop)
• libf/dyn3d/dudv2.F (modified) (1 prop)
• libf/dyn3d/dump2d.F (modified) (1 prop)
• libf/dyn3d/dynetat0.F (modified) (1 prop)
• libf/dyn3d/dynredem.F (modified) (1 prop)
• libf/dyn3d/ener.h (modified) (1 prop)
• libf/dyn3d/enercin.F (modified) (1 prop)
• libf/dyn3d/exner_hyb.F (modified) (1 prop)
• libf/dyn3d/exner_milieu.F (modified) (1 prop)
• libf/dyn3d/extrapol.F (modified) (1 prop)
• libf/dyn3d/flumass.F (modified) (1 prop)
• libf/dyn3d/fluxstokenc.F (modified) (1 prop)
• libf/dyn3d/friction.F (modified) (1 prop)
• libf/dyn3d/fxhyp.F (modified) (1 prop)
• libf/dyn3d/fxy.F (modified) (1 prop)
• libf/dyn3d/fxyhyper.F (modified) (1 prop)
• libf/dyn3d/fxysinus.F (modified) (1 prop)
• libf/dyn3d/fyhyp.F (modified) (1 prop)
• libf/dyn3d/gcm.F (modified) (1 prop)
• libf/dyn3d/geopot.F (modified) (1 prop)
• libf/dyn3d/getparam.F90 (modified) (1 prop)
• libf/dyn3d/gr_dyn_fi.F (modified) (1 prop)
• libf/dyn3d/gr_ecrit_fi.F (modified) (1 prop)
• libf/dyn3d/gr_fi_dyn.F (modified) (1 prop)
• libf/dyn3d/gr_int_dyn.F (modified) (1 prop)
• libf/dyn3d/gr_u_scal.F (modified) (1 prop)
• libf/dyn3d/gr_v_scal.F (modified) (1 prop)
• libf/dyn3d/grad.F (modified) (1 prop)
• libf/dyn3d/gradiv.F (modified) (1 prop)
• libf/dyn3d/gradiv2.F (modified) (1 prop)
• libf/dyn3d/gradsdef.h (modified) (1 prop)
• libf/dyn3d/grid_atob.F (modified) (1 prop)
• libf/dyn3d/grid_noro.F (modified) (1 prop)
• libf/dyn3d/grilles_gcm_netcdf.F (modified) (1 prop)
• libf/dyn3d/grilles_gcm_netcdf_sub.F90 (modified) (1 prop)
• libf/dyn3d/groupe.F (modified) (1 prop)
• libf/dyn3d/groupeun.F (modified) (1 prop)
• libf/dyn3d/guide_mod.F90 (modified) (15 diffs, 1 prop)
• libf/dyn3d/heavyside.F (modified) (1 prop)
• libf/dyn3d/infotrac.F90 (modified) (1 prop)
• libf/dyn3d/ini_paramLMDZ_dyn.h (modified) (1 prop)
• libf/dyn3d/iniacademic.F90 (modified) (1 prop)
• libf/dyn3d/iniconst.F90 (modified) (1 prop)
• libf/dyn3d/inidissip.F90 (modified) (1 prop)
• libf/dyn3d/inigeom.F (modified) (1 prop)
• libf/dyn3d/inigrads.F (modified) (1 diff, 1 prop)
• libf/dyn3d/iniinterp_horiz.F (modified) (1 prop)
• libf/dyn3d/iniprint.h (modified) (1 prop)
• libf/dyn3d/initial0.F (modified) (1 prop)
• libf/dyn3d/integrd.F (modified) (1 prop)
• libf/dyn3d/inter_barxy_m.F90 (modified) (1 prop)
• libf/dyn3d/interp_horiz.F (modified) (1 prop)
• libf/dyn3d/interpost.F (modified) (1 prop)
• libf/dyn3d/interpre.F (modified) (1 prop)
• libf/dyn3d/invert_lat.F90 (modified) (1 prop)
• libf/dyn3d/juldate.F (modified) (1 prop)
• libf/dyn3d/laplacien.F (modified) (1 prop)
• libf/dyn3d/laplacien_gam.F (modified) (1 prop)
• libf/dyn3d/laplacien_rot.F (modified) (1 prop)
• libf/dyn3d/laplacien_rotgam.F (modified) (1 prop)
• libf/dyn3d/leapfrog.F (modified) (1 prop)
• libf/dyn3d/limx.F (modified) (1 prop)
• libf/dyn3d/limy.F (modified) (1 prop)
• libf/dyn3d/limz.F (modified) (1 prop)
• libf/dyn3d/logic.h (modified) (1 prop)
• libf/dyn3d/massbar.F (modified) (1 prop)
• libf/dyn3d/massbarxy.F (modified) (1 prop)
• libf/dyn3d/massdair.F (modified) (1 prop)
• libf/dyn3d/minmax.F (modified) (1 prop)
• libf/dyn3d/minmax2.F (modified) (1 prop)
• libf/dyn3d/mod_const_mpi.F90 (modified) (1 prop)
• libf/dyn3d/nxgrad.F (modified) (1 prop)
• libf/dyn3d/nxgrad_gam.F (modified) (1 prop)
• libf/dyn3d/nxgradst.F (modified) (1 prop)
• libf/dyn3d/nxgraro2.F (modified) (1 prop)
• libf/dyn3d/nxgrarot.F (modified) (1 prop)
• libf/dyn3d/paramet.h (modified) (1 prop)
• libf/dyn3d/pbar.F (modified) (1 prop)
• libf/dyn3d/pentes_ini.F (modified) (1 prop)
• libf/dyn3d/ppm3d.F (modified) (1 prop)
• libf/dyn3d/prather.F (modified) (1 prop)
• libf/dyn3d/pres2lev_mod.F90 (modified) (1 prop)
• libf/dyn3d/pression.F (modified) (1 prop)
• libf/dyn3d/profvert.def (modified) (1 prop)
• libf/dyn3d/psextbar.F (modified) (1 prop)
• libf/dyn3d/q_sat.F (modified) (1 prop)
• libf/dyn3d/qminimum.F (modified) (1 prop)
• libf/dyn3d/ran1.F (modified) (1 prop)
• libf/dyn3d/rotat.F (modified) (1 prop)
• libf/dyn3d/rotat_nfil.F (modified) (1 prop)
• libf/dyn3d/rotatf.F (modified) (1 prop)
• libf/dyn3d/rotatst.F (modified) (1 prop)
• libf/dyn3d/serre.h (modified) (1 prop)
• libf/dyn3d/sort.F (modified) (1 prop)
• libf/dyn3d/sortvarc.F (modified) (1 prop)
• libf/dyn3d/sortvarc0.F (modified) (1 prop)
• libf/dyn3d/startvar.F90 (modified) (1 prop)
• libf/dyn3d/sw_case_williamson91_6.F (modified) (1 prop)
• libf/dyn3d/temps.h (modified) (1 prop)
• libf/dyn3d/test_period.F (modified) (1 prop)
• libf/dyn3d/tetaleveli1j.F (modified) (1 prop)
• libf/dyn3d/tetaleveli1j1.F (modified) (1 prop)
• libf/dyn3d/top_bound.F (modified) (1 prop)
• libf/dyn3d/tourabs.F (modified) (1 prop)
• libf/dyn3d/tourpot.F (modified) (1 prop)
• libf/dyn3d/traceurpole.F (modified) (1 prop)
• libf/dyn3d/tracstoke.h (modified) (1 prop)
• libf/dyn3d/ugeostr.F90 (modified) (1 prop)
• libf/dyn3d/vitvert.F (modified) (1 prop)
• libf/dyn3d/vlsplt.F (modified) (1 prop)
• libf/dyn3d/vlspltqs.F (modified) (1 prop)
• libf/dyn3d/wrgrads.F (modified) (1 prop)
• libf/dyn3d/write_grads_dyn.h (modified) (1 prop)
• libf/dyn3d/write_paramLMDZ_dyn.h (modified) (1 prop)

LMDZ5/branches/testing/libf/dyn3d/guide_mod.F90

@@ -143 +143 @@
     ncidpl=-99
     if (guide_modele) then
-       if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
+       if (ncidpl.eq.-99) then
+          rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
+          if (rcod.NE.NF_NOERR) THEN
+             print *,'Guide: probleme -> pas de fichier apbp.nc'
+             CALL abort_gcm(modname,abort_message,1)
+          endif
+       endif
     else
          if (guide_u) then
-           if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
+           if (ncidpl.eq.-99) then
+               rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
+               if (rcod.NE.NF_NOERR) THEN
+                  print *,'Guide: probleme -> pas de fichier u.nc'
+                  CALL abort_gcm(modname,abort_message,1)
+               endif
+           endif
          elseif (guide_v) then
-           if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
+           if (ncidpl.eq.-99) then
+               rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
+               if (rcod.NE.NF_NOERR) THEN
+                  print *,'Guide: probleme -> pas de fichier v.nc'
+                  CALL abort_gcm(modname,abort_message,1)
+               endif
+           endif
          elseif (guide_T) then
-           if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
+           if (ncidpl.eq.-99) then
+               rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
+               if (rcod.NE.NF_NOERR) THEN
+                  print *,'Guide: probleme -> pas de fichier T.nc'
+                  CALL abort_gcm(modname,abort_message,1)
+               endif
+           endif
          elseif (guide_Q) then
-           if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
+           if (ncidpl.eq.-99) then
+               rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
+               if (rcod.NE.NF_NOERR) THEN
+                  print *,'Guide: probleme -> pas de fichier hur.nc'
+                  CALL abort_gcm(modname,abort_message,1)
+               endif
+           endif
          endif
     endif 
@@ -990 +1020 @@
     INTEGER               :: status,rcode
 
+    CHARACTER (len = 80)   :: abort_message
+    CHARACTER (len = 20)   :: modname = 'guide_read'
 ! -----------------------------------------------------------------
 ! Premier appel: initialisation de la lecture des fichiers
@@ -998 +1030 @@
 ! Niveaux de pression si non constants
          if (guide_modele) then
-             print *,'Lecture du guidage sur niveaux mod�le'
+             print *,'Lecture du guidage sur niveaux modele'
              rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable AP, fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable BP, fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidpl,varidap',ncidpl,varidap
          endif
@@ -1007 +1051 @@
          if (guide_u) then
              rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier u.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable UWND, fichier u.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidu,varidu',ncidu,varidu
              if (ncidpl.eq.-99) ncidpl=ncidu
@@ -1014 +1066 @@
          if (guide_v) then
              rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier v.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable VWND, fichier v.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidv,varidv',ncidv,varidv
              if (ncidpl.eq.-99) ncidpl=ncidv
@@ -1021 +1081 @@
          if (guide_T) then
              rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier T.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable AIR, fichier T.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidT,varidT',ncidt,varidt
              if (ncidpl.eq.-99) ncidpl=ncidt
@@ -1028 +1096 @@
          if (guide_Q) then
              rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier hur.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable RH, fichier hur.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidQ,varidQ',ncidQ,varidQ
              if (ncidpl.eq.-99) ncidpl=ncidQ
@@ -1035 +1111 @@
          if ((guide_P).OR.(guide_modele)) then
              rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier ps.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidps, 'SP', varidps)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable SP, fichier ps.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidps,varidps',ncidps,varidps
          endif
@@ -1173 +1257 @@
     INTEGER               :: i
 
+    CHARACTER (len = 80)   :: abort_message
+    CHARACTER (len = 20)   :: modname = 'guide_read2D'
 ! -----------------------------------------------------------------
 ! Premier appel: initialisation de la lecture des fichiers
@@ -1181 +1267 @@
 ! Niveaux de pression si non constants
          if (guide_modele) then
-             print *,'Lecture du guidage sur niveaux mod�le'
+             print *,'Lecture du guidage sur niveaux modele'
              rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable AP, fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable BP, fichier apbp.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidpl,varidap',ncidpl,varidap
          endif
@@ -1190 +1288 @@
          if (guide_u) then
              rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier u.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable UWND, fichier u.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidu,varidu',ncidu,varidu
              if (ncidpl.eq.-99) ncidpl=ncidu
@@ -1197 +1303 @@
          if (guide_v) then
              rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier v.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable VWND, fichier v.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidv,varidv',ncidv,varidv
              if (ncidpl.eq.-99) ncidpl=ncidv
@@ -1204 +1318 @@
          if (guide_T) then
              rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier T.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable AIR, fichier T.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidT,varidT',ncidt,varidt
              if (ncidpl.eq.-99) ncidpl=ncidt
@@ -1211 +1333 @@
          if (guide_Q) then
              rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier hur.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable RH, fichier hur.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidQ,varidQ',ncidQ,varidQ
              if (ncidpl.eq.-99) ncidpl=ncidQ
@@ -1218 +1348 @@
          if ((guide_P).OR.(guide_modele)) then
              rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de fichier ps.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              rcode = nf90_inq_varid(ncidps, 'SP', varidps)
+             IF (rcode.NE.NF_NOERR) THEN
+              print *,'Guide: probleme -> pas de variable SP, fichier ps.nc'
+              CALL abort_gcm(modname,abort_message,1)
+             ENDIF
              print*,'ncidps,varidps',ncidps,varidps
          endif

LMDZ5/branches/testing/libf/dyn3d/inigrads.F

@@ -13 +13 @@
       real xmin,xmax,ymin,ymax
 
-      character file*10,titlel*40
+      character(len=*),intent(in) :: file
+      character(len=*),intent(in) :: titlel
 
 #include "gradsdef.h"