Changeset 2987 for branches

Timestamp:

2011-10-25T09:18:32+02:00 (13 years ago)

Author:

cetlod

Message:

branch dev_LOCEAN_2011:minor changes in PISCES

Location:

branches/2011/dev_LOCEAN_2011/NEMOGCM

Files:

• CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces (modified) (3 diffs)
• CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_pisces (modified) (3 diffs)
• NEMO/TOP_SRC/PISCES/p4zflx.F90 (modified) (3 diffs)

branches/2011/dev_LOCEAN_2011/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces

@@ -28 +28 @@
 !              !              !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
    sn_patm     = 'presatm'    ,     -1            , 'patm'     ,  .true.      , .true. ,   'yearly'  , ''       , ''
-   cn_dir      = './'      !  root directory for the location of the dynamical files
+   cn_dir      = './'     !  root directory for the location of the dynamical files
+!
+   ln_presatm  = .true.   ! constant atmopsheric pressure (F) or from a file (T)
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampisbio     !   biological parameters
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   nrdttrc    =  1       ! time step frequency for biology
-   wsbio      =  2.      ! POC sinking speed
-   xkmort     =  1.E-7   ! half saturation constant for mortality
-   ferat3     =  10.E-6  ! Fe/C in zooplankton 
-   wsbio2     =  30.     ! Big particles sinking speed
+   nrdttrc    =  1        ! time step frequency for biology
+   wsbio      =  2.       ! POC sinking speed
+   xkmort     =  1.E-7    ! half saturation constant for mortality
+   ferat3     =  10.E-6   ! Fe/C in zooplankton 
+   wsbio2     =  30.      ! Big particles sinking speed
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
@@ -67 +69 @@
    excret     =  0.05     ! excretion ratio of phytoplankton
    excret2    =  0.05     ! excretion ratio of diatoms
-   ln_newprod =  .FALSE.  ! Enable new parame. of production (T/F) 
+   ln_newprod =  .false.  ! Enable new parame. of production (T/F) 
    bresp      =  0.00333  ! Basal respiration rate
    chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton
@@ -158 +160 @@
    cn_dir      = './'      !  root directory for the location of the dynamical files
    ln_dust     =  .true.   ! boolean for dust input from the atmosphere
-   ln_river    =  .false.   ! boolean for river input of nutrients
+   ln_river    =  .true.   ! boolean for river input of nutrients
    ln_ndepo    =  .true.   ! boolean for atmospheric deposition of N
    ln_ironsed  =  .true.   ! boolean for Fe input from sediments

branches/2011/dev_LOCEAN_2011/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_pisces

@@ -28 +28 @@
 !              !              !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
    sn_patm     = 'presatm'    ,     -1            , 'patm'     ,  .true.      , .true. ,   'yearly'  , ''       , ''
-   cn_dir      = './'      !  root directory for the location of the dynamical files
+   cn_dir      = './'     !  root directory for the location of the dynamical files
+!
+   ln_presatm  = .true.   ! constant atmopsheric pressure (F) or from a file (T)
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampisbio     !   biological parameters
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   nrdttrc    =  4       ! time step frequency for biology
-   wsbio      =  2.      ! POC sinking speed
-   xkmort     =  1.E-7   ! half saturation constant for mortality
-   ferat3     =  10.E-6  ! Fe/C in zooplankton 
-   wsbio2     =  30.     ! Big particles sinking speed
+   nrdttrc    =  4        ! time step frequency for biology
+   wsbio      =  2.       ! POC sinking speed
+   xkmort     =  1.E-7    ! half saturation constant for mortality
+   ferat3     =  10.E-6   ! Fe/C in zooplankton 
+   wsbio2     =  30.      ! Big particles sinking speed
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
@@ -67 +69 @@
    excret     =  0.05     ! excretion ratio of phytoplankton
    excret2    =  0.05     ! excretion ratio of diatoms
-   ln_newprod =  .FALSE.  ! Enable new parame. of production (T/F) 
+   ln_newprod =  .false.  ! Enable new parame. of production (T/F) 
    bresp      =  0.00333  ! Basal respiration rate
    chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton
@@ -158 +160 @@
    cn_dir      = './'      !  root directory for the location of the dynamical files
    ln_dust     =  .true.   ! boolean for dust input from the atmosphere
-   ln_river    =  .false.   ! boolean for river input of nutrients
+   ln_river    =  .true.   ! boolean for river input of nutrients
    ln_ndepo    =  .true.   ! boolean for atmospheric deposition of N
    ln_ironsed  =  .true.   ! boolean for Fe input from sediments

branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zflx.F90

@@ -46 +46 @@
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: atcco2h, years
    INTEGER  :: nmaxrec, numco2
+
+   !                                         !!* nampisatm namelist (Atmospheric PRessure) *
+   LOGICAL, PUBLIC ::   ln_presatm = .true.  !: ref. pressure: global mean Patm (F) or a constant (F)
 
    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
       !!
-      NAMELIST/nampisatm/ sn_patm, cn_dir
+      NAMELIST/nampisatm/ ln_presatm, sn_patm, cn_dir
 
       !                                         ! -------------------- !
@@ -344 +347 @@
          READ  ( numnatp, nampisatm ) 
          !
-         ALLOCATE( sf_patm(1), STAT=ierr )           !* allocate and fill sf_patm (forcing structure) with sn_patm
-         IF( ierr > 0 )   CALL ctl_stop( 'STOP', 'p4z_flx: unable to allocate sf_patm structure' )
-         !
-         CALL fld_fill( sf_patm, (/ sn_patm /), cn_dir, 'p4z_flx', 'Atmospheric pressure ', 'nampisatm' )
-                                ALLOCATE( sf_patm(1)%fnow(jpi,jpj,1)   )
-         IF( sn_patm%ln_tint )  ALLOCATE( sf_patm(1)%fdta(jpi,jpj,1,2) )
-         !
-      ENDIF
-      !
-      CALL fld_read( kt, 1, sf_patm )               !* input Patm provided at kt + 1/2
-      patm(:,:) = sf_patm(1)%fnow(:,:,1)                        ! atmospheric pressure
-
+         !
+         IF(lwp) THEN                                 !* control print
+            WRITE(numout,*)
+            WRITE(numout,*) '   Namelist nampisatm : Atmospheric Pressure as external forcing'
+            WRITE(numout,*) '      constant atmopsheric pressure (F) or from a file (T)  ln_presatm = ', ln_presatm
+            WRITE(numout,*)
+         ENDIF
+         !
+         IF( ln_presatm ) THEN
+            ALLOCATE( sf_patm(1), STAT=ierr )           !* allocate and fill sf_patm (forcing structure) with sn_patm
+            IF( ierr > 0 )   CALL ctl_stop( 'STOP', 'p4z_flx: unable to allocate sf_patm structure' )
+            !
+            CALL fld_fill( sf_patm, (/ sn_patm /), cn_dir, 'p4z_flx', 'Atmospheric pressure ', 'nampisatm' )
+                                   ALLOCATE( sf_patm(1)%fnow(jpi,jpj,1)   )
+            IF( sn_patm%ln_tint )  ALLOCATE( sf_patm(1)%fdta(jpi,jpj,1,2) )
+         ENDIF
+         !                                         
+         IF( .NOT.ln_presatm )   patm(:,:) = 1.e0    ! Initialize patm if no reading from a file
+         !
+      ENDIF
+      !
+      IF( ln_presatm ) THEN
+         CALL fld_read( kt, 1, sf_patm )               !* input Patm provided at kt + 1/2
+         patm(:,:) = sf_patm(1)%fnow(:,:,1)                        ! atmospheric pressure
+      ENDIF
+      !
    END SUBROUTINE p4z_patm