source: CONFIG/trunk/SCRIPT/REGRID_forcage/DYN2PHY/phygrid_ncfile.f @ 1356

      SUBROUTINE phygrid_ncfile (filename)
!--------------------------------------------------------------
!     Prepare a netCDF file for lmdz physics grid
!
!     Field units: in this file used by the GCM all units are
!     in kg-species/m2/s. The fields are modified accordingly.
!--------------------------------------------------------------
      USE sflx
      USE final_grid_lmdz
      USE grid, ONLY: nrecs

      IMPLICIT NONE

      INCLUDE 'netcdf.inc'

!--------------------------------------------------------------
!       ... Dummy arguments
!--------------------------------------------------------------
      CHARACTER(LEN=*), INTENT(in) :: filename

!--------------------------------------------------------------
!       ... Local variables
!--------------------------------------------------------------
      CHARACTER(LEN=100) :: name
      INTEGER :: len_name
      INTEGER :: iret
      INTEGER :: ncid
      INTEGER :: i

!     ... Dimension ids
      INTEGER :: vector_dim
      INTEGER :: time_dim
      INTEGER :: lev_dim

!     ... Dimension lengths
      INTEGER, PARAMETER :: time_len   =  NF_UNLIMITED
      INTEGER, PARAMETER :: month_len = 12

!     ... Variable ids
      INTEGER :: vector_id
      INTEGER :: time_id
      INTEGER :: lev_id
      INTEGER :: so4_id

      INTEGER, PARAMETER :: vector_rank = 1
      INTEGER, PARAMETER :: time_rank = 1
      INTEGER, PARAMETER :: lev_rank = 1

      INTEGER :: vector_dims(vector_rank)
      INTEGER :: time_dims(time_rank)
      INTEGER :: time_start(time_rank), time_count(time_rank)
      INTEGER :: lev_dims(lev_rank)

      REAL    :: month(month_len)
      REAL    :: vector (klon)
      REAL    :: field  (klon, month_len)

      INTEGER, PARAMETER :: so4_rank = 2
      INTEGER            :: so4_dims(so4_rank)

      DATA month /15.,45.,75.,105.,135.,165.,195.,225.,255.,285.,315.,345./

      DO i = 1, klon
        vector(i) = FLOAT(i)
      END DO

!     ... Enter define mode
      PRINT *, '     Enter netcdf define mode ...'
      iret = nf_create(filename, NF_CLOBBER, ncid)
      CALL check_err(iret)

!     ... Define dimensions
      iret = nf_def_dim(ncid, 'vector', klon, vector_dim)
      CALL check_err(iret)
      iret = nf_def_dim(ncid, 'time', NF_UNLIMITED, time_dim)
      CALL check_err(iret)
      iret = nf_def_dim(ncid, 'lev', nlev_len, lev_dim)
      CALL check_err(iret)

      vector_dims(1) = vector_dim
      iret = nf_def_var(ncid, 'vector', NF_FLOAT, vector_rank, vector_dims, vector_id)
      CALL check_err(iret)
      time_dims(1) = time_dim
      iret = nf_def_var(ncid, 'time', NF_FLOAT, time_rank, time_dims, time_id)
      CALL check_err(iret)
      lev_dims(1) = lev_dim
      iret = nf_def_var(ncid, 'lev', NF_FLOAT, lev_rank, lev_dims, lev_id)
      CALL check_err(iret)

      name = 'LMDz surface emissions'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, NF_GLOBAL, 'title', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'S. Szopa'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, NF_GLOBAL, 'source', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'physics grid point'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, vector_id, 'long_name', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'index'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, vector_id, 'units', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'days since 0000-01-01'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, time_id, 'units', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'Time'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, time_id, 'long_name', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)

!     ... Fields
      so4_dims(3) = time_dim
      so4_dims(2) = lev_dim
      so4_dims(1) = vector_dim

      iret = nf_def_var(ncid, 'so4', NF_FLOAT, so4_rank, so4_dims, so4_id)
      CALL check_err(iret)
      name = 'vmr'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, so4_id, 'units', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)
      name = 'SO4'
      len_name = LEN_TRIM(ADJUSTL(name))
      iret = nf_put_att_text(ncid, so4_id, 'long_name', len_name, name(:LEN_TRIM(name)))
      CALL check_err(iret)


!     ... Leave define mode
      PRINT *, '     Leave netcdf define mode.'
      iret = nf_enddef(ncid)
      CALL check_err(iret)

      PRINT *, '     Store variables ...'
!     ... Store time
      time_start(1) = 1
      time_count(1) = month_len
      iret = nf_put_vara_real(ncid, time_id, time_start, time_count, month)
      CALL check_err(iret)

!     ... store vector
      iret = nf_put_var_real(ncid, vector_id, vector)
      CALL check_err(iret)
!     ... store lev
      iret = nf_put_var_real(ncid, lev_id, nlev)
      CALL check_err(iret)

!     ... store fields


      iret = nf_put_var_real(ncid, so4_id, so4_phy_1x1)
      CALL check_err(iret)


      PRINT *, '     Done phy file.'
      iret = nf_close(ncid)
      CALL check_err(iret)

      END SUBROUTINE phygrid_ncfile