Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/conf_gcm.F90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/conf_gcm.F90
r5113 r5116 102 102 ! adjust=y not implemented in case of OpenMP threads... 103 103 !$OMP PARALLEL 104 IF ((OMP_GET_NUM_THREADS()>1).and.adjust) then105 write(lunout,*)'conf_gcm: Error, adjust should be set to n' &104 IF ((OMP_GET_NUM_THREADS()>1).and.adjust) THEN 105 WRITE(lunout,*)'conf_gcm: Error, adjust should be set to n' & 106 106 ,' when running with OpenMP threads' 107 107 abort_message = 'Wrong value for adjust' … … 436 436 ! ......... ( modif le 17/04/96 ) ......... 437 437 438 test_etatinit: IF (.not. etatinit) then438 test_etatinit: IF (.not. etatinit) THEN 439 439 !Config Key = clon 440 440 !Config Desc = centre du zoom, longitude … … 454 454 455 455 IF( ABS(clat - clatt)>= 0.001 ) THEN 456 write(lunout,*)'conf_gcm: La valeur de clat passee par run.def', &456 WRITE(lunout,*)'conf_gcm: La valeur de clat passee par run.def', & 457 457 ' est differente de celle lue sur le fichier start ' 458 458 CALL abort_gcm("conf_gcm","stopped",1) … … 468 468 469 469 IF( ABS(grossismx - grossismxx)>= 0.001 ) THEN 470 write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', &470 WRITE(lunout,*)'conf_gcm: La valeur de grossismx passee par ', & 471 471 'run.def est differente de celle lue sur le fichier start ' 472 472 CALL abort_gcm("conf_gcm","stopped",1) … … 482 482 483 483 IF( ABS(grossismy - grossismyy)>= 0.001 ) THEN 484 write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', &484 WRITE(lunout,*)'conf_gcm: La valeur de grossismy passee par ', & 485 485 'run.def est differente de celle lue sur le fichier start ' 486 486 CALL abort_gcm("conf_gcm","stopped",1) … … 488 488 489 489 IF( grossismx<1. ) THEN 490 write(lunout,*) &490 WRITE(lunout,*) & 491 491 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' 492 492 CALL abort_gcm("conf_gcm","stopped",1) … … 496 496 497 497 IF( grossismy<1. ) THEN 498 write(lunout,*) &498 WRITE(lunout,*) & 499 499 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' 500 500 CALL abort_gcm("conf_gcm","stopped",1) … … 503 503 ENDIF 504 504 505 write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay505 WRITE(lunout,*)'conf_gcm: alphax alphay',alphax,alphay 506 506 507 507 ! alphax et alphay sont les anciennes formulat. des grossissements … … 517 517 IF( .NOT.fxyhypb ) THEN 518 518 IF( fxyhypbb ) THEN 519 write(lunout,*)' ******** PBS DANS CONF_GCM ******** '520 write(lunout,*)' *** fxyhypb lu sur le fichier start est ', &519 WRITE(lunout,*)' ******** PBS DANS CONF_GCM ******** ' 520 WRITE(lunout,*)' *** fxyhypb lu sur le fichier start est ', & 521 521 'F alors qu il est T sur run.def ***' 522 522 CALL abort_gcm("conf_gcm","stopped",1) … … 524 524 ELSE 525 525 IF( .NOT.fxyhypbb ) THEN 526 write(lunout,*)' ******** PBS DANS CONF_GCM ******** '527 write(lunout,*)' *** fxyhypb lu sur le fichier start est ', &526 WRITE(lunout,*)' ******** PBS DANS CONF_GCM ******** ' 527 WRITE(lunout,*)' *** fxyhypb lu sur le fichier start est ', & 528 528 'T alors qu il est F sur run.def **** ' 529 529 CALL abort_gcm("conf_gcm","stopped",1) … … 541 541 IF( fxyhypb ) THEN 542 542 IF( ABS(dzoomx - dzoomxx)>= 0.001 ) THEN 543 write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ', &543 WRITE(lunout,*)'conf_gcm: La valeur de dzoomx passee par ', & 544 544 'run.def est differente de celle lue sur le fichier start ' 545 545 CALL abort_gcm("conf_gcm","stopped",1) … … 557 557 IF( fxyhypb ) THEN 558 558 IF( ABS(dzoomy - dzoomyy)>= 0.001 ) THEN 559 write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ', &559 WRITE(lunout,*)'conf_gcm: La valeur de dzoomy passee par ', & 560 560 'run.def est differente de celle lue sur le fichier start ' 561 561 CALL abort_gcm("conf_gcm","stopped",1) … … 572 572 IF( fxyhypb ) THEN 573 573 IF( ABS(taux - tauxx)>= 0.001 ) THEN 574 write(lunout,*)'conf_gcm: La valeur de taux passee par ', &574 WRITE(lunout,*)'conf_gcm: La valeur de taux passee par ', & 575 575 'run.def est differente de celle lue sur le fichier start ' 576 576 CALL abort_gcm("conf_gcm","stopped",1) … … 587 587 IF( fxyhypb ) THEN 588 588 IF( ABS(tauy - tauyy)>= 0.001 ) THEN 589 write(lunout,*)'conf_gcm: La valeur de tauy passee par ', &589 WRITE(lunout,*)'conf_gcm: La valeur de tauy passee par ', & 590 590 'run.def est differente de celle lue sur le fichier start ' 591 591 CALL abort_gcm("conf_gcm","stopped",1) … … 607 607 IF( .NOT.ysinus ) THEN 608 608 IF( ysinuss ) THEN 609 write(lunout,*)' ******** PBS DANS CONF_GCM ******** '610 write(lunout,*)' *** ysinus lu sur le fichier start est F', &609 WRITE(lunout,*)' ******** PBS DANS CONF_GCM ******** ' 610 WRITE(lunout,*)' *** ysinus lu sur le fichier start est F', & 611 611 ' alors qu il est T sur run.def ***' 612 612 CALL abort_gcm("conf_gcm","stopped",1) … … 614 614 ELSE 615 615 IF( .NOT.ysinuss ) THEN 616 write(lunout,*)' ******** PBS DANS CONF_GCM ******** '617 write(lunout,*)' *** ysinus lu sur le fichier start est T', &616 WRITE(lunout,*)' ******** PBS DANS CONF_GCM ******** ' 617 WRITE(lunout,*)' *** ysinus lu sur le fichier start est T', & 618 618 ' alors qu il est F sur run.def **** ' 619 619 CALL abort_gcm("conf_gcm","stopped",1) … … 697 697 CALL getin('ok_dyn_xios',ok_dyn_xios) 698 698 699 write(lunout,*)' #########################################'700 write(lunout,*)' Configuration des parametres du gcm: '701 write(lunout,*)' planet_type = ', planet_type702 write(lunout,*)' calend = ', calend703 write(lunout,*)' dayref = ', dayref704 write(lunout,*)' anneeref = ', anneeref705 write(lunout,*)' nday = ', nday706 write(lunout,*)' day_step = ', day_step707 write(lunout,*)' iperiod = ', iperiod708 write(lunout,*)' nsplit_phys = ', nsplit_phys709 write(lunout,*)' iconser = ', iconser710 write(lunout,*)' iecri = ', iecri711 write(lunout,*)' periodav = ', periodav712 write(lunout,*)' output_grads_dyn = ', output_grads_dyn713 write(lunout,*)' dissip_period = ', dissip_period714 write(lunout,*)' lstardis = ', lstardis715 write(lunout,*)' nitergdiv = ', nitergdiv716 write(lunout,*)' nitergrot = ', nitergrot717 write(lunout,*)' niterh = ', niterh718 write(lunout,*)' tetagdiv = ', tetagdiv719 write(lunout,*)' tetagrot = ', tetagrot720 write(lunout,*)' tetatemp = ', tetatemp721 write(lunout,*)' coefdis = ', coefdis722 write(lunout,*)' purmats = ', purmats723 write(lunout,*)' read_start = ', read_start724 write(lunout,*)' iflag_phys = ', iflag_phys725 write(lunout,*)' iphysiq = ', iphysiq726 write(lunout,*)' clonn = ', clonn727 write(lunout,*)' clatt = ', clatt728 write(lunout,*)' grossismx = ', grossismx729 write(lunout,*)' grossismy = ', grossismy730 write(lunout,*)' fxyhypbb = ', fxyhypbb731 write(lunout,*)' dzoomxx = ', dzoomxx732 write(lunout,*)' dzoomy = ', dzoomyy733 write(lunout,*)' tauxx = ', tauxx734 write(lunout,*)' tauyy = ', tauyy735 write(lunout,*)' offline = ', offline736 write(lunout,*)' offline_time = ', offline_time737 write(lunout,*)' type_trac = ', type_trac738 write(lunout,*)' ok_dynzon = ', ok_dynzon739 write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins740 write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave741 write(lunout,*)' ok_dyn_xios = ', ok_dyn_xios742 write(lunout,*)' adv_qsat_liq = ', adv_qsat_liq699 WRITE(lunout,*)' #########################################' 700 WRITE(lunout,*)' Configuration des parametres du gcm: ' 701 WRITE(lunout,*)' planet_type = ', planet_type 702 WRITE(lunout,*)' calend = ', calend 703 WRITE(lunout,*)' dayref = ', dayref 704 WRITE(lunout,*)' anneeref = ', anneeref 705 WRITE(lunout,*)' nday = ', nday 706 WRITE(lunout,*)' day_step = ', day_step 707 WRITE(lunout,*)' iperiod = ', iperiod 708 WRITE(lunout,*)' nsplit_phys = ', nsplit_phys 709 WRITE(lunout,*)' iconser = ', iconser 710 WRITE(lunout,*)' iecri = ', iecri 711 WRITE(lunout,*)' periodav = ', periodav 712 WRITE(lunout,*)' output_grads_dyn = ', output_grads_dyn 713 WRITE(lunout,*)' dissip_period = ', dissip_period 714 WRITE(lunout,*)' lstardis = ', lstardis 715 WRITE(lunout,*)' nitergdiv = ', nitergdiv 716 WRITE(lunout,*)' nitergrot = ', nitergrot 717 WRITE(lunout,*)' niterh = ', niterh 718 WRITE(lunout,*)' tetagdiv = ', tetagdiv 719 WRITE(lunout,*)' tetagrot = ', tetagrot 720 WRITE(lunout,*)' tetatemp = ', tetatemp 721 WRITE(lunout,*)' coefdis = ', coefdis 722 WRITE(lunout,*)' purmats = ', purmats 723 WRITE(lunout,*)' read_start = ', read_start 724 WRITE(lunout,*)' iflag_phys = ', iflag_phys 725 WRITE(lunout,*)' iphysiq = ', iphysiq 726 WRITE(lunout,*)' clonn = ', clonn 727 WRITE(lunout,*)' clatt = ', clatt 728 WRITE(lunout,*)' grossismx = ', grossismx 729 WRITE(lunout,*)' grossismy = ', grossismy 730 WRITE(lunout,*)' fxyhypbb = ', fxyhypbb 731 WRITE(lunout,*)' dzoomxx = ', dzoomxx 732 WRITE(lunout,*)' dzoomy = ', dzoomyy 733 WRITE(lunout,*)' tauxx = ', tauxx 734 WRITE(lunout,*)' tauyy = ', tauyy 735 WRITE(lunout,*)' offline = ', offline 736 WRITE(lunout,*)' offline_time = ', offline_time 737 WRITE(lunout,*)' type_trac = ', type_trac 738 WRITE(lunout,*)' ok_dynzon = ', ok_dynzon 739 WRITE(lunout,*)' ok_dyn_ins = ', ok_dyn_ins 740 WRITE(lunout,*)' ok_dyn_ave = ', ok_dyn_ave 741 WRITE(lunout,*)' ok_dyn_xios = ', ok_dyn_xios 742 WRITE(lunout,*)' adv_qsat_liq = ', adv_qsat_liq 743 743 else 744 744 !Config Key = clon … … 775 775 776 776 IF( grossismx<1. ) THEN 777 write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismx < 1 . *** '777 WRITE(lunout,*) 'conf_gcm: ***ATTENTION !! grossismx < 1 . *** ' 778 778 CALL abort_gcm("conf_gcm","stopped",1) 779 779 ELSE … … 782 782 783 783 IF( grossismy<1. ) THEN 784 write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** '784 WRITE(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** ' 785 785 CALL abort_gcm("conf_gcm","stopped",1) 786 786 ELSE … … 788 788 ENDIF 789 789 790 write(lunout,*) 'conf_gcm: alphax alphay ',alphax,alphay790 WRITE(lunout,*) 'conf_gcm: alphax alphay ',alphax,alphay 791 791 792 792 ! alphax et alphay sont les anciennes formulat. des grossissements … … 907 907 CALL getin('use_filtre_fft',use_filtre_fft) 908 908 IF (use_filtre_fft .AND. grossismx /= 1.0) THEN 909 write(lunout,*)'WARNING !!! '910 write(lunout,*)"A zoom in longitude is not compatible", &909 WRITE(lunout,*)'WARNING !!! ' 910 WRITE(lunout,*)"A zoom in longitude is not compatible", & 911 911 " with the FFT filter ", & 912 912 "---> FFT filter deactivated" … … 970 970 CALL getin('read_orop',read_orop) 971 971 972 write(lunout,*)' #########################################'973 write(lunout,*)' Configuration des parametres de cel0_limit: '974 write(lunout,*)' planet_type = ', planet_type975 write(lunout,*)' calend = ', calend976 write(lunout,*)' dayref = ', dayref977 write(lunout,*)' anneeref = ', anneeref978 write(lunout,*)' nday = ', nday979 write(lunout,*)' day_step = ', day_step980 write(lunout,*)' iperiod = ', iperiod981 write(lunout,*)' iconser = ', iconser982 write(lunout,*)' iecri = ', iecri983 write(lunout,*)' periodav = ', periodav984 write(lunout,*)' output_grads_dyn = ', output_grads_dyn985 write(lunout,*)' dissip_period = ', dissip_period986 write(lunout,*)' lstardis = ', lstardis987 write(lunout,*)' nitergdiv = ', nitergdiv988 write(lunout,*)' nitergrot = ', nitergrot989 write(lunout,*)' niterh = ', niterh990 write(lunout,*)' tetagdiv = ', tetagdiv991 write(lunout,*)' tetagrot = ', tetagrot992 write(lunout,*)' tetatemp = ', tetatemp993 write(lunout,*)' coefdis = ', coefdis994 write(lunout,*)' purmats = ', purmats995 write(lunout,*)' read_start = ', read_start996 write(lunout,*)' iflag_phys = ', iflag_phys997 write(lunout,*)' iphysiq = ', iphysiq998 write(lunout,*)' clon = ', clon999 write(lunout,*)' clat = ', clat1000 write(lunout,*)' grossismx = ', grossismx1001 write(lunout,*)' grossismy = ', grossismy1002 write(lunout,*)' fxyhypb = ', fxyhypb1003 write(lunout,*)' dzoomx = ', dzoomx1004 write(lunout,*)' dzoomy = ', dzoomy1005 write(lunout,*)' taux = ', taux1006 write(lunout,*)' tauy = ', tauy1007 write(lunout,*)' offline = ', offline1008 write(lunout,*)' type_trac = ', type_trac1009 write(lunout,*)' ok_dynzon = ', ok_dynzon1010 write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins1011 write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave1012 write(lunout,*)' ok_dyn_xios = ', ok_dyn_xios1013 write(lunout,*)' use_filtre_fft = ', use_filtre_fft1014 write(lunout,*)' use_mpi_alloc = ', use_mpi_alloc1015 write(lunout,*)' ok_strato = ', ok_strato1016 write(lunout,*)' ok_gradsfile = ', ok_gradsfile1017 write(lunout,*)' ok_limit = ', ok_limit1018 write(lunout,*)' ok_etat0 = ', ok_etat01019 write(lunout,*)' ok_guide = ', ok_guide1020 write(lunout,*)' read_orop = ', read_orop972 WRITE(lunout,*)' #########################################' 973 WRITE(lunout,*)' Configuration des parametres de cel0_limit: ' 974 WRITE(lunout,*)' planet_type = ', planet_type 975 WRITE(lunout,*)' calend = ', calend 976 WRITE(lunout,*)' dayref = ', dayref 977 WRITE(lunout,*)' anneeref = ', anneeref 978 WRITE(lunout,*)' nday = ', nday 979 WRITE(lunout,*)' day_step = ', day_step 980 WRITE(lunout,*)' iperiod = ', iperiod 981 WRITE(lunout,*)' iconser = ', iconser 982 WRITE(lunout,*)' iecri = ', iecri 983 WRITE(lunout,*)' periodav = ', periodav 984 WRITE(lunout,*)' output_grads_dyn = ', output_grads_dyn 985 WRITE(lunout,*)' dissip_period = ', dissip_period 986 WRITE(lunout,*)' lstardis = ', lstardis 987 WRITE(lunout,*)' nitergdiv = ', nitergdiv 988 WRITE(lunout,*)' nitergrot = ', nitergrot 989 WRITE(lunout,*)' niterh = ', niterh 990 WRITE(lunout,*)' tetagdiv = ', tetagdiv 991 WRITE(lunout,*)' tetagrot = ', tetagrot 992 WRITE(lunout,*)' tetatemp = ', tetatemp 993 WRITE(lunout,*)' coefdis = ', coefdis 994 WRITE(lunout,*)' purmats = ', purmats 995 WRITE(lunout,*)' read_start = ', read_start 996 WRITE(lunout,*)' iflag_phys = ', iflag_phys 997 WRITE(lunout,*)' iphysiq = ', iphysiq 998 WRITE(lunout,*)' clon = ', clon 999 WRITE(lunout,*)' clat = ', clat 1000 WRITE(lunout,*)' grossismx = ', grossismx 1001 WRITE(lunout,*)' grossismy = ', grossismy 1002 WRITE(lunout,*)' fxyhypb = ', fxyhypb 1003 WRITE(lunout,*)' dzoomx = ', dzoomx 1004 WRITE(lunout,*)' dzoomy = ', dzoomy 1005 WRITE(lunout,*)' taux = ', taux 1006 WRITE(lunout,*)' tauy = ', tauy 1007 WRITE(lunout,*)' offline = ', offline 1008 WRITE(lunout,*)' type_trac = ', type_trac 1009 WRITE(lunout,*)' ok_dynzon = ', ok_dynzon 1010 WRITE(lunout,*)' ok_dyn_ins = ', ok_dyn_ins 1011 WRITE(lunout,*)' ok_dyn_ave = ', ok_dyn_ave 1012 WRITE(lunout,*)' ok_dyn_xios = ', ok_dyn_xios 1013 WRITE(lunout,*)' use_filtre_fft = ', use_filtre_fft 1014 WRITE(lunout,*)' use_mpi_alloc = ', use_mpi_alloc 1015 WRITE(lunout,*)' ok_strato = ', ok_strato 1016 WRITE(lunout,*)' ok_gradsfile = ', ok_gradsfile 1017 WRITE(lunout,*)' ok_limit = ', ok_limit 1018 WRITE(lunout,*)' ok_etat0 = ', ok_etat0 1019 WRITE(lunout,*)' ok_guide = ', ok_guide 1020 WRITE(lunout,*)' read_orop = ', read_orop 1021 1021 ENDIF test_etatinit 1022 1022
Note: See TracChangeset
for help on using the changeset viewer.