program convert_lai  implicit none  include 'netcdf.inc'!-----------------------------------------------------------------! make lai-sai-hgt dataset!-----------------------------------------------------------------  integer, parameter :: r8 = selected_real_kind(12)! File specific settings  integer, parameter :: nlon = 720   !input grid : longitude points  integer, parameter :: nlat = 360   !input grid : latitude  points  integer, parameter :: numpft = 16  !number of plant types  real(r8) landmask(nlon,nlat)       !fraction of land  real(r8) mlai (nlon,nlat,0:numpft) !monthly lai in  real(r8) msai (nlon,nlat,0:numpft) !monthly sai in  real(r8) mhgtt(nlon,nlat,0:numpft) !monthly height (top) in  real(r8) mhgtb(nlon,nlat,0:numpft) !monthly height (bottom) in  real(r8) readdum  integer :: dimlon_id               !netCDF dimension id  integer :: dimlat_id               !netCDF dimension id  integer :: dimpft_id               !netCDF dimension id  integer :: dimtim_id               !netCDF dimension id  integer :: lon_id                  !1d longitude array id  integer :: lat_id                  !1d latitude array id  integer :: longxy_id               !2d longitude array id  integer :: latixy_id               !2d latitude array id  integer :: edgen_id                !northern edge of grid (edge(1)) id  integer :: edgee_id                !eastern  edge of grid (edge(2)) id  integer :: edges_id                !southern edge of grid (edge(3)) id  integer :: edgew_id                !western  edge of grid (edge(4)) id  integer :: landmask_id             !landmask id  integer :: mlai_id                 !monthly mlai id  integer :: msai_id                 !monthly msai id    integer :: mhgtt_id                !monthly mghtt id  integer :: mhgtb_id                !monthly mhgtb id    integer :: ntim                    !month time index  integer :: i,j,l                   !indices  integer :: ndata = 1               !input unit  integer :: ndata2 = 2              !input unit  integer :: ncid                    !netCDF file id  integer :: beg4d(4),len4d(4)       !netCDF edge  integer :: dim1_id(1)              !netCDF dimension id for 1-d variables  integer :: dim2_id(2)              !netCDF dimension id for 2-d variables  integer :: dim4_id(4)              !netCDF dimension id for 4-d variables  integer :: status                  !status  character(len=80) :: filei, filei2, fileo  !file names  character(len=80) :: name,unit     !netCDF attributes  real(r8) :: lon(nlon)              !longitude dimension array (1d)  real(r8) :: lat(nlat)              !latitude dimension array (1d)   real(r8) :: longxy(nlon,nlat)      !longitude dimension array (2d)          real(r8) :: latixy(nlon,nlat)      !longitude dimension array (2d)  real(r8) :: edge(4)                !N,E,S,W edges of grid  real(r8) :: dx,dy                  !grid increments!-----------------------------------------------------------------! Determine output filename  fileo = '/ptmp/slevis/lsmv2_2/input/0.5x0.5/new/mksrf_lai.nc'!-----------------------------------------------------------------! Determine input grid!-----------------------------------------------------------------! -----------------------------------------------------------------! Data are 0.5 x 0.5 degree, stored in latitude bands,! from south to north. In a given latitude band, data begin! at the dateline and proceed eastward. So first data! point (x(1,1)) is a box centered at 89.75S, 179.75W!!   89.5S  ---------------------!          |         |         |!          |    x    |    x    |!          |  (1,1)  |  (2,1)  |!          |         |         |!   90.0S  ---------------------!        180.W     179.5W    179.0W! -----------------------------------------------------------------  ! Define North, East, South, West edges of grid    edge(1) =   90.  edge(2) =  180.  edge(3) =  -90.  edge(4) = -180.  ! Make latitudes and longitudes at center of grid cell    dx = (edge(2)-edge(4)) / nlon  dy = (edge(1)-edge(3)) / nlat    do j = 1, nlat     do i = 1, nlon        latixy(i,j) = (edge(3)+dy/2.) + (j-1)*dy        longxy(i,j) = (edge(4)+dx/2.) + (i-1)*dx     end do  end do    lat(:) = latixy(1,:)  lon(:) = longxy(:,1)!-----------------------------------------------------------------! Create netcdf output file!-----------------------------------------------------------------! Open netcdf file    call wrap_create (fileo, nf_clobber, ncid)  call wrap_put_att_text (ncid, nf_global, 'data_type', 'webb_soil_texture_data')  ! Define dimensions    call wrap_def_dim (ncid, 'lon' , nlon        , dimlon_id)  call wrap_def_dim (ncid, 'lat' , nlat        , dimlat_id)  call wrap_def_dim (ncid, 'pft' , numpft+1    , dimpft_id)  call wrap_def_dim (ncid, 'time', nf_unlimited, dimtim_id)  ! Define input file independent variables     name = 'lon'  unit = 'degrees east'  dim1_id(1) = dimlon_id  call wrap_def_var (ncid,'LON', nf_float, 1, dim1_id, lon_id)  call wrap_put_att_text (ncid, lon_id, 'long_name', name)  call wrap_put_att_text (ncid, lon_id, 'units'    , unit)    name = 'lat'  unit = 'degrees north'  dim1_id(1) = dimlat_id  call wrap_def_var (ncid,'LAT', nf_float, 1, dim1_id, lat_id)  call wrap_put_att_text (ncid, lat_id, 'long_name', name)  call wrap_put_att_text (ncid, lat_id, 'units'    , unit)    name = 'longitude-2d'  unit = 'degrees east'  dim2_id(1) = dimlon_id  dim2_id(2) = dimlat_id  call wrap_def_var (ncid, 'LONGXY', nf_float, 2, dim2_id, longxy_id)  call wrap_put_att_text (ncid, longxy_id, 'long_name', name)  call wrap_put_att_text (ncid, longxy_id, 'units'    , unit)  name = 'latitude-2d'  unit = 'degrees north'  dim2_id(1) = dimlon_id  dim2_id(2) = dimlat_id  call wrap_def_var (ncid, 'LATIXY', nf_float, 2, dim2_id, latixy_id)  call wrap_put_att_text (ncid, latixy_id, 'long_name', name)  call wrap_put_att_text (ncid, latixy_id, 'units'    , unit)  name = 'northern edge of surface grid'  unit = 'degrees north'  call wrap_def_var (ncid, 'EDGEN', nf_float, 0, 0, edgen_id)  call wrap_put_att_text (ncid, edgen_id, 'long_name', name)  call wrap_put_att_text (ncid, edgen_id, 'units'    , unit)    name = 'eastern edge of surface grid'  unit = 'degrees east'  call wrap_def_var (ncid, 'EDGEE', nf_float, 0, 0, edgee_id)  call wrap_put_att_text (ncid, edgee_id, 'long_name', name)  call wrap_put_att_text (ncid, edgee_id, 'units'    , unit)    name = 'southern edge of surface grid'  unit = 'degrees north'  call wrap_def_var (ncid, 'EDGES', nf_float, 0, 0, edges_id)  call wrap_put_att_text (ncid, edges_id, 'long_name', name)  call wrap_put_att_text (ncid, edges_id, 'units'    , unit)    name = 'western edge of surface grid'  unit = 'degrees east'  call wrap_def_var (ncid, 'EDGEW', nf_float, 0, 0, edgew_id)  call wrap_put_att_text (ncid, edgew_id, 'long_name', name)  call wrap_put_att_text (ncid, edgew_id, 'units'    , unit)     ! Define input file specific variables    name = 'land mask'  unit = 'unitless'  dim2_id(1) = dimlon_id  dim2_id(2) = dimlat_id  call wrap_def_var (ncid, 'LANDMASK', nf_float, 2, dim2_id, landmask_id)  call wrap_put_att_text (ncid, landmask_id, 'long_name', name)  call wrap_put_att_text (ncid, landmask_id, 'units'    , unit)  dim4_id(1) = dimlon_id  dim4_id(2) = dimlat_id  dim4_id(3) = dimpft_id  dim4_id(4) = dimtim_id  name = 'monthly leaf area index'  unit = 'unitless'  call wrap_def_var (ncid ,'MONTHLY_LAI', nf_float, 4, dim4_id, mlai_id)  call wrap_put_att_text (ncid, mlai_id, 'long_name', name)  call wrap_put_att_text (ncid, mlai_id, 'units'    , unit)       name = 'monthly stem area index'  unit = 'unitless'  call wrap_def_var (ncid ,'MONTHLY_SAI', nf_float, 4, dim4_id, msai_id)  call wrap_put_att_text (ncid, msai_id, 'long_name', name)  call wrap_put_att_text (ncid, msai_id, 'units'    , unit)       name = 'monthly height top'  unit = 'meters'  call wrap_def_var (ncid ,'MONTHLY_HEIGHT_TOP', nf_float, 4, dim4_id, mhgtt_id)  call wrap_put_att_text (ncid, mhgtt_id, 'long_name', name)  call wrap_put_att_text (ncid, mhgtt_id, 'units'    , unit)    name = 'monthly height bottom'  unit = 'meters'  call wrap_def_var (ncid ,'MONTHLY_HEIGHT_BOT', nf_float, 4, dim4_id, mhgtb_id)  call wrap_put_att_text (ncid, mhgtb_id, 'long_name', name)  call wrap_put_att_text (ncid, mhgtb_id, 'units'    , unit)       status = nf_enddef(ncid)! Write grid variables  call wrap_put_var_realx (ncid, lon_id     , lon)  call wrap_put_var_realx (ncid, lat_id     , lat)  call wrap_put_var_realx (ncid, longxy_id  , longxy)  call wrap_put_var_realx (ncid, latixy_id  , latixy)