CDO

Hi all,

I am trying to remap data from MOZART hybrid levels to ECMWF with remapeta. After going through the documentation step by step I still seem to miss something because the program aborts with error:

cdo remapeta (Abort): Surface pressure not found!

I have surface pressure in my file but maybe it doesn't have the required variable name (I tried with PS, ps and pres so far).

Can anybody give me some guidance on this?

This is how I ran the program:

cdo remapeta,myvct.txt in.nc out.nc

with myvct.txt from the cdo documentation:

    0       0.00000000000000000       0.00000000000000000 
    1    2000.00000000000000000       0.00000000000000000 
    2    4000.00000000000000000       0.00000000000000000 
    3    6046.10937500000000000       0.00033899326808751 
    4    8267.92968750000000000       0.00335718691349030 
    5   10609.51171875000000000       0.01307003945112228 
    6   12851.10156250000000000       0.03407714888453484 
    7   14698.50000000000000000       0.07064980268478394 
    8   15861.12890625000000000       0.12591671943664551 
    9   16116.23828125000000000       0.20119541883468628 
   10   15356.92187500000000000       0.29551959037780762 
   11   13621.46093750000000000       0.40540921688079834 
   12   11101.55859375000000000       0.52493220567703247 
   13    8127.14453125000000000       0.64610791206359863 
   14    5125.14062500000000000       0.75969839096069336 
   15    2549.96899414062500000       0.85643762350082397 
   16     783.19506835937500000       0.92874687910079956 
   17       0.00000000000000000       0.97298520803451538 
   18       0.00000000000000000       0.99228149652481079 
   19       0.00000000000000000       1.00000000000000000

ncdump -h in.nc looks like this:

dimensions:
        lon = 38 ;
        lat = 26 ;
        mlev = 64 ;
        ilev = 65 ;
variables:
        float lon(lon) ;
                lon:long_name = "longitude" ;
                lon:units = "degrees_east" ;
                lon:standard_name = "longitude" ;
                lon:axis = "X" ;
        float lat(lat) ;
                lat:long_name = "latitude" ;
                lat:units = "degrees_north" ;
                lat:standard_name = "latitude" ;
                lat:axis = "Y" ;
        int mlev(mlev) ;
                mlev:long_name = "hybrid level at layer midpoints" ;
                mlev:standard_name = "hybrid_sigma_pressure" ;
                mlev:units = "level" ;
                mlev:positive = "down" ;
                mlev:formula = "hyam hybm (mlev=hyam+hybm*aps)" ;
                mlev:formula_terms = "ap: hyam b: hybm ps: aps" ;
                mlev:borders = "ilev" ;
        int ilev(ilev) ;
                ilev:long_name = "hybrid level at layer interfaces" ;
                ilev:standard_name = "hybrid_sigma_pressure" ;
                ilev:units = "level" ;
                ilev:positive = "down" ;
                ilev:formula = "hyai hybi (ilev=hyai+hybi*aps)" ;
                ilev:formula_terms = "ap: hyai b: hybi ps: aps" ;
        double hyai(ilev) ;
                hyai:long_name = "hybrid A coefficient at layer interfaces" ;
                hyai:units = "Pa" ;
        double hybi(ilev) ;
                hybi:long_name = "hybrid B coefficient at layer interfaces" ;
                hybi:units = "1" ;
        double hyam(mlev) ;
                hyam:long_name = "hybrid A coefficient at layer midpoints" ;
                hyam:units = "Pa" ;
        double hybm(mlev) ;
                hybm:long_name = "hybrid B coefficient at layer midpoints" ;
                hybm:units = "1" ;
        float DMS_VMR_inst(mlev, lat, lon) ;
                DMS_VMR_inst:units = "VMR" ;
        float NO3_VMR_inst(mlev, lat, lon) ;
                NO3_VMR_inst:units = "VMR" ;
        float PS(lat, lon) ;
                PS:units = "PA" ;

Help is greatly appreciated!

Thanks,
Andreas