output doc

[gitpeterwind]

Hi @avaldebe ( CC @mifads ) Do we have a documentation for outputting the data? (config_Output.nml), I didn't see any in the user guide. Specifically, is it possible to output a group in 3D? (it didn't work when I tried).

Also I wanted to add some special things in the "outputting syntax" chapter. I have put it in Log.changes 24 May for now.

[avaldebe]

It is undocumented, as far as I remember. You can look at the CAMS50 output namelist on the research repo for an example of 3D output.

[gitpeterwind]

OK, thanks. Do we have a separate documentation for our internal use? For 3D it worked for 'ug', but not for 'ppb'. (CAMS50 had only ug). I think I can do with ug for now.

[avaldebe]

We do not have a separate documentation for internal use. We struggle to maintain one documentation as it is. Please use the docs branch in this repository for documenting the 3D output, as it will be the official documentation for the next open source release.

As for the supported units, the 3D output should work for any unit supported by Group_Units (Units_ml.f90). Can you share (on this tread) the 3D entries on your OutputConcs_config namelist? And the error message?

[gitpeterwind]

this line: 'OX' ,'ppb' ,'3d' ,'AIR_CONCS','GROUP','Y',! give this error: USET XNCOL FOUND 6 60 CH4 4 60 *** NOT FOUND 3 28D3_ppb_OX;Class:3D_PPB_GROUP STOP-ALL ERROR: Derived 3D class NOT FOUND 3 28D3_ppb_OX3D_PPB_GROUP application called MPI_Abort(MPI_COMM_WORLD, 9) - process 15 forrtl: error (69): process interrupted (SIGINT)

When 'ppb' is replaced by 'ug', it works fine

[avaldebe]

I see, the "3D_PPB_GROUP" case is missing in in Derived (Derived_ml.90). I'll see if it can be easily added.

[avaldebe]

I committed the patch to our research repo. It is only a single line change (my text editor removed some extra trailing spaces). Run a single day and could output D3_ug_OX and D3_ppb_OX.

Please note that by default all model levels are outputted following the model level idexing, ie from model top to bottom. This can be modified in OutputSize_config namelist. The following example, from the CAMS50 output configuration, selects 8 levels indexed from the model surface:

&OutputSize_config
  num_lev3d =8,
  lev3d=1,2,3,4,6,9,10,12,
  lev3d_from_surface=T,
&end

[gitpeterwind]

Thanks, worked perfectly now! Yes, I output only the 7 lowest layers, using the "lev3d" settings.