Changeset 2987 for branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO
- Timestamp:
- 2011-10-25T09:18:32+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zflx.F90
r2977 r2987 46 46 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: atcco2h, years 47 47 INTEGER :: nmaxrec, numco2 48 49 ! !!* nampisatm namelist (Atmospheric PRessure) * 50 LOGICAL, PUBLIC :: ln_presatm = .true. !: ref. pressure: global mean Patm (F) or a constant (F) 48 51 49 52 REAL(wp) , ALLOCATABLE, SAVE, DIMENSION(:,:) :: patm ! atmospheric pressure at kt [N/m2] … … 329 332 TYPE(FLD_N) :: sn_patm ! informations about the fields to be read 330 333 !! 331 NAMELIST/nampisatm/ sn_patm, cn_dir334 NAMELIST/nampisatm/ ln_presatm, sn_patm, cn_dir 332 335 333 336 ! ! -------------------- ! … … 344 347 READ ( numnatp, nampisatm ) 345 348 ! 346 ALLOCATE( sf_patm(1), STAT=ierr ) !* allocate and fill sf_patm (forcing structure) with sn_patm 347 IF( ierr > 0 ) CALL ctl_stop( 'STOP', 'p4z_flx: unable to allocate sf_patm structure' ) 348 ! 349 CALL fld_fill( sf_patm, (/ sn_patm /), cn_dir, 'p4z_flx', 'Atmospheric pressure ', 'nampisatm' ) 350 ALLOCATE( sf_patm(1)%fnow(jpi,jpj,1) ) 351 IF( sn_patm%ln_tint ) ALLOCATE( sf_patm(1)%fdta(jpi,jpj,1,2) ) 352 ! 353 ENDIF 354 ! 355 CALL fld_read( kt, 1, sf_patm ) !* input Patm provided at kt + 1/2 356 patm(:,:) = sf_patm(1)%fnow(:,:,1) ! atmospheric pressure 357 349 ! 350 IF(lwp) THEN !* control print 351 WRITE(numout,*) 352 WRITE(numout,*) ' Namelist nampisatm : Atmospheric Pressure as external forcing' 353 WRITE(numout,*) ' constant atmopsheric pressure (F) or from a file (T) ln_presatm = ', ln_presatm 354 WRITE(numout,*) 355 ENDIF 356 ! 357 IF( ln_presatm ) THEN 358 ALLOCATE( sf_patm(1), STAT=ierr ) !* allocate and fill sf_patm (forcing structure) with sn_patm 359 IF( ierr > 0 ) CALL ctl_stop( 'STOP', 'p4z_flx: unable to allocate sf_patm structure' ) 360 ! 361 CALL fld_fill( sf_patm, (/ sn_patm /), cn_dir, 'p4z_flx', 'Atmospheric pressure ', 'nampisatm' ) 362 ALLOCATE( sf_patm(1)%fnow(jpi,jpj,1) ) 363 IF( sn_patm%ln_tint ) ALLOCATE( sf_patm(1)%fdta(jpi,jpj,1,2) ) 364 ENDIF 365 ! 366 IF( .NOT.ln_presatm ) patm(:,:) = 1.e0 ! Initialize patm if no reading from a file 367 ! 368 ENDIF 369 ! 370 IF( ln_presatm ) THEN 371 CALL fld_read( kt, 1, sf_patm ) !* input Patm provided at kt + 1/2 372 patm(:,:) = sf_patm(1)%fnow(:,:,1) ! atmospheric pressure 373 ENDIF 374 ! 358 375 END SUBROUTINE p4z_patm 359 376
Note: See TracChangeset
for help on using the changeset viewer.