33 |
use netcdf95, only: nf95_get_att, nf95_get_var, nf95_inq_varid, & |
use netcdf95, only: nf95_get_att, nf95_get_var, nf95_inq_varid, & |
34 |
nf95_inquire_variable, NF95_OPEN |
nf95_inquire_variable, NF95_OPEN |
35 |
|
|
36 |
REAL, intent(out):: pctsrf(klon, nbsrf) |
REAL, intent(out):: pctsrf(:, :) ! (klon, nbsrf) |
37 |
REAL, intent(out):: ftsol(klon, nbsrf) |
REAL, intent(out):: ftsol(klon, nbsrf) |
38 |
REAL, intent(out):: ftsoil(klon, nsoilmx, nbsrf) |
REAL, intent(out):: ftsoil(klon, nsoilmx, nbsrf) |
39 |
REAL, intent(out):: qsurf(klon, nbsrf) |
REAL, intent(out):: qsurf(klon, nbsrf) |
70 |
integer, intent(out):: ncid_startphy |
integer, intent(out):: ncid_startphy |
71 |
|
|
72 |
! Local: |
! Local: |
|
REAL fractint(klon) |
|
73 |
INTEGER varid, ndims |
INTEGER varid, ndims |
74 |
INTEGER ierr, i |
INTEGER ierr, i |
75 |
|
|
110 |
|
|
111 |
! Verification de l'adequation entre le masque et les sous-surfaces |
! Verification de l'adequation entre le masque et les sous-surfaces |
112 |
|
|
|
fractint = pctsrf(:, is_ter) + pctsrf(:, is_lic) |
|
113 |
DO i = 1 , klon |
DO i = 1 , klon |
114 |
IF ( abs(fractint(i) - masque(i) ) > EPSFRA ) THEN |
IF (abs(pctsrf(i, is_ter) + pctsrf(i, is_lic) - masque(i)) > EPSFRA) THEN |
115 |
print *, 'phyetat0: attention fraction terre pas ', & |
print *, & |
116 |
'coherente ', i, masque(i), pctsrf(i, is_ter), pctsrf(i, is_lic) |
'phyetat0: pctsrf does not agree with masque for continents', & |
117 |
|
i, masque(i), pctsrf(i, is_ter), pctsrf(i, is_lic) |
118 |
ENDIF |
ENDIF |
119 |
END DO |
END DO |
120 |
fractint = pctsrf(:, is_oce) + pctsrf(:, is_sic) |
|
121 |
DO i = 1 , klon |
DO i = 1 , klon |
122 |
IF ( abs( fractint(i) - (1. - masque(i))) > EPSFRA ) THEN |
IF (abs(pctsrf(i, is_oce) + pctsrf(i, is_sic) - (1. - masque(i))) & |
123 |
print *, 'phyetat0 attention fraction ocean pas ', & |
> EPSFRA) THEN |
124 |
'coherente ', i, masque(i) , pctsrf(i, is_oce), pctsrf(i, is_sic) |
print *, 'phyetat0: pctsrf does not agree with masque for ocean ', & |
125 |
|
'and sea-ice', i, masque(i) , pctsrf(i, is_oce), & |
126 |
|
pctsrf(i, is_sic) |
127 |
ENDIF |
ENDIF |
128 |
END DO |
END DO |
129 |
|
|