/[lmdze]/trunk/dyn3d/start_init_phys_m.f
ViewVC logotype

Diff of /trunk/dyn3d/start_init_phys_m.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/start_init_phys_m.f90 revision 49 by guez, Wed Aug 24 11:43:14 2011 UTC trunk/dyn3d/start_init_phys_m.f revision 99 by guez, Wed Jul 2 18:39:15 2014 UTC
# Line 14  CONTAINS Line 14  CONTAINS
14      use dimens_m, only: iim, jjm      use dimens_m, only: iim, jjm
15      use gr_int_dyn_m, only: gr_int_dyn      use gr_int_dyn_m, only: gr_int_dyn
16      use inter_barxy_m, only: inter_barxy      use inter_barxy_m, only: inter_barxy
     use jumble, only: find_longitude, find_latitude  
17      use netcdf, only: nf90_nowrite      use netcdf, only: nf90_nowrite
18      use netcdf95, only: nf95_open, nf95_close, nf95_get_var, nf95_inq_varid, &      use netcdf95, only: nf95_open, nf95_close, nf95_get_var, nf95_inq_varid, &
19           nf95_gw_var           nf95_gw_var, find_coord
20      use nr_util, only: assert, pi      use nr_util, only: assert, pi
21    
22      REAL, intent(out):: tsol_2d(:, :), qsol_2d(:, :) ! (iim + 1, jjm + 1)      REAL, intent(out):: tsol_2d(:, :), qsol_2d(:, :) ! (iim + 1, jjm + 1)
# Line 41  CONTAINS Line 40  CONTAINS
40    
41      call nf95_open('ECPHY.nc', nf90_nowrite, ncid)      call nf95_open('ECPHY.nc', nf90_nowrite, ncid)
42    
43      call find_longitude(ncid, varid=varid)      call find_coord(ncid, varid=varid, std_name="longitude")
44      call nf95_gw_var(ncid, varid, lon_ini)      call nf95_gw_var(ncid, varid, lon_ini)
45      lon_ini = lon_ini * pi / 180. ! convert to rad      lon_ini = lon_ini * pi / 180. ! convert to rad
46      iml_phys = size(lon_ini)      iml_phys = size(lon_ini)
47    
48      call find_latitude(ncid, varid=varid)      call find_coord(ncid, varid=varid, std_name="latitude")
49      call nf95_gw_var(ncid, varid, lat_ini)      call nf95_gw_var(ncid, varid, lat_ini)
50      lat_ini = lat_ini * pi / 180. ! convert to rad      lat_ini = lat_ini * pi / 180. ! convert to rad
51      jml_phys = size(lat_ini)      jml_phys = size(lat_ini)
# Line 57  CONTAINS Line 56  CONTAINS
56      ALLOCATE(lon_rad(iml_phys))      ALLOCATE(lon_rad(iml_phys))
57      ALLOCATE(lat_rad(jml_phys))      ALLOCATE(lat_rad(jml_phys))
58    
59      ! We get the two standard variables      ! Surface temperature:
     ! 'ST': surface temperature  
60      call nf95_inq_varid(ncid, 'ST', varid)      call nf95_inq_varid(ncid, 'ST', varid)
61      call nf95_get_var(ncid, varid, var_ana)      call nf95_get_var(ncid, varid, var_ana)
62      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
# Line 66  CONTAINS Line 64  CONTAINS
64           rlatv, tmp_var)           rlatv, tmp_var)
65      tsol_2d = gr_int_dyn(tmp_var)      tsol_2d = gr_int_dyn(tmp_var)
66    
67      ! Soil moisture      ! Soil moisture:
68      call nf95_inq_varid(ncid, 'CDSW', varid)      call nf95_inq_varid(ncid, 'CDSW', varid)
69      call nf95_get_var(ncid, varid, var_ana)      call nf95_get_var(ncid, varid, var_ana)
70      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)
Legend:
Removed from v.49  
changed lines
  Added in v.99
  ViewVC Help
Powered by ViewVC 1.1.21