Dear Sir:
I would like to run the CLM5-BGC-crop at a high resolution (0.1x0.1) over the user-defined region (eg., USA). I am able to successfully create the relevant global surface dataset at 0.1x0.1 degree following the thread I posted before (high resolution (0.1 degree) CLM5-BGC-crop regional simulation). Then I used
"getregional_datasets.pl" to subset the data over my region of interest.
My create newcase command is as follows:
create_newcase --case BgcCrop2000_BANGhigh --compset 2000_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV \
--res CLM_USRDAT --user-mods-dir $MYDATA_DIR --machine svante --compiler intel \
--run-unsupported
I am able to compile and submit the case. But I got the following running error message:
"surfdata/fatmgrid lon/lat mismatch error",
I went to the log file (see the attachment) and looked at the relevant source code related to the error message:
eal(r8) :: rmaxlon,rmaxlat ! local min/max vars
rmaxlon = 0.0_r8
rmaxlat = 0.0_r8
do n = begg,endg
if (ldomain%lonc(n)-surfdata_domain%lonc(n) > 300.) then
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)-360._r8))
elseif (ldomain%lonc(n)-surfdata_domain%lonc(n) < -300.) then
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)+360._r8))
else
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)))
endif
rmaxlat = max(rmaxlat,abs(ldomain%latc(n)-surfdata_domain%latc(n)))
enddo
if (rmaxlon > 0.001_r8 .or. rmaxlat > 0.001_r8) then
write(iulog,*)' ERROR: surfdata/fatmgrid lon/lat mismatch error', rmaxlon,rmaxlat
call endrun(msg=errMsg(sourcefile, __LINE__))
end if
I don't think I should relax the criterion "if (rmaxlon > 0.001_r8 .or. rmaxlat > 0.001_r8)" to prevent the error. If not, how to solve the problem?
Thanks for your help.
Xiang
I would like to run the CLM5-BGC-crop at a high resolution (0.1x0.1) over the user-defined region (eg., USA). I am able to successfully create the relevant global surface dataset at 0.1x0.1 degree following the thread I posted before (high resolution (0.1 degree) CLM5-BGC-crop regional simulation). Then I used
"getregional_datasets.pl" to subset the data over my region of interest.
My create newcase command is as follows:
create_newcase --case BgcCrop2000_BANGhigh --compset 2000_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV \
--res CLM_USRDAT --user-mods-dir $MYDATA_DIR --machine svante --compiler intel \
--run-unsupported
I am able to compile and submit the case. But I got the following running error message:
"surfdata/fatmgrid lon/lat mismatch error",
I went to the log file (see the attachment) and looked at the relevant source code related to the error message:
eal(r8) :: rmaxlon,rmaxlat ! local min/max vars
rmaxlon = 0.0_r8
rmaxlat = 0.0_r8
do n = begg,endg
if (ldomain%lonc(n)-surfdata_domain%lonc(n) > 300.) then
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)-360._r8))
elseif (ldomain%lonc(n)-surfdata_domain%lonc(n) < -300.) then
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)+360._r8))
else
rmaxlon = max(rmaxlon,abs(ldomain%lonc(n)-surfdata_domain%lonc(n)))
endif
rmaxlat = max(rmaxlat,abs(ldomain%latc(n)-surfdata_domain%latc(n)))
enddo
if (rmaxlon > 0.001_r8 .or. rmaxlat > 0.001_r8) then
write(iulog,*)' ERROR: surfdata/fatmgrid lon/lat mismatch error', rmaxlon,rmaxlat
call endrun(msg=errMsg(sourcefile, __LINE__))
end if
I don't think I should relax the criterion "if (rmaxlon > 0.001_r8 .or. rmaxlat > 0.001_r8)" to prevent the error. If not, how to solve the problem?
Thanks for your help.
Xiang