matthias_demuzere@ees_kuleuven_be
New Member
Dear, I am currently deriving the surfdata file for CLM4.5 for a limited area domain on a 0.025x0.025° resolution.I was able to go through all steps, and am currently performing the mkmapdata < namelist step.All goes well until I reach the entry mksrf_fch4 = '/scratch/leuven/300/vsc30061/cesm/inputdata/lnd/clm2/rawdata/mksrf_ch4inversion_360x720_cruncep_simyr2000.c130322.nc'At this step, the routine fails with the following error: Open CH4 parameter file: /scratch/leuven/300/vsc30061/cesm/inputdata/lnd/clm2/rawdata/mksrf_ch4inversion _360x720_cruncep_simyr2000.c130322.nc max_bad_r8 ERROR: f0 = 1.00000000000000 greater than 1.00000000000000 at 13540 max_bad_r8 ERROR: f0 = 1.00000000000000 greater than 1.00000000000000 at 27716 So it seems that during a checking routine, some pixels have values outside a predefined range?What would be your suggestion to go forward with this? In attach you can also find the .log file. Unfortunately, the error itself is not caught in this file, only the already processed entries are listed. [Edit]As a test, I changed max_valid_f0 from 1._r8 to 1.0000005_r8 in routine mkCH4inversionMod.F90.Doing so makes that the mkmapdata routine continues. But now I get another error: Special land units near 100%, but not quite for n,suma = 10620 100.000000000000 Adjusting special land units to 100% pctpft < 0.0 = -2.220446049250313E-016 suma, pcturb_excess, sumpft = 0.000000000000000E+000 72.1264010429353 72.1264010429353 pctpft < 0.0 = -1.421085471520200E-014 suma, pcturb_excess, sumpft = 0.000000000000000E+000 72.1264010429353 72.1264010429353 abort: This error seems to originate from mksurfdat.F90. But I am a bit confused about the values?For clarity I added the log file of this edit, indicated by the _1 at the end of the filename. Thanks.Regards,Matthias