90 #if ! defined in_arpege
95 CHARACTER(*),
INTENT(in) :: &
97 TYPE(t_dom),
DIMENSION(np),
INTENT(in) :: &
99 TYPE(t_sit),
DIMENSION(nt,np),
INTENT(in) :: &
102 LOGICAL,
DIMENSION(np) :: &
105 zlatn0,zlats0,zehn,zehs,zshn,zshs,zvhn,zvhs
106 REAL,
DIMENSION(np) :: &
114 IF ( nprinto>=2 )
THEN
120 WRITE(noutlu,*)
' **** gltools_chkglo_r ****'
132 zfsit(:) = sum( tpsit(:,:)%fsi, dim=1 )
134 zhsiw(:) = sum( tpsit(:,:)%fsi*tpsit(:,:)%hsi, dim=1 )
141 WHERE( tpdom(:)%lat>zlatn0 )
145 WHERE( tpdom(:)%lat<zlats0 )
157 zehn = sum( tpdom(:)%srf, mask=(ghnorth(:).AND.zfsit(:)>xfsic) ) / 1.e+12
158 zehs = sum( tpdom(:)%srf, mask=(ghsouth(:).AND.zfsit(:)>xfsic) ) / 1.e+12
164 zshn = sum( tpdom(:)%srf*zfsit(:), mask=ghnorth(:) ) / 1.e+12
165 zshs = sum( tpdom(:)%srf*zfsit(:), mask=ghsouth(:) ) / 1.e+12
171 zvhn = sum( tpdom(:)%srf*zhsiw(:), mask=ghnorth(:) ) / 1.e+12
172 zvhs = sum( tpdom(:)%srf*zhsiw(:), mask=ghsouth(:) ) / 1.e+12
179 #if ! defined in_arpege
188 WRITE(noutlu,*)
' North South'
189 WRITE(noutlu,1000) zshn,zshs
190 WRITE(noutlu,1100) zehn,zehs
191 WRITE(noutlu,1200) zvhn,zvhs
202 1000
FORMAT( 5x,
"Ice surface (SISH.SIG)",2(4x,f9.5) )
203 1100
FORMAT( 5x,
"Ice extent (SIEH.SIG)",2(4x,f9.5) )
204 1200
FORMAT( 5x,
"Ice volume (SIVH.SIG)",2(4x,f9.5) )