| 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:: zmasq(KLON) ! fraction of land |
REAL, save, protected:: zmasq(KLON) ! fraction of land |
| 13 |
|
|
| 14 |
private klon |
private klon |
| 15 |
|
|
| 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 |
Parent Directory
| 
Revision Log
| 
Patch