81 ( kinit,omsg,tpdom,tptfl,tpsit,pice_a,pemps_a,psalt_a,psalf_a)
90 INTEGER,
INTENT(in) :: &
92 CHARACTER(*),
INTENT(in) :: &
94 TYPE(t_dom),
DIMENSION(nx,ny),
INTENT(in) :: &
96 TYPE(t_tfl),
DIMENSION(nx,ny),
INTENT(in) :: &
98 TYPE(t_sit),
DIMENSION(nt,nx,ny),
INTENT(in) :: &
100 REAL,
INTENT(inout) :: &
101 pemps_a, pice_a, psalt_a, psalf_a
103 REAL,
DIMENSION(nx,ny) :: &
104 zice, zemps, zsalt, zsalf
106 zice_a, zemps_a, zdemps, zdmice, zaddterm_a, zsalt_a, zsalf_a, zdmsalt, zdsalf
110 zice(:,:) = rhoice*sum( tpsit(:,:,:)%fsi*tpsit(:,:,:)%hsi, dim=1 )
111 zice_a =
glt_avg(tpdom, zice(:,:), 1)
112 zemps(:,:) = tptfl(:,:)%cio
113 zemps_a =
glt_avg(tpdom, zemps(:,:), 1)
115 zsalt(:,:) = rhoice*sum( tpsit(:,:,:)%fsi*tpsit(:,:,:)%hsi*tpsit(:,:,:)%ssi, dim=1 )*1.e-3
116 zsalt_a =
glt_avg(tpdom, zsalt(:,:), 1)
118 zsalf(:,:) = tptfl(:,:)%sio
119 zsalf_a =
glt_avg(tpdom, zsalf(:,:), 1)
121 zdmice = ( zice_a - pice_a) / dtt
122 zdemps = zemps_a - pemps_a
123 zdmsalt = ( zsalt_a - psalt_a) / dtt
124 zdsalf = zsalf_a - psalf_a
127 '--------------------------------------------------------------------'
129 WRITE(noutlu,*)
' Change in ice content :', zdmice
130 WRITE(noutlu,*)
' Change in emps :', zdemps
131 WRITE(noutlu,*)
' BILAN DMICE-EMP :', zdmice + zdemps
132 WRITE(noutlu,*) omsg,
' Salt Content :', zsalt_a
133 WRITE(noutlu,*)
' Change in Salt Content :', zdmsalt
134 WRITE(noutlu,*)
' Change in salt flux :', zdsalf
135 WRITE(noutlu,*)
' BILAN Dsalt_content/salflx :', zdmsalt + zdsalf
139 WRITE(noutlu,*) omsg,
' Salt Content :', zsalt_a
140 WRITE(noutlu,*) omsg,
' Mass Content :', zice_a
subroutine glt_updice(kinit, omsg, tpdom, tptfl, tpsit, pice_a, pemps_a, psalt_a, psalf_a)