Commit eb3930d2 authored by dumoda01's avatar dumoda01

Ecriture dun fichier se sortie pour les simulations en batch

parent 66a36157
......@@ -38,6 +38,7 @@
! !USES:
! default: all is private.
use bio_var
use output
private
!
! !PUBLIC MEMBER FUNCTIONS:
......@@ -115,11 +116,8 @@
REALTYPE :: b0 = 0.0
REALTYPE :: vp = 1.5
REALTYPE :: alpha = 0.065
!CHG1 Intensité optimale (I_opt))
REALTYPE :: I_opt = 10.0
!CHG1 Paramere d'inibition (inib)
REALTYPE :: inib = 0.05
!CHG2 Paramere de flottaison du phytoplancton
REALTYPE :: theta = 0.0
REALTYPE :: w_pmin = -0.06
REALTYPE :: w_pmax = -0.38
......@@ -174,6 +172,7 @@
! !INPUT PARAMETERS:
integer, intent(in) :: namlst
character(len=*), intent(in) :: fname
character(len=20) :: pfile
integer, intent(in) :: unit
!
! !REVISION HISTORY:
......@@ -199,20 +198,35 @@
numcc=numc
! Print some parameter values in standard output
LEVEL3 'FASHAM parameters ...'
! and save them in a separate file [out_fn]_fasham.par
pfile = trim(out_fn) // '_fasham.par'
open(10,status='unknown',action='write',file=pfile)
LEVEL3 'FASHAM parameters saved in ', pfile
write(*,900) ' vp = ',vp
write(10,901) vp
write(*,900) ' alpha = ',alpha
write(10,901) alpha
write(*,900) ' inib = ',inib
write(10,901) inib
write(*,900) ' k1 = ',k1
write(10,901) k1
write(*,900) ' k2 = ',k2
write(10,901) k2
write(*,900) ' w_p = ',w_p
write(10,901) w_p
write(*,900) ' theta = ',theta
write(10,901) theta
write(*,900) ' w_pmin = ',w_pmin
write(10,901) w_pmin
write(*,900) ' w_pmax = ',w_pmax
write(10,901) w_pmax
write(*,900) ' aa = ',aa
write(10,901) aa
write(*,900) ' g2 = ',g2
write(10,901) g2
900 format (A,F8.5)
900 format (a,f8.5)
901 format (f8.5)
! Conversion from day to second
vp = vp /secs_pr_day
......
......@@ -26,6 +26,7 @@
! !USES:
! default: all is private.
use bio_var
use output
private
!
! !PUBLIC MEMBER FUNCTIONS:
......@@ -136,6 +137,7 @@
! !INPUT PARAMETERS:
integer, intent(in) :: namlst
character(len=*), intent(in) :: fname
character(len=20) :: pfile
integer, intent(in) :: unit
!
! !REVISION HISTORY:
......@@ -160,16 +162,27 @@
numcc=numc
! Print some parameter values in standard output
LEVEL3 'NPZD parameters ...'
! and save them in a separate file [out_fn]_npzd.par
pfile = trim(out_fn) // '_npzd.par'
open(11,status='unknown',action='write',file=pfile)
LEVEL3 'NPZD parameters saved in ', pfile
write(*,900) ' rmax = ',rmax
write(11,901) rmax
write(*,900) ' iv = ',iv
write(11,901) iv
write(*,900) ' alpha = ',alpha
write(11,901) alpha
write(*,900) ' inib = ',inib
write(11,901) inib
write(*,900) ' k1 = ',k1
write(11,901) k1
write(*,900) ' aa = ',aa
write(11,901) aa
write(*,900) ' g2 = ',g2
write(11,901) g2
900 format (A,F8.5)
900 format (a,f8.5)
901 format (f8.5)
! Conversion from day to second
rpn = rpn /secs_pr_day
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment