Changeset 473 for trunk/NEMO/OPA_SRC/TRA/trabbc.F90
- Timestamp:
- 2006-05-11T17:04:37+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/trabbc.F90
r457 r473 137 137 !!---------------------------------------------------------------------- 138 138 !! * Modules used 139 USE io ipsl139 USE iom 140 140 141 141 !! * local declarations 142 CHARACTER (len=32) :: clname143 142 INTEGER :: ji, jj ! dummy loop indices 144 INTEGER :: inum = 11 ! temporary logical unit 145 INTEGER :: itime ! temporary integers 146 REAL(wp) :: zdate0, zdt ! temporary scalars 147 REAL(wp), DIMENSION(1) :: zdept ! temporary workspace 148 REAL(wp), DIMENSION(jpidta,jpjdta) :: & 149 zlamt, zphit, zdta ! temporary workspace 143 INTEGER :: inum ! temporary logical unit 150 144 151 145 NAMELIST/nambbc/ngeo_flux, ngeo_flux_const … … 188 182 CASE ( 2 ) ! variable geothermal heat flux 189 183 ! read the geothermal fluxes in mW/m2 190 clname = 'geothermal_heating' 191 #if defined key_agrif 192 if ( .NOT. Agrif_Root() ) then 193 clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 194 endif 195 #endif 196 itime = 1 197 zlamt(:,:) = 0. 198 zphit(:,:) = 0. 199 IF(lwp) WRITE(numout,*) ' *** variable geothermal heat flux read in ', clname, ' file' 200 CALL restini( clname, jpidta, jpjdta, zlamt, zphit, 1, zdept , 'NONE', & 201 & itime, zdate0, zdt, inum, domain_id=nidom ) 202 CALL restget( inum, 'heatflow', jpidta, jpjdta, 1, itime, .FALSE., zdta ) 203 DO jj = 1, nlcj 204 DO ji = 1, nlci 205 qgh_trd(ji,jj) = zdta(mig(ji),mjg(jj)) 206 END DO 207 END DO 208 209 CALL restclo( inum ) 184 185 IF(lwp) WRITE(numout,*) ' *** variable geothermal heat flux' 186 CALL iom_open ( 'geothermal_heating.nc', inum ) 187 CALL iom_get ( inum, jpdom_data, 'heatflow', qgh_trd ) 188 CALL iom_close (inum) 189 210 190 qgh_trd(:,:) = qgh_trd(:,:) * 1.e-3 ! conversion in W/m2 211 191 212 192 CASE DEFAULT 213 IF(lwp) WRITE(numout,cform_err) 214 IF(lwp) WRITE(numout,*) ' bad flag value for ngeo_flux = ', ngeo_flux 215 nstop = nstop + 1 216 193 WRITE(ctmp1,*) ' bad flag value for ngeo_flux = ', ngeo_flux 194 CALL ctl_stop( ctmp1 ) 217 195 END SELECT 218 196
Note: See TracChangeset
for help on using the changeset viewer.