source: LMDZ6/branches/SETHET_DECOUPLE/libf/phymar/iotd_ini.F90

      SUBROUTINE iotd_ini(fichnom,iim,jjm,llm,prlonv,prlatu,pcoordv)
      IMPLICIT NONE

!=======================================================================
!
!   Auteur:  L. Fairhead  ,  P. Le Van, Y. Wanherdrick, F. Forget
!   -------
!
!   Objet:
!   ------
!
!   'Initialize' the diagfi.nc file: write down dimensions as well
!   as time-independent fields (e.g: geopotential, mesh area, ...)
!
!=======================================================================
!-----------------------------------------------------------------------
!   Declarations:
!   -------------

#include "netcdf.inc"
#include "iotd.h"

!   Arguments:
!   ----------

      integer iim,jjm,llm
      real prlonv(iim),prlatu(jjm),pcoordv(llm),timestep
      INTEGER id_FOCE

      integer corner(4),edges(4),ndim
      real  px(1000)
      character (len=10) :: nom

!   Local:
!   ------
      INTEGER ierr

      integer :: nvarid
      integer, dimension(2) :: id  
      integer :: varid

      character*10 fichnom
      real*4 rlonv(iim),rlatu(jjm),coordv(llm)

      real pi

      print*,'INIIO prlonv ',prlonv
      imax=iim
      jmax=jjm
      lmax=llm

      rlonv=prlonv
      rlatu=prlatu
      coordv=pcoordv

!-----------------------------------------------------------------------


      pi=2.*asin(1.)

! Define dimensions
    
         ! Create the NetCDF file
         ierr=NF_CREATE(fichnom, NF_CLOBBER, nid)
         ! Define the 'Time' dimension
         ierr=nf_def_dim(nid,"Time",NF_UNLIMITED,dim_coord(4))
         ! Define the 'Time' variable
         ierr=NF_DEF_VAR(nid, "Time", NF_FLOAT, 1, dim_coord(4),varid)
!        ! Add a long_name attribute
!        ierr=NF_PUT_ATT_TEXT(nid, varid, "long_name",4,"Time")
!        ! Add a units attribute
         ierr=NF_PUT_ATT_TEXT(nid, varid,'units',29,"days since 0000-00-0 00:00:00")
         ! Switch out of NetCDF Define mode

      ierr=NF_DEF_DIM(nid, "longitude", iim, dim_coord(1))
      ierr=NF_DEF_DIM(nid, "latitude", jjm, dim_coord(2))
      ierr=NF_DEF_DIM(nid, "altitude", llm, dim_coord(3))


      ierr=NF_ENDDEF(nid)
!
!  Contol parameters for this run
! --------------------------

      ierr=NF_REDEF(nid)
      ierr=NF_DEF_VAR(nid,"longitude", NF_FLOAT, 1, dim_coord(1),nvarid)
!     ierr=NF_PUT_ATT_TEXT(nid,nvarid,"long_name", 14,
!    .      "East longitude")
!     ierr=NF_PUT_ATT_TEXT(nid,nvarid,'units',12,"degrees_east")
      ierr=NF_ENDDEF(nid)
      ierr=NF_PUT_VAR_REAL(nid,nvarid,rlonv)
       print*,ierr

! --------------------------
      ierr=NF_REDEF(nid)
      ierr=NF_DEF_VAR(nid, "latitude", NF_FLOAT, 1, dim_coord(2),nvarid)
!     ierr=NF_PUT_ATT_TEXT(nid,nvarid,'units',13,"degrees_north")
!     ierr=NF_PUT_ATT_TEXT(nid,nvarid,"long_name", 14,"North latitude")
      ierr=NF_ENDDEF(nid)
      ierr=NF_PUT_VAR_REAL(nid,nvarid,rlatu)
!
! --------------------------
      ierr=NF_REDEF(nid)
      ierr=NF_DEF_VAR(nid, "altitude", NF_FLOAT, 1,dim_coord(3),nvarid)
      ierr=NF_PUT_ATT_TEXT(nid,nvarid,"long_name",10,"pseudo-alt")
!     ierr=NF_PUT_ATT_TEXT(nid,nvarid,'units',2,"km")
      if ( pcoordv(2)>pcoordv(1) ) then
         ierr=NF_PUT_ATT_TEXT(nid,nvarid,"long_name",10,"pseudo-alt")
         ierr=NF_PUT_ATT_TEXT(nid,nvarid,'positive',2,"up")
      else
         ierr=NF_PUT_ATT_TEXT(nid,nvarid,"long_name",8,"pressure")
         ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',4,"down")
      endif
      ierr=NF_ENDDEF(nid)

      ierr=NF_PUT_VAR_REAL(nid,nvarid,coordv)
!
      END