Changeset 2141 for LMDZ5/trunk/libf/dyn3dmem
- Timestamp:
- Nov 10, 2014, 6:32:56 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F
r2111 r2141 2 2 ! $Id$ 3 3 ! 4 c 5 c 4 ! 5 ! 6 6 SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 ) 7 c 7 ! 8 8 USE control_mod 9 9 #ifdef CPP_IOIPSL … … 21 21 22 22 IMPLICIT NONE 23 c-----------------------------------------------------------------------24 cAuteurs : L. Fairhead , P. Le Van .25 c 26 cArguments :27 c 28 ctapedef :29 cetatinit : = TRUE , on ne compare pas les valeurs des para-30 c-metres du zoom avec celles lues sur le fichier start .31 cclesphy0 : sortie .32 c 23 !----------------------------------------------------------------------- 24 ! Auteurs : L. Fairhead , P. Le Van . 25 ! 26 ! Arguments : 27 ! 28 ! tapedef : 29 ! etatinit : = TRUE , on ne compare pas les valeurs des para- 30 ! -metres du zoom avec celles lues sur le fichier start . 31 ! clesphy0 : sortie . 32 ! 33 33 LOGICAL etatinit 34 34 INTEGER tapedef … … 37 37 PARAMETER( longcles = 20 ) 38 38 REAL clesphy0( longcles ) 39 c 40 cDeclarations :41 c--------------39 ! 40 ! Declarations : 41 ! -------------- 42 42 #include "dimensions.h" 43 43 #include "paramet.h" … … 51 51 ! #include "clesphys.h" 52 52 #include "iniprint.h" 53 c 54 c 55 clocal:56 c------53 ! 54 ! 55 ! local: 56 ! ------ 57 57 58 58 CHARACTER ch1*72,ch2*72,ch3*72,ch4*12 … … 66 66 integer,external :: OMP_GET_NUM_THREADS 67 67 #endif 68 c 69 c-------------------------------------------------------------------70 c 71 c......... Version du 29/04/97 ..........72 c 73 cNouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,74 ctetatemp ajoutes pour la dissipation .75 c 76 cAutre parametre ajoute en fin de liste de tapedef : ** fxyhypb **77 c 78 cSi fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.79 cSinon , choix de fxynew , a derivee sinusoidale ..80 c 81 c...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou82 cLIMIT_LMD pour l'initialisation de start.dat (dic) et83 cde limit.dat ( dic) ...........84 cSinon etatinit = . FALSE .85 c 86 cDonc etatinit = .F. si on veut comparer les valeurs de grossismx ,87 cgrossismy,clon,clat, fxyhypb lues sur le fichier start avec88 ccelles passees par run.def , au debut du gcm, apres l'appel a89 clectba .90 cCes parmetres definissant entre autres la grille et doivent etre91 cpareils et coherents , sinon il y aura divergence du gcm .92 c 93 c-----------------------------------------------------------------------94 cinitialisations:95 c----------------96 68 ! 69 ! ------------------------------------------------------------------- 70 ! 71 ! ......... Version du 29/04/97 .......... 72 ! 73 ! Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot, 74 ! tetatemp ajoutes pour la dissipation . 75 ! 76 ! Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb ** 77 ! 78 ! Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb. 79 ! Sinon , choix de fxynew , a derivee sinusoidale .. 80 ! 81 ! ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou 82 ! LIMIT_LMD pour l'initialisation de start.dat (dic) et 83 ! de limit.dat ( dic) ........... 84 ! Sinon etatinit = . FALSE . 85 ! 86 ! Donc etatinit = .F. si on veut comparer les valeurs de grossismx , 87 ! grossismy,clon,clat, fxyhypb lues sur le fichier start avec 88 ! celles passees par run.def , au debut du gcm, apres l'appel a 89 ! lectba . 90 ! Ces parmetres definissant entre autres la grille et doivent etre 91 ! pareils et coherents , sinon il y aura divergence du gcm . 92 ! 93 !----------------------------------------------------------------------- 94 ! initialisations: 95 ! ---------------- 96 97 97 !Config Key = lunout 98 98 !Config Desc = unite de fichier pour les impressions … … 134 134 CALL getin('prt_level',prt_level) 135 135 136 c-----------------------------------------------------------------------137 cParametres de controle du run:138 c-----------------------------------------------------------------------136 !----------------------------------------------------------------------- 137 ! Parametres de controle du run: 138 !----------------------------------------------------------------------- 139 139 !Config Key = planet_type 140 140 !Config Desc = planet type ("earth", "mars", "venus", ...) … … 266 266 CALL getin('dissip_period',dissip_period) 267 267 268 ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...269 ccc268 !cc .... P. Le Van , modif le 29/04/97 .pour la dissipation ... 269 !cc 270 270 271 271 !Config Key = lstardis … … 409 409 CALL getin('iphysiq', iphysiq) 410 410 411 if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", 412 $ "iphysiq must be a multiple of iperiod", 1) 413 411 414 !Config Key = ip_ebil_dyn 412 415 !Config Desc = PRINT level for energy conserv. diag. … … 422 425 423 426 424 ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...425 c......... ( modif le 17/04/96 ) .........426 c 427 !cc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... 428 ! ......... ( modif le 17/04/96 ) ......... 429 ! 427 430 IF( etatinit ) GO TO 100 428 431 … … 443 446 CALL getin('clat',clatt) 444 447 445 c 446 c 448 ! 449 ! 447 450 IF( ABS(clat - clatt).GE. 0.001 ) THEN 448 451 write(lunout,*)'conf_gcm: La valeur de clat passee par run.def', … … 498 501 499 502 write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay 500 c 501 calphax et alphay sont les anciennes formulat. des grossissements502 c 503 c 503 ! 504 ! alphax et alphay sont les anciennes formulat. des grossissements 505 ! 506 ! 504 507 505 508 !Config Key = fxyhypb … … 526 529 ENDIF 527 530 ENDIF 528 c 531 ! 529 532 !Config Key = dzoomx 530 533 !Config Desc = extension en longitude … … 589 592 ENDIF 590 593 591 cc594 !c 592 595 IF( .NOT.fxyhypb ) THEN 593 596 … … 617 620 ENDIF 618 621 ENDIF ! of IF( .NOT.fxyhypb ) 619 c 622 ! 620 623 !Config Key = offline 621 624 !Config Desc = Nouvelle eau liquide … … 722 725 723 726 RETURN 724 c...............................................725 c 727 ! ............................................... 728 ! 726 729 100 CONTINUE 727 730 !Config Key = clon … … 767 770 768 771 IF( grossismy.LT.1. ) THEN 769 write(lunout,*) 770 & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' 772 write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** ' 771 773 STOP 772 774 ELSE … … 775 777 776 778 write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay 777 c 778 calphax et alphay sont les anciennes formulat. des grossissements779 c 780 c 779 ! 780 ! alphax et alphay sont les anciennes formulat. des grossissements 781 ! 782 ! 781 783 782 784 !Config Key = fxyhypb … … 826 828 ysinus = .TRUE. 827 829 CALL getin('ysinus',ysinus) 828 c 830 ! 829 831 !Config Key = offline 830 832 !Config Desc = Nouvelle eau liquide … … 891 893 !Config Def = false 892 894 !Config Help = permet d'activer l'utilisation des FFT pour effectuer 893 !Config le filtrage aux poles. 895 !Config le filtrage aux poles. 894 896 use_filtre_fft=.FALSE. 895 897 CALL getin('use_filtre_fft',use_filtre_fft) … … 1001 1003 write(lunout,*)' ok_limit = ', ok_limit 1002 1004 write(lunout,*)' ok_etat0 = ', ok_etat0 1003 c 1005 ! 1004 1006 RETURN 1005 1007 END
Note: See TracChangeset
for help on using the changeset viewer.