4 |
|
|
5 |
IMPLICIT none |
IMPLICIT none |
6 |
|
|
7 |
REAL, save:: rlat(klon), rlon(klon) |
REAL, save, protected:: rlat(klon), rlon(klon) |
8 |
! latitude and longitude of a point of the scalar grid identified |
! latitude and longitude of a point of the scalar grid identified |
9 |
! by a simple index, in degrees |
! by a simple index, in degrees |
10 |
|
|
11 |
integer, save:: itau_phy |
integer, save, protected:: itau_phy |
12 |
|
REAL, save, protected:: zmasq(KLON) ! fraction of land |
13 |
|
|
14 |
private klon |
private klon |
15 |
|
|
26 |
! Objet : lecture de l'état initial pour la physique |
! Objet : lecture de l'état initial pour la physique |
27 |
|
|
28 |
USE conf_gcm_m, ONLY: raz_date |
USE conf_gcm_m, ONLY: raz_date |
29 |
use dimphy, only: zmasq, klev |
use dimphy, only: klev |
30 |
USE dimsoil, ONLY : nsoilmx |
USE dimsoil, ONLY : nsoilmx |
31 |
USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf |
USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf |
32 |
use netcdf, only: nf90_global, nf90_inq_varid, NF90_NOERR, NF90_NOWRITE |
use netcdf, only: nf90_global, nf90_inq_varid, NF90_NOERR, NF90_NOWRITE |
357 |
|
|
358 |
END SUBROUTINE phyetat0 |
END SUBROUTINE phyetat0 |
359 |
|
|
360 |
|
!********************************************************************* |
361 |
|
|
362 |
|
subroutine phyetat0_new |
363 |
|
|
364 |
|
use nr_util, only: pi |
365 |
|
|
366 |
|
use dimensions, only: iim, jjm |
367 |
|
use dynetat0_m, only: rlatu, rlonv |
368 |
|
use grid_change, only: dyn_phy |
369 |
|
USE start_init_orog_m, only: mask |
370 |
|
|
371 |
|
!------------------------------------------------------------------------- |
372 |
|
|
373 |
|
rlat(1) = 90. |
374 |
|
rlat(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.) * 180. / pi |
375 |
|
! (with conversion to degrees) |
376 |
|
rlat(klon) = - 90. |
377 |
|
|
378 |
|
rlon(1) = 0. |
379 |
|
rlon(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.) * 180. / pi |
380 |
|
! (with conversion to degrees) |
381 |
|
rlon(klon) = 0. |
382 |
|
|
383 |
|
zmasq = pack(mask, dyn_phy) |
384 |
|
itau_phy = 0 |
385 |
|
|
386 |
|
end subroutine phyetat0_new |
387 |
|
|
388 |
end module phyetat0_m |
end module phyetat0_m |