--- trunk/libf/dyn3d/start_init_phys_m.f90	2011/08/24 11:43:14	49
+++ trunk/dyn3d/start_init_phys_m.f	2014/07/02 18:39:15	99
@@ -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)
@@ -57,8 +56,7 @@
     ALLOCATE(lon_rad(iml_phys))
     ALLOCATE(lat_rad(jml_phys))
 
-    ! We get the two standard variables
-    ! 'ST': surface temperature
+    ! Surface temperature:
     call nf95_inq_varid(ncid, 'ST', varid)
     call nf95_get_var(ncid, varid, var_ana)
     CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
@@ -66,7 +64,7 @@
          rlatv, tmp_var) 
     tsol_2d = gr_int_dyn(tmp_var)
 
-    ! Soil moisture
+    ! Soil moisture:
     call nf95_inq_varid(ncid, 'CDSW', varid)
     call nf95_get_var(ncid, varid, var_ana)
     CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)