Changeset 5093 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common
- Timestamp:
- Jul 21, 2024, 1:07:18 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/dyn3d_common
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/comdissnew.h
r1952 r5093 5 5 ! ATTENTION!!!!: ce fichier include est compatible format fixe/format libre 6 6 ! veillez à n'utiliser que des ! pour les commentaires 7 ! et à bien positionner les & des lignes de continuation 7 ! et à bien positionner les & des lignes de continuation 8 8 ! (les placer en colonne 6 et en colonne 73) 9 9 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F
r5086 r5093 24 24 c 25 25 c======================================================================= 26 c Discretisation verticale en coordonn �e hybride (ou sigma)26 c Discretisation verticale en coordonnée hybride (ou sigma) 27 27 c 28 28 c======================================================================= … … 208 208 c Calcul au milieu des couches : 209 209 c WARNING : le choix de placer le milieu des couches au niveau de 210 c pression interm �diaire est arbitraire et pourrait etre modifi�.210 c pression intermédiaire est arbitraire et pourrait etre modifié. 211 211 c Le calcul du niveau pour la derniere couche 212 212 c (on met la meme distance (en log pression) entre P(llm) 213 213 c et P(llm -1) qu'entre P(llm-1) et P(llm-2) ) est 214 c Specifique. Ce choix est sp �cifi�ici ET dans exner_milieu.F214 c Specifique. Ce choix est spécifié ici ET dans exner_milieu.F 215 215 216 216 DO l = 1, llm-1 … … 284 284 c L'objectif est de calculer newsig telle que 285 285 c (1 -pa/preff)*exp(1-1./newsig**2)+(pa/preff)*newsig = sig 286 c Cela ne se r �soud pas analytiquement:287 c => on r �soud par iterration bourrine286 c Cela ne se résoud pas analytiquement: 287 c => on résoud par iterration bourrine 288 288 c ---------------------------------------------- 289 289 c Information : where exp(1-1./x**2) become << x … … 320 320 newsig=(X2+newsig)*0.5 321 321 end if 322 c Test : on arete lorsque on approxime sig � moins de 0.01 m pr�s322 c Test : on arete lorsque on approxime sig à moins de 0.01 m près 323 323 c (en pseudo altitude) : 324 324 IF(abs(10.*log(F))<1.E-5) goto 999 -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F
r5086 r5093 41 41 42 42 c CALCUL DE LA PRESSION DE SURFACE 43 c Les coefficients ap et bp sont pass �s en common44 c Calcul de la pression au sol en mb optimis �e pour43 c Les coefficients ap et bp sont passés en common 44 c Calcul de la pression au sol en mb optimisée pour 45 45 c la vectorialisation 46 46 … … 67 67 c RECONSTRUCTION DES CHAMPS CONTRAVARIANTS 68 68 c Le programme ppm3d travaille avec les composantes 69 c de vitesse et pas les flux, on doit donc passer de l'un �l'autre70 c Dans le m �me temps, on fait le changement d'orientation du vent en v69 c de vitesse et pas les flux, on doit donc passer de l'un à l'autre 70 c Dans le même temps, on fait le changement d'orientation du vent en v 71 71 do l=1,llm 72 72 do j=1,jjm … … 99 99 100 100 c INVERSION DES NIVEAUX 101 c le programme ppm3d travaille avec une 3 �me coordonn�e invers�e par rapport101 c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport 102 102 c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface 103 c On passe donc des niveaux du LMDZ �ceux de Lin103 c On passe donc des niveaux du LMDZ à ceux de Lin 104 104 105 105 do l=1,llm+1 -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.F
r5079 r5093 68 68 implicit none 69 69 70 c rajout de d �clarations70 c rajout de déclarations 71 71 c integer Jmax,kmax,ndt0,nstep,k,j,i,ic,l,js,jn,imh,iad,jad,krd 72 72 c integer iu,iiu,j2,jmr,js0,jt … … 769 769 DO k=1,NLAY 770 770 DO I=1,IMR 771 c j=1 c'est le p �le Sud, j=JNP c'est le p�le Nord771 c j=1 c'est le pôle Sud, j=JNP c'est le pôle Nord 772 772 Q(I, 2,k,IC) = Q(I, 1,k,IC) 773 773 Q(I,JMR,k,IC) = Q(I,JNP,k,IC) -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.F
r2622 r5093 28 28 29 29 30 c On impose une seule valeur au pôle Sud j=jjm+1=jjp1 30 c On impose une seule valeur au pôle Sud j=jjm+1=jjp1 31 31 sommemasses=0 32 32 sommemqs=0
Note: See TracChangeset
for help on using the changeset viewer.