Changeset 2718 for trunk/NEMOGCM/TOOLS/WEIGHTS/nocsutil/scripgrid_mod.F90

Timestamp:

2011-04-04T12:03:24+02:00 (13 years ago)

Author:

sga

Message:

NEMO trunk: small corrections to the weights generation code in the TOOLS directory for nemo 3.3

• typos in README file corrected
• scripgrid.F90 now allows any input nemo grid filename (but lat/lon names have to begin with glam/gphi)
• scripinterp.F90 now actually does something

File:

• trunk/NEMOGCM/TOOLS/WEIGHTS/nocsutil/scripgrid_mod.F90 (modified) (1 diff)

trunk/NEMOGCM/TOOLS/WEIGHTS/nocsutil/scripgrid_mod.F90

@@ -92 +92 @@
     call release_unit(iunit)
 
-    if (trim(nemo_file) .eq. 'coordinates.nc') then
-
-      ! set up the names of the corner variables for a given input
-      ! the offset represents what needs to be added to (i,j) to get to the correct 
-      ! element in the corner arrays to correspond to the point northeast of the center
-      if (nemo_lon(5:5) == "t") then
-        corn_lon = "glamf"
-        corn_lat = "gphif"
-        offset = (/ 0,0 /)
-      else if (nemo_lon(5:5) == "u") then
-        corn_lon = "glamv"
-        corn_lat = "gphiv"
-        offset = (/ 1,0 /)
-      else if (nemo_lon(5:5) == "v") then
-        corn_lon = "glamu"
-        corn_lat = "gphiu"
-        offset = (/ 0,1 /)
-      else
-        write(6,*) 'unknown nemo_lon name'
-        stop
-      endif
-
-      write(6,*) "processing " // trim(nemo_file)
-      call convertNEMO(nemo_file, nemo_lon, nemo_lat, corn_lon, corn_lat, &
-                       offset, nemogrid_file)
-
+    if (nemo_lon(1:4) .ne. 'glam' .or. nemo_lat(1:4) .ne. 'gphi') then
+      write(6,*) 'lon name does not start with "glam" or lat name does not start with "gphi"'
+      stop
+    endif
+
+    ! set up the names of the corner variables for a given input
+    ! the offset represents what needs to be added to (i,j) to get to the correct 
+    ! element in the corner arrays to correspond to the point northeast of the center
+    if (nemo_lon(5:5) == "t") then
+      corn_lon = "glamf"
+      corn_lat = "gphif"
+      offset = (/ 0,0 /)
+    else if (nemo_lon(5:5) == "u") then
+      corn_lon = "glamv"
+      corn_lat = "gphiv"
+      offset = (/ 1,0 /)
+    else if (nemo_lon(5:5) == "v") then
+      corn_lon = "glamu"
+      corn_lat = "gphiu"
+      offset = (/ 0,1 /)
     else
-
-      write(6,*) "processing regular grid " // trim(nemo_file)
-      call convertFLUX(nemo_file, nemo_lon, nemo_lat, &
-                       nemo_mask, nemo_mask_value, nemogrid_file)
-
+      write(6,*) 'unknown nemo_lon name'
+      stop
     endif
+
+    write(6,*) "processing " // trim(nemo_file)
+    call convertNEMO(nemo_file, nemo_lon, nemo_lat, corn_lon, corn_lat, &
+                     offset, nemogrid_file)
 
     write(6,*) "processing regular grid"