How to change the onean-land distribution in CESM

[xianglin]

Hi all,
Recently I've been trying to use CESM 1.2.2 to modify some of the land and sea distributions. I am having some problems running gen_domain. Could you please give some suggestions?
./gen_domain -m /glade/u/home/xianglinren/cesmtry/map_gx1PT_TO_fv19_25_aave.231031.nc -o gx1PT -l fv19_25

fmap = /glade/u/home/xianglinren/cesmtry/map_gx1PT_TO_fv19_25_aave.231031.nc


fn1_out_ocn= domain.ocn.gx1PT.231031.nc


fn2_out_lnd= domain.lnd.fv19_25_gx1PT.231031.nc


fn2_out_ocn= domain.ocn.fv19_25_gx1PT.231031.nc


usercomment= null


eps = 9.999999960041972E-013


fminval= 1.000000047497451E-003


fmaxval= 1.00000000000000


set_fv_pole_yc = 0





input SCRIP data...


pole_fix = F





input file =


/glade/u/home/xianglinren/cesmtry/map_gx1PT_TO_fv19_25_aave.231031.nc


open /glade/u/home/xianglinren/cesmtry/map_gx1PT_TO_fv19_25_aave.231031.nc


domain_a = /glade/u/home/xianglinren/cesmtry/gx1PT_2_231031.nc


domain_b = /glade/u/home/xianglinren/cesmtry/fv1.9x2.5_090205.nc


grid_file_ocn= /glade/u/home/xianglinren/cesmtry/gx1PT_2_231031.nc


grid_file_atm= /glade/u/home/xianglinren/cesmtry/fv1.9x2.5_090205.nc


n,nv,ni,nj,na,ns= 122880 4 122880 1 122880


159921





output domain data...


nf = 1


create domain.ocn.gx1PT.231031.nc


netcdf error with rcode = 13 error =


Permission denied





forrtl: severe (174): SIGSEGV, segmentation fault occurred


Image PC Routine Line Source


gen_domain 00000000007CF6D1 Unknown Unknown Unknown


gen_domain 00000000007CD80B Unknown Unknown Unknown


gen_domain 00000000007A79F4 Unknown Unknown Unknown


gen_domain 00000000007A7806 Unknown Unknown Unknown


gen_domain 000000000078CA79 Unknown Unknown Unknown


gen_domain 0000000000748476 Unknown Unknown Unknown


libpthread-2.22.s 00007FD376E3BC00 Unknown Unknown Unknown


libc-2.22.so 00007FD376B82DE1 __strncat_chk Unknown Unknown


gen_domain 000000000074770C Unknown Unknown Unknown


gen_domain 00000000004055ED fmain_IP_check_re 480 gen_domain.F90


gen_domain 000000000040C21B fmain_IP_gen_doma 438 gen_domain.F90


gen_domain 0000000000404B56 MAIN__ 110 gen_domain.F90


gen_domain 000000000040419E Unknown Unknown Unknown


libc-2.22.so 00007FD376AA8A35 __libc_start_main Unknown Unknown


gen_domain 00000000004040A9 Unknown Unknown Unknown


Best,
Xianglin

 

[islas]

Hi Xianglin, do you need to be using CESM1? It is easier to do this in the latest CESM2 tags and if you want to use those I can point to some help on that but those instructions won't work for CESM1.

 

[xianglin]

Hi Islas,

Thanks for the reply. It would be better if it could be done with CESM1. I can also use CESM2.

My purpose is to make part of the ocean into land.
I first modified the kmt as well as the region_mask file, and next generated the corresponding ocn_atm mapping file.
On this basis I tried to run for a month with the results shown. It seems that pop didn't run successfully, but there is also a problem here, there is no log generated in the corresponding folder.


Just wondering what I should do for this situation?

Xianglin

 

[islas]

Ok, maybe someone else on this forum can help with a CESM1 problem. I'm not familiar with how to do this with CESM1. For the more recent tags of CESM2 we are working on instructions for how to do this but the way it works is rather different from CESM1. These pages are still under construction. But here The fsurdat and mask mesh modifier tools for land surface modifications | Community Earth System Model is some description and under "How to use these tools" theres a clickable link for an example of how to fill in the indian ocean with land. I don't think these will help you with your CESM1 problem. But if you decide to move to the most recent tags of CESM2, these will be useful.

 

[slevis]

I hope this points you in the correct direction:

In the error message that you posted, I see

create domain.ocn.gx1PT.231031.nc
netcdf error with rcode = 13 error =
Permission denied

which suggests to me that you may need to ask the file (or directory) owner for read (or write) access.

 

[erik]

I need to point out that at this time CESM1.2.2 is only supported for Paleo climate simulations -- which sounds like what you are doing here. See Support Policy | Community Earth System Model for support details.

Since, this is more of a Paleo Climate issue I'm moving this to that forum.

 

[Siva]

For the gen_domain issue, as suggested in the earlier response, it must be a problem with the user permissions or providing the right path for the input/output files.

Supposing you've solved this already, as you mentioned a different issue in the second post. A failed run at the very beginning could mean a problem with the input files as you're trying to change the kmt, region mask, domain, mapping files, etc... You might want to check these files for inconsistencies.