Land and Ocean Masks

[rasmus@uoregon_edu]

I'm trying to build a regional land model for a high resolution 10x10km grid on our local machines and I believe I have sorted out most of the issues. When I use mknoocnmap.pl to create a SCRIP grid everything works fine.  However, when I create my own SCRIP grid with equal area cells I run into a problem: I've messed up the land/ocean mask somehow and the land model is being run over the Pacific so my masks (perhaps only one) seem to be inverted.  But when I look at the surfdata file the mask seems to be correct (i.e., land where land should be). When I look at other SCRIP grids (e.g., SCRIPgrid_0.23x0.31_nomask_c110308.nc) the grid_imask values are all 1.0.  Perhaps I should set the grid_imask values to 1 in the SCRIPgrid_10x10km_westernUS_nomask_c150812.nc file I created. Thanks in advance,Craig

 

[chakravorty_aniket@gmail_com]

Hello,I am also trying to build a regional land surface model for a high resolution at 25x25km grid.However, I am not familiar with the use of mknoocnmap.pl script for regional simulations. The tutorial and README only talks about the procedure for point simulation. Could you tell me the syntax (with the all the necessary arguments) of mknoocnmap.pl for a regional simulation, as an example.I would really appreciate the help Aniket Chakravorty 

 

[rasmus@uoregon_edu]

I think I've figured out how to solve most of our issues.  You can run mknooconmap as shown below for a regular lat/long grid.  We are running only CLM so I don't think the distinction between points and a regional simulation is meaningful as there is no atmospheric model that needs communication between cells.  In fact this creates a SCRIPgrid file so it can be also be used when running an atmospheric simulation as well (I think).  Here is command I used for a 211x232 "point" grid:./mknoocnmap.pl -v -p 39.1,248.5 -n 211x232_westernUS -nx 211 -ny 232 -dx 30.0 -dy 22.0However, this really is no-ocean as the mask is set to 1 everywhere.  Since a portion of this grid is over the Pacific, I thought that this would create problems but I think mksurfdata.pl works the mask out as I assume the input data to mksurfdata.pl has the correct ocean/land masks.We ended up creating our own grid programatically in the end.  I can send you a copy of the ncl scripts if you'd like.  Although I didn't write them (I would have just used Fortran or C as I don't know ncl very well).We had wanted to use a projected (equal area) grid but had to give up in the end because I couldn't figure out how to get CESM to use a projected grid correctly without hacking code (I was willing but we decided it was probably easier and safer to map back to an equal area grid after running CESM).Craig Rasmussen

 

[tariqmme23@hotmail_com]

Hello,i am trying to generating surface dataset for a single point simulations using CLM4.5 but did not get success.i use the following command to generate SCRIP Files./mknoocnmap.pl -p 25.67,68.72 -n test_tndjm
and all detail about the running of this command is here.
Create SCRIP grid and mapping files for a single-pointCreate land SCRIP gridfileExecute: env S_LAT=25.62 W_LON=68.67 N_LAT=25.72 E_LON=68.77 NX=1 NY=1 PTNAME=test_tndjm   GRIDFILE=/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_nomask_c160906.nc ncl mkscripgrid.ncl Copyright (C) 1995-2015 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.3.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details.warning:stringtointeger: a bad value was passed; input strings must contain numeric digits, replacing with missing valuefatal:fspan: number of elements parameter is less-than-or-equal-to one, can't continue fatal:["Execute.c":8575]:Execute: Error occurred at or near line 99 in file mkscripgrid.ncl
Create ocean SCRIP gridfileExecute: env S_LAT=25.62 W_LON=68.67 N_LAT=25.72 E_LON=68.77 NX=1 NY=1 PTNAME=test_tndjm   GRIDFILE=/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_noocean_c160906.nc IMASK=0 ncl mkscripgrid.ncl Copyright (C) 1995-2015 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.3.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details.fatal:fspan: number of elements parameter is less-than-or-equal-to one, can't continue fatal:["Execute.c":8575]:Execute: Error occurred at or near line 99 in file mkscripgrid.ncl
Create unity mapping file between the two gridfileExecute: env GRIDFILE1=/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_noocean_c160906.nc GRIDFILE2=/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_nomask_c160906.nc MAPFILE=map_test_tndjm_noocean_to_test_tndjm_nomask_aave_da_160906.nc  ncl /public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapdata/mkunitymap.ncl Copyright (C) 1995-2015 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.3.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. fatal:["FileSupport.c":3731]:_NclFindFileExt: Requested file or does not exist
fatal:["FileSupport.c":4159]:(/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_noocean_c160906.nc) has no file extension, can't determine type of file to openfatal:["FileSupport.c":3731]:_NclFindFileExt: Requested file or does not exist fatal:["FileSupport.c":4159]:(/public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapgrids/SCRIPgrid_test_tndjm_nomask_c160906.nc) has no file extension, can't determine type of file to openfatal:file (nca) isn't definedfatal:["Execute.c":8575]:Execute: Error occurred at or near line 64 in file /public/xzhtest/tariq/cesm1_2_0/models/lnd/clm/tools/shared/mkmapdata/mkunitymap.ncl   Successfully created grid/mapping files for single-point
but at the end there is no SCRIP GRIDFILES generated
your help will be highly appreciated
RegardsTariq