Commit 0ffc2cb7 authored by Miguel Tremblay's avatar Miguel Tremblay
Browse files

Removal of unused thermodynamical functions and their variables.

parent f0e8121b
......@@ -188,9 +188,8 @@
$ KNAMS, STLO, KARMAN, RIC, CHLC, CHLF, DG2RAD
REAL TTT,PRS,QQQ,EEE,TVI,QST
REAL FOEW,FODLE,FOQST,FODQS,FOEFQ,FOQFE,FOTVT,FOTTV,FOHR
REAL FOEWA,FODLA,FOQSA,FODQA,FOHRA
REAL TTT,QQQ
REAL FOTVT
LOGICAL bSilent
COMMON /SILENT/ bSilent
......@@ -205,52 +204,16 @@
* * REVISION 01 - MAI 94 - N. BRUNET
* NOUVELLE VERSION POUR FAIBLES PRESSIONS
*
* FONCTION DE TENSION DE VAPEUR SATURANTE (TETENS) - EW OU EI SELON TT
FOEW(TTT) = 610.78*EXP( AMIN1(SIGN(17.269,TTT-TRPL),SIGN
W (21.875,TTT-TRPL))*ABS(TTT-TRPL)/(TTT-35.86+AMAX1(0.,SIGN
W (28.2,TRPL-TTT))))
*
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW (OU LN EI)
FODLE(TTT)=(4097.93+AMAX1(0.,SIGN(1709.88,TRPL-TTT)))
W /((TTT-35.86+AMAX1(0.,SIGN(28.2,TRPL-TTT)))*(TTT-35.86+AMAX1(0.
W ,SIGN(28.2,TRPL-TTT))))
*
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE (QSAT)
FOQST(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEW(TTT))-EPS2)
*
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQS(QST,TTT)=QST*(1.+DELTA*QST)*FODLE(TTT)
* QST EST LA SORTIE DE FOQST
*
* FONCTION CALCULANT TENSION VAP (EEE) FN DE HUM SP (QQQ) ET PRS
FOEFQ(QQQ,PRS) = AMIN1(PRS,(QQQ*PRS) / (EPS1 + EPS2*QQQ))
*
* FONCTION CALCULANT HUM SP (QQQ) DE TENS. VAP (EEE) ET PRES (PRS)
FOQFE(EEE,PRS) = AMIN1(1.,EPS1*EEE / (PRS-EPS2*EEE))
*
* FONCTION CALCULANT TEMP VIRT. (TVI) DE TEMP (TTT) ET HUM SP (QQQ)
FOTVT(TTT,QQQ) = TTT * (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT TTT DE TEMP VIRT. (TVI) ET HUM SP (QQQ)
FOTTV(TVI,QQQ) = TVI / (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT HUM REL DE HUM SP (QQQ), TEMP (TTT) ET PRES (PRS)
* HR = E/ESAT
FOHR(QQQ,TTT,PRS) = AMIN1(PRS,FOEFQ(QQQ,PRS)) / FOEW(TTT)
*
* LES 5 FONCTIONS SUIVANTES SONT VALIDES DANS LE CONTEXTE OU ON
* NE DESIRE PAS TENIR COMPTE DE LA PHASE GLACE DANS LES CALCULS
* DE SATURATION.
* FONCTION DE VAPEUR SATURANTE (TETENS)
FOEWA(TTT)=610.78*EXP(17.269*(TTT-TRPL)/(TTT-35.86))
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW
FODLA(TTT)=17.269*(TRPL-35.86)/(TTT-35.86)**2
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE
FOQSA(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEWA(TTT))-EPS2)
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQA(QST,TTT)=QST*(1.+DELTA*QST)*FODLA(TTT)
* FONCTION CALCULANT L'HUMIDITE RELATIVE
FOHRA(QQQ,TTT,PRS)=AMIN1(PRS,FOEFQ(QQQ,PRS))/FOEWA(TTT)
*
***
*
......
......@@ -74,69 +74,8 @@
REAL PADDING
COMMON /BUFFER_SIZE/ DTMAX, Nl, DT, PADDING, n
REAL TTT,PRS,QQQ,EEE,TVI,QST
REAL FOEW,FODLE,FOQST,FODQS,FOEFQ,FOQFE,FOTVT,FOTTV,FOHR
REAL FOEWA,FODLA,FOQSA,FODQA,FOHRA
INTRINSIC SIGN
***
*
* DEFINITION DES FONCTIONS THERMODYNAMIQUES DE BASE
* POUR LES CONSTANTES, UTILISER LE COMMON /CTESDYN/
* NOTE: TOUTES LES FONCTIONS TRAVAILLENT AVEC LES UNITES S.I.
* I.E. TTT EN DEG K, PRS EN PA, QQQ EN KG/KG
* *** N. BRUNET - MAI 90 ***
* * REVISION 01 - MAI 94 - N. BRUNET
* NOUVELLE VERSION POUR FAIBLES PRESSIONS
*
* FONCTION DE TENSION DE VAPEUR SATURANTE (TETENS) - EW OU EI SELON TT
FOEW(TTT) = 610.78*EXP( AMIN1(SIGN(17.269,TTT-TRPL),SIGN
W (21.875,TTT-TRPL))*ABS(TTT-TRPL)/(TTT-35.86+AMAX1(0.,SIGN
W (28.2,TRPL-TTT))))
*
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW (OU LN EI)
FODLE(TTT)=(4097.93+AMAX1(0.,SIGN(1709.88,TRPL-TTT)))
W /((TTT-35.86+AMAX1(0.,SIGN(28.2,TRPL-TTT)))*(TTT-35.86+AMAX1(0.
W ,SIGN(28.2,TRPL-TTT))))
*
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE (QSAT)
FOQST(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEW(TTT))-EPS2)
*
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQS(QST,TTT)=QST*(1.+DELTA*QST)*FODLE(TTT)
* QST EST LA SORTIE DE FOQST
*
* FONCTION CALCULANT TENSION VAP (EEE) FN DE HUM SP (QQQ) ET PRS
FOEFQ(QQQ,PRS) = AMIN1(PRS,(QQQ*PRS) / (EPS1 + EPS2*QQQ))
*
* FONCTION CALCULANT HUM SP (QQQ) DE TENS. VAP (EEE) ET PRES (PRS)
FOQFE(EEE,PRS) = AMIN1(1.,EPS1*EEE / (PRS-EPS2*EEE))
*
* FONCTION CALCULANT TEMP VIRT. (TVI) DE TEMP (TTT) ET HUM SP (QQQ)
FOTVT(TTT,QQQ) = TTT * (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT TTT DE TEMP VIRT. (TVI) ET HUM SP (QQQ)
FOTTV(TVI,QQQ) = TVI / (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT HUM REL DE HUM SP (QQQ), TEMP (TTT) ET PRES (PRS)
* HR = E/ESAT
FOHR(QQQ,TTT,PRS) = AMIN1(PRS,FOEFQ(QQQ,PRS)) / FOEW(TTT)
*
* LES 5 FONCTIONS SUIVANTES SONT VALIDES DANS LE CONTEXTE OU ON
* NE DESIRE PAS TENIR COMPTE DE LA PHASE GLACE DANS LES CALCULS
* DE SATURATION.
* FONCTION DE VAPEUR SATURANTE (TETENS)
FOEWA(TTT)=610.78*EXP(17.269*(TTT-TRPL)/(TTT-35.86))
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW
FODLA(TTT)=17.269*(TRPL-35.86)/(TTT-35.86)**2
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE
FOQSA(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEWA(TTT))-EPS2)
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQA(QST,TTT)=QST*(1.+DELTA*QST)*FODLA(TTT)
* FONCTION CALCULANT L'HUMIDITE RELATIVE
FOHRA(QQQ,TTT,PRS)=AMIN1(PRS,FOEFQ(QQQ,PRS))/FOEWA(TTT)
*
if(silent) then
bSilent = .true.
......
......@@ -177,10 +177,8 @@
$ KNAMS, STLO, KARMAN, RIC, CHLC, CHLF, DG2RAD
REAL TTT,PRS,QQQ,EEE,TVI,QST
REAL FOEW,FODLE,FOQST,FODQS,FOEFQ,FOQFE,FOTVT,FOTTV,FOHR
REAL FOEWA,FODLA,FOQSA,FODQA,FOHRA
REAL TTT,QQQ
REAL FOTVT
LOGICAL bSilent
COMMON /SILENT/ bSilent
......@@ -194,52 +192,9 @@
* * REVISION 01 - MAI 94 - N. BRUNET
* NOUVELLE VERSION POUR FAIBLES PRESSIONS
*
* FONCTION DE TENSION DE VAPEUR SATURANTE (TETENS) - EW OU EI SELON TT
FOEW(TTT) = 610.78*EXP( AMIN1(SIGN(17.269,TTT-TRPL),SIGN
W (21.875,TTT-TRPL))*ABS(TTT-TRPL)/(TTT-35.86+AMAX1(0.,SIGN
W (28.2,TRPL-TTT))))
*
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW (OU LN EI)
FODLE(TTT)=(4097.93+AMAX1(0.,SIGN(1709.88,TRPL-TTT)))
W /((TTT-35.86+AMAX1(0.,SIGN(28.2,TRPL-TTT)))*(TTT-35.86+AMAX1(0.
W ,SIGN(28.2,TRPL-TTT))))
*
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE (QSAT)
FOQST(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEW(TTT))-EPS2)
*
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQS(QST,TTT)=QST*(1.+DELTA*QST)*FODLE(TTT)
* QST EST LA SORTIE DE FOQST
*
* FONCTION CALCULANT TENSION VAP (EEE) FN DE HUM SP (QQQ) ET PRS
FOEFQ(QQQ,PRS) = AMIN1(PRS,(QQQ*PRS) / (EPS1 + EPS2*QQQ))
*
* FONCTION CALCULANT HUM SP (QQQ) DE TENS. VAP (EEE) ET PRES (PRS)
FOQFE(EEE,PRS) = AMIN1(1.,EPS1*EEE / (PRS-EPS2*EEE))
*
* FONCTION CALCULANT TEMP VIRT. (TVI) DE TEMP (TTT) ET HUM SP (QQQ)
FOTVT(TTT,QQQ) = TTT * (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT TTT DE TEMP VIRT. (TVI) ET HUM SP (QQQ)
FOTTV(TVI,QQQ) = TVI / (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT HUM REL DE HUM SP (QQQ), TEMP (TTT) ET PRES (PRS)
* HR = E/ESAT
FOHR(QQQ,TTT,PRS) = AMIN1(PRS,FOEFQ(QQQ,PRS)) / FOEW(TTT)
*
* LES 5 FONCTIONS SUIVANTES SONT VALIDES DANS LE CONTEXTE OU ON
* NE DESIRE PAS TENIR COMPTE DE LA PHASE GLACE DANS LES CALCULS
* DE SATURATION.
* FONCTION DE VAPEUR SATURANTE (TETENS)
FOEWA(TTT)=610.78*EXP(17.269*(TTT-TRPL)/(TTT-35.86))
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW
FODLA(TTT)=17.269*(TRPL-35.86)/(TTT-35.86)**2
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE
FOQSA(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEWA(TTT))-EPS2)
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQA(QST,TTT)=QST*(1.+DELTA*QST)*FODLA(TTT)
* FONCTION CALCULANT L'HUMIDITE RELATIVE
FOHRA(QQQ,TTT,PRS)=AMIN1(PRS,FOEFQ(QQQ,PRS))/FOEWA(TTT)
*
***
*
......
......@@ -433,9 +433,8 @@
$ KNAMS, STLO, KARMAN, RIC, CHLC, CHLF, DG2RAD
*
*
REAL TTT,PRS,QQQ,EEE,TVI,QST
REAL FOEW,FODLE,FOQST,FODQS,FOEFQ,FOQFE,FOTVT,FOTTV,FOHR
REAL FOEWA,FODLA,FOQSA,FODQA,FOHRA
REAL TTT,PRS
REAL FOEW,FOQST
*
INTRINSIC SIGN
......@@ -452,48 +451,10 @@
FOEW(TTT) = 610.78*EXP( AMIN1(SIGN(17.269,TTT-TRPL),SIGN
W (21.875,TTT-TRPL))*ABS(TTT-TRPL)/(TTT-35.86+AMAX1(0.,SIGN
W (28.2,TRPL-TTT))))
*
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW (OU LN EI)
FODLE(TTT)=(4097.93+AMAX1(0.,SIGN(1709.88,TRPL-TTT)))
W /((TTT-35.86+AMAX1(0.,SIGN(28.2,TRPL-TTT)))*(TTT-35.86+AMAX1(0.
W ,SIGN(28.2,TRPL-TTT))))
*
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE (QSAT)
FOQST(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEW(TTT))-EPS2)
*
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQS(QST,TTT)=QST*(1.+DELTA*QST)*FODLE(TTT)
* QST EST LA SORTIE DE FOQST
*
* FONCTION CALCULANT TENSION VAP (EEE) FN DE HUM SP (QQQ) ET PRS
FOEFQ(QQQ,PRS) = AMIN1(PRS,(QQQ*PRS) / (EPS1 + EPS2*QQQ))
*
* FONCTION CALCULANT HUM SP (QQQ) DE TENS. VAP (EEE) ET PRES (PRS)
FOQFE(EEE,PRS) = AMIN1(1.,EPS1*EEE / (PRS-EPS2*EEE))
*
* FONCTION CALCULANT TEMP VIRT. (TVI) DE TEMP (TTT) ET HUM SP (QQQ)
FOTVT(TTT,QQQ) = TTT * (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT TTT DE TEMP VIRT. (TVI) ET HUM SP (QQQ)
FOTTV(TVI,QQQ) = TVI / (1.0 + DELTA*QQQ)
*
* FONCTION CALCULANT HUM REL DE HUM SP (QQQ), TEMP (TTT) ET PRES (PRS)
* HR = E/ESAT
FOHR(QQQ,TTT,PRS) = AMIN1(PRS,FOEFQ(QQQ,PRS)) / FOEW(TTT)
*
* LES 5 FONCTIONS SUIVANTES SONT VALIDES DANS LE CONTEXTE OU ON
* NE DESIRE PAS TENIR COMPTE DE LA PHASE GLACE DANS LES CALCULS
* DE SATURATION.
* FONCTION DE VAPEUR SATURANTE (TETENS)
FOEWA(TTT)=610.78*EXP(17.269*(TTT-TRPL)/(TTT-35.86))
* FONCTION CALCULANT LA DERIVEE SELON T DE LN EW
FODLA(TTT)=17.269*(TRPL-35.86)/(TTT-35.86)**2
* FONCTION CALCULANT L'HUMIDITE SPECIFIQUE SATURANTE
FOQSA(TTT,PRS)=EPS1/(AMAX1(1.,PRS/FOEWA(TTT))-EPS2)
* FONCTION CALCULANT LA DERIVEE DE QSAT SELON T
FODQA(QST,TTT)=QST*(1.+DELTA*QST)*FODLA(TTT)
* FONCTION CALCULANT L'HUMIDITE RELATIVE
FOHRA(QQQ,TTT,PRS)=AMIN1(PRS,FOEFQ(QQQ,PRS))/FOEWA(TTT)
*
***
*
......
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