convert_mksrf.F90

 

filenameresolution
Input: topography/ice 10min
Output:    
  mksrf_lanwat.myrun.nc 0.5o
  mksrf_pft.myrun.nc 0.5o
  mksrf_glacier.myrun.nc 0.5o

 

Usage: 

  1.  cp convert_mksrf.template convert_mksrf.template.myrun 
  2.  modify convert_mksrf.template.myrun to point to your input files 
  3. cp convert_mksrf.template.myrun convert_mksrf.F90
  4.  compile (gmake)
  5.  execute:  ./convert_mksrf

 

Input File format for convert_mksrf.F90

 

        float TOP(lat, lon) ;

                TOP:units = "meter" ;

                TOP:long_name = "10-min elevation from USGS 30-sec dataset" ;

        float lat(lat) ;

                lat:long_name = "lat" ;

                lat:units = "degrees north" ;

        float lon(lon) ;

                lon:units = "degrees east" ;

                lon:long_name = "lon" ;

        float ICE(lat, lon) ;

                ICE:units = "%" ;

                ICE:long_name = "Fraction of Grid Cell Covered with Glacier" ;

                ICE:valid_min = 0.f ;

                ICE:valid_max = 100.f ;

        float LANDFRAC(lat, lon) ;

                LANDFRAC:units = "meter" ;

                LANDFRAC:long_name = "10-min elevation from USGS 30-sec dataset" ;

                LANDFRAC:valid_max = 1 ;

                LANDFRAC:valid_min = 0 ;

        float LANDMASK(lat, lon) ;

                LANDMASK:units = "meter" ;

                LANDMASK:long_name = "10-min elevation from USGS 30-sec dataset" ;

                LANDMASK:valid_max = 1 ;

                LANDMASK:valid_min = 0 ;