michael_veres@yahoo_com
New Member
Hello all,
I am running CESM1.0.4 and using CAM4 physics, CLM and DOCN using a highly idealized landmask and land fraction. I have been able to get the model to run, but I have noticed that the reference height values are behaving strangely. Anytime LANDFRAC is between 0 and 1, it appears to be scaling the values. The following reference height variables are a small printout of a few TREFHT values and the corresponding LANDFRAC.
(7372) LANDFRAC = 0; TREFHT = 278.667
(7373) LANDFRAC = 0; TREFHT = 278.588
(7374) LANDFRAC = 0.754609; TREFHT = 209.526
(7375) LANDFRAC = 1; TREFHT = 277.95
(7376) LANDFRAC = 1; TREFHT = 277.843
(7377) LANDFRAC = 0.84133; TREFHT = 233.3
(7378 ) LANDFRAC = 0.0947464; TREFHT = 26.0883
(7379) LANDFRAC = 0; TREFHT = 277.458
(7380) LANDFRAC = 0; TREFHT = 277.711
(7381) LANDFRAC = 0; TREFHT = 278.157
It is (apparently) calculating the all-ocean and all-land points correctly, but the others are relatively proportional to their LANDFRAC. Does anybody know what could cause this? I also get the same issue with QREFHT. All other (non-reference height) variables (including the 3D temperature) appear fine. The only changes to the source code that I made was to modify the dry air mass quantity in inidat.F90. I did not make any changes regarding the reference height calculations.
I don't know if this is a bug or if I'm misunderstanding something. Any insight on this would be appreciated. I am using Bluefire, so I can provide a file path for the data files for NCAR users.
mveres
I am running CESM1.0.4 and using CAM4 physics, CLM and DOCN using a highly idealized landmask and land fraction. I have been able to get the model to run, but I have noticed that the reference height values are behaving strangely. Anytime LANDFRAC is between 0 and 1, it appears to be scaling the values. The following reference height variables are a small printout of a few TREFHT values and the corresponding LANDFRAC.
(7372) LANDFRAC = 0; TREFHT = 278.667
(7373) LANDFRAC = 0; TREFHT = 278.588
(7374) LANDFRAC = 0.754609; TREFHT = 209.526
(7375) LANDFRAC = 1; TREFHT = 277.95
(7376) LANDFRAC = 1; TREFHT = 277.843
(7377) LANDFRAC = 0.84133; TREFHT = 233.3
(7378 ) LANDFRAC = 0.0947464; TREFHT = 26.0883
(7379) LANDFRAC = 0; TREFHT = 277.458
(7380) LANDFRAC = 0; TREFHT = 277.711
(7381) LANDFRAC = 0; TREFHT = 278.157
It is (apparently) calculating the all-ocean and all-land points correctly, but the others are relatively proportional to their LANDFRAC. Does anybody know what could cause this? I also get the same issue with QREFHT. All other (non-reference height) variables (including the 3D temperature) appear fine. The only changes to the source code that I made was to modify the dry air mass quantity in inidat.F90. I did not make any changes regarding the reference height calculations.
I don't know if this is a bug or if I'm misunderstanding something. Any insight on this would be appreciated. I am using Bluefire, so I can provide a file path for the data files for NCAR users.
mveres