Changeset 8624a75 in flexpart.git for flexpart_code/grib2nc4/processmetfields.F90

Timestamp:

Aug 24, 2017, 7:07:08 PM (7 years ago)

Author:

Don Morton <Don.Morton@…>

Branches:

FPv9.3.2, grib2nc4_repair

Children:

857dfd0

Parents:

4c0504c

Message:

Enhancements to FPv9.3.2

Documented in Ticket #182 (as well as CTBTO ticket 357)

File:

• flexpart_code/grib2nc4/processmetfields.F90 (modified) (4 diffs)

flexpart_code/grib2nc4/processmetfields.F90

@@ -1 @@
-subroutine processmetfields(ind,metdata_format)
+subroutine processmetfields(ind,metdata_format,coordinates_provided, coordX, coordY)
+
   !                       i     o
   !*****************************************************************************
@@ -37 +38 @@
   integer :: dumpData
   character(len=512):: fpfname, dumpPath, filename
+  integer :: coordX, coordY
+  logical :: coordinates_provided
 
 
@@ -107 +110 @@
              call calcpar_ecmwf(memind(1),uuh,vvh,pvh) 
              call calcpar_nests(memind(1),uuhn,vvhn,pvhn,metdata_format)
-             call verttransform_ecmwf(memind(1),uuh,vvh,wwh,pvh)
+             if ( coordinates_provided ) then
+                 call verttransform_grib2nc4_ecmwf(memind(1),uuh,vvh,wwh,pvh,coordX,coordY)
+             else
+                 call verttransform_grib2nc4_ecmwf(memind(1),uuh,vvh,wwh,pvh,-1,-1)
+             endif
              call verttransform_nests(memind(1),uuhn,vvhn,wwhn,pvhn)
              memtime(1)=wftime(ind)
@@ -116 +123 @@
              call calcpar_gfs(memind(1),uuh,vvh,pvh)
              call calcpar_nests(memind(1),uuhn,vvhn,pvhn,metdata_format)
-             call verttransform_gfs(memind(1),uuh,vvh,wwh,pvh)
+             if ( coordinates_provided ) then
+                 call verttransform_grib2nc4_gfs(memind(1),uuh,vvh,wwh,pvh,coordX,coordY)
+             else
+                 call verttransform_grib2nc4_gfs(memind(1),uuh,vvh,wwh,pvh,-1,-1)
+             endif
              call verttransform_nests(memind(1),uuhn,vvhn,wwhn,pvhn)
              memtime(1)=wftime(ind)