--- trunk/libf/dyn3d/start_init_phys_m.f90	2011/08/24 11:43:14	49
+++ trunk/dyn3d/start_init_phys_m.f	2014/03/05 14:57:53	82
@@ -14,10 +14,9 @@
     use dimens_m, only: iim, jjm
     use gr_int_dyn_m, only: gr_int_dyn
     use inter_barxy_m, only: inter_barxy
-    use jumble, only: find_longitude, find_latitude
     use netcdf, only: nf90_nowrite
     use netcdf95, only: nf95_open, nf95_close, nf95_get_var, nf95_inq_varid, &
-         nf95_gw_var
+         nf95_gw_var, find_coord
     use nr_util, only: assert, pi
 
     REAL, intent(out):: tsol_2d(:, :), qsol_2d(:, :) ! (iim + 1, jjm + 1)
@@ -41,12 +40,12 @@
 
     call nf95_open('ECPHY.nc', nf90_nowrite, ncid)
 
-    call find_longitude(ncid, varid=varid)
+    call find_coord(ncid, varid=varid, std_name="longitude")
     call nf95_gw_var(ncid, varid, lon_ini)
     lon_ini = lon_ini * pi / 180. ! convert to rad
     iml_phys = size(lon_ini)
 
-    call find_latitude(ncid, varid=varid)
+    call find_coord(ncid, varid=varid, std_name="latitude")
     call nf95_gw_var(ncid, varid, lat_ini)
     lat_ini = lat_ini * pi / 180. ! convert to rad
     jml_phys = size(lat_ini)