Changeset 888 for trunk/NEMO/LIM_SRC_2/limwri_dimg_2.h90
- Timestamp:
- 2008-04-11T19:05:03+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/limwri_dimg_2.h90
r823 r888 2 2 !!---------------------------------------------------------------------- 3 3 !! LIM 2.0, UCL-LOCEAN-IPSL (2005) 4 !! $ Header$4 !! $ Id: $ 5 5 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 6 6 !!---------------------------------------------------------------------- … … 82 82 83 83 zsto = rdt_ice 84 zout = nwrite * rdt_ice / n fice84 zout = nwrite * rdt_ice / nn_fsbc 85 85 zsec = 0. 86 86 niter = 0 … … 106 106 zcmo(ji,jj,5) = sist (ji,jj) 107 107 zcmo(ji,jj,6) = fbif (ji,jj) 108 zcmo(ji,jj,7) = zindb * ( u _ice(ji,jj ) * tmu(ji,jj ) + u_ice(ji+1,jj ) * tmu(ji+1,jj ) &109 + u_ice(ji,jj+1) * tmu(ji,jj+1) + u_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &108 zcmo(ji,jj,7) = zindb * ( ui_ice(ji,jj ) * tmu(ji,jj ) + ui_ice(ji+1,jj ) * tmu(ji+1,jj ) & 109 & + ui_ice(ji,jj+1) * tmu(ji,jj+1) + ui_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & 110 110 / ztmu 111 111 112 zcmo(ji,jj,8) = zindb * ( v _ice(ji,jj ) * tmu(ji,jj ) + v_ice(ji+1,jj ) * tmu(ji+1,jj ) &113 + v_ice(ji,jj+1) * tmu(ji,jj+1) + v_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &112 zcmo(ji,jj,8) = zindb * ( vi_ice(ji,jj ) * tmu(ji,jj ) + vi_ice(ji+1,jj ) * tmu(ji+1,jj ) & 113 & + vi_ice(ji,jj+1) * tmu(ji,jj+1) + vi_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & 114 114 / ztmu 115 zcmo(ji,jj,9) = sst_ io(ji,jj)116 zcmo(ji,jj,10) = sss_ io(ji,jj)117 118 zcmo(ji,jj,11) = fnsolar(ji,jj) + fsolar(ji,jj)119 zcmo(ji,jj,12) = fsolar(ji,jj)120 zcmo(ji,jj,13) = fnsolar(ji,jj)115 zcmo(ji,jj,9) = sst_m(ji,jj) 116 zcmo(ji,jj,10) = sss_m(ji,jj) 117 118 zcmo(ji,jj,11) = qns(ji,jj) + qsr(ji,jj) 119 zcmo(ji,jj,12) = qsr(ji,jj) 120 zcmo(ji,jj,13) = qns(ji,jj) 121 121 ! See thersf for the coefficient 122 zcmo(ji,jj,14) = - fsalt(ji,jj) * rday * ( sss_io(ji,jj) + epsi16 ) / soce123 zcmo(ji,jj,15) = gtaux(ji,jj)124 zcmo(ji,jj,16) = gtauy(ji,jj)125 zcmo(ji,jj,17) = ( 1.0 - frld(ji,jj) ) * qsr_ice (ji,jj) + frld(ji,jj) * qsr_oce(ji,jj)126 zcmo(ji,jj,18) = ( 1.0 - frld(ji,jj) ) * qnsr_ice(ji,jj) + frld(ji,jj) * qnsr_oce(ji,jj)122 zcmo(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 123 zcmo(ji,jj,15) = utaui_ice(ji,jj) 124 zcmo(ji,jj,16) = vtaui_ice(ji,jj) 125 zcmo(ji,jj,17) = qsr_ice(ji,jj) 126 zcmo(ji,jj,18) = qns_ice(ji,jj) 127 127 zcmo(ji,jj,19) = sprecip(ji,jj) 128 128 END DO … … 132 132 nmoyice = nmoyice + 1 133 133 ! compute mean value if it is time to write on file 134 IF ( MOD(kt+n fice-1-nit000+1,nwrite) == 0 ) THEN134 IF ( MOD(kt+nn_fsbc-1-nit000+1,nwrite) == 0 ) THEN 135 135 rcmoy(:,:,:) = rcmoy(:,:,:) / FLOAT(nmoyice) 136 136 #else 137 IF ( MOD(kt-n fice-1-nit000+1,nwrite) == 0 ) THEN137 IF ( MOD(kt-nn_fsbc-1-nit000+1,nwrite) == 0 ) THEN 138 138 ! case of instantaneaous output rcmoy(:,:, 1:jpnoumax ) = 0.e0 139 139 DO jj = 2 , jpjm1 … … 149 149 rcmoy(ji,jj,5) = sist (ji,jj) 150 150 rcmoy(ji,jj,6) = fbif (ji,jj) 151 rcmoy(ji,jj,7) = zindb * ( u _ice(ji,jj ) * tmu(ji,jj ) + u_ice(ji+1,jj ) * tmu(ji+1,jj ) &152 + u_ice(ji,jj+1) * tmu(ji,jj+1) + u_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &151 rcmoy(ji,jj,7) = zindb * ( ui_ice(ji,jj ) * tmu(ji,jj ) + ui_ice(ji+1,jj ) * tmu(ji+1,jj ) & 152 & + ui_ice(ji,jj+1) * tmu(ji,jj+1) + ui_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & 153 153 / ztmu 154 154 155 rcmoy(ji,jj,8) = zindb * ( v _ice(ji,jj ) * tmu(ji,jj ) + v_ice(ji+1,jj ) * tmu(ji+1,jj ) &156 + v_ice(ji,jj+1) * tmu(ji,jj+1) + v_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &155 rcmoy(ji,jj,8) = zindb * ( vi_ice(ji,jj ) * tmu(ji,jj ) + vi_ice(ji+1,jj ) * tmu(ji+1,jj ) & 156 & + vi_ice(ji,jj+1) * tmu(ji,jj+1) + vi_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & 157 157 / ztmu 158 rcmoy(ji,jj,9) = sst_ io(ji,jj)159 rcmoy(ji,jj,10) = sss_ io(ji,jj)160 161 rcmoy(ji,jj,11) = fnsolar(ji,jj) + fsolar(ji,jj)162 rcmoy(ji,jj,12) = fsolar(ji,jj)163 rcmoy(ji,jj,13) = fnsolar(ji,jj)158 rcmoy(ji,jj,9) = sst_m(ji,jj) 159 rcmoy(ji,jj,10) = sss_m(ji,jj) 160 161 rcmoy(ji,jj,11) = qns(ji,jj) + qsr(ji,jj) 162 rcmoy(ji,jj,12) = qsr(ji,jj) 163 rcmoy(ji,jj,13) = qns(ji,jj) 164 164 ! See thersf for the coefficient 165 rcmoy(ji,jj,14) = - fsalt(ji,jj) * rday * ( sss_io(ji,jj) + epsi16 ) / soce166 rcmoy(ji,jj,15) = gtaux(ji,jj)167 rcmoy(ji,jj,16) = gtauy(ji,jj)168 rcmoy(ji,jj,17) = ( 1.0 - frld(ji,jj) ) * qsr_ice (ji,jj) + frld(ji,jj) * qsr_oce(ji,jj)169 rcmoy(ji,jj,18) = ( 1.0 - frld(ji,jj) ) * qnsr_ice(ji,jj) + frld(ji,jj) * qnsr_oce(ji,jj)165 rcmoy(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 166 rcmoy(ji,jj,15) = utaui_ice(ji,jj) 167 rcmoy(ji,jj,16) = vtaui_ice(ji,jj) 168 rcmoy(ji,jj,17) = qsr_ice(ji,jj) 169 rcmoy(ji,jj,18) = qns_ice(ji,jj) 170 170 rcmoy(ji,jj,19) = sprecip(ji,jj) 171 171 END DO … … 201 201 rcmoy(:,:,:) = 0.0 202 202 nmoyice = 0 203 END IF ! MOD(kt+n fice-1-nit000+1, nwrite == 0 ) !203 END IF ! MOD(kt+nn_fsbc-1-nit000+1, nwrite == 0 ) ! 204 204 205 205 END SUBROUTINE lim_wri_2
Note: See TracChangeset
for help on using the changeset viewer.