Hard-coded levels in CALIPSO phase diagnostics and lmd_stats.

[GoogleCodeExporter]

E-mail from Gregory:

Christine Nam, who is working with Johannes Quaas in Hamburg, pointed out an 
issue in the lidar simulator routine “lmd_stats_ipsl.F90”

There is a grid box dimension, called Ntemps, that is hard coded in the 
lmd_stats_ipsl.F90 routine. The value of this dimension, called is 40.
I used 40 because it was the same value as Nlevels and it allowed me to record 
the netcdf files containing the variable using Ntemps, without adding a new 
grid box dimension in COSP.
However, when the USE_VGRID option is set to false and you use your own 
vertical grid, which is different from 40, the code may crash.

Also there is another issue of the same kind. I used hard coded levels in the 
phase part of the routine:
667     do nlev=Nlevels,18,-1  ! from 19.2km until 8.16km
780      do nlev=17,1,-1  ! from 8.16km until 0km

Again with the original grid it will cause issues.

I don’t work at the LMD anymore so I can’t submit through the google code 
project but I’ve enclosed the new routine to this email.
The change is minor. Instead of using the level 18, I put an if condition on 
the pressure level (350hPa) corresponding to 8.16km (=level 18).

Unfortunately, I don’t know how to add the Ntemps dimension to the grid box 
(gbx) value in COSP. Yet, that would be necessary to completely solve the 
problem. “lidarcldtmp" variable should  be defined with Ntemps instead of 
Nlevels. Would you have time to look at it?

Thanks,
Greg

Original issue reported on code.google.com by alejandr...@metoffice.gov.uk on 17 Oct 2014 at 12:19

Attachments:

• lmd_ipsl_stats.F90