Changeset 2641 for LMDZ5/branches/testing/libf/dyn3dpar
- Timestamp:
- Sep 29, 2016, 11:26:46 PM (8 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 deleted
- 50 edited
- 2 copied
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2596-2640
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/dyn3dpar/addfi_p.F
r1999 r2641 48 48 #include "dimensions.h" 49 49 #include "paramet.h" 50 #include "comconst.h"51 50 #include "comgeom.h" 52 #include "serre.h"53 51 c 54 52 c Arguments : -
LMDZ5/branches/testing/libf/dyn3dpar/advect_new_p.F
r1910 r2641 6 6 USE parallel_lmdz 7 7 USE write_field_p 8 USE comconst_mod, ONLY: daysec 9 USE logic_mod, ONLY: conser 10 8 11 IMPLICIT NONE 9 12 c======================================================================= … … 28 31 #include "dimensions.h" 29 32 #include "paramet.h" 30 #include "comconst.h"31 #include "comvert.h"32 33 #include "comgeom.h" 33 #include "logic.h"34 #include "ener.h"35 34 36 35 c Arguments: … … 117 116 118 117 uav(ij,l)=0.25*(ucov(ij,l)+ucov(ij-iip1,l)) 119 . 118 . +0.25*(ucov(ij+iip1,l)+ucov(ij,l)) 120 119 ENDDO 121 120 … … 250 249 DO l=1,llm 251 250 DO ij=ijb,ije-1 252 253 251 du(ij,l)=du(ij,l)+du2(ij,l)-du1(ij,l) 252 ENDDO 254 253 255 254 DO ij = ijb+iip1-1, ije, iip1 … … 266 265 DO l=1,llm 267 266 DO ij=ijb,ije 268 269 267 dv(ij,l)=dv(ij,l)+dv2(ij,l)-dv1(ij,l) 268 ENDDO 270 269 ENDDO 271 270 c$OMP END DO NOWAIT … … 276 275 DO l=1,llm 277 276 DO ij=ijb,ije 278 279 277 dteta(ij,l)=dteta(ij,l)+dteta2(ij,l)-dteta1(ij,l) 278 ENDDO 280 279 ENDDO 281 280 c$OMP END DO NOWAIT -
LMDZ5/branches/testing/libf/dyn3dpar/advect_p.F
r1910 r2641 5 5 USE parallel_lmdz 6 6 USE write_field_p 7 USE comconst_mod, ONLY: daysec 8 USE logic_mod, ONLY: conser 9 7 10 IMPLICIT NONE 8 11 c======================================================================= … … 27 30 #include "dimensions.h" 28 31 #include "paramet.h" 29 #include "comconst.h"30 #include "comvert.h"31 32 #include "comgeom.h" 32 #include "logic.h"33 #include "ener.h"34 33 35 34 c Arguments: … … 86 85 87 86 uav(ij,l)=0.25*(ucov(ij,l)+ucov(ij-iip1,l)) 88 . 87 . +0.25*(ucov(ij+iip1,l)+ucov(ij,l)) 89 88 ENDDO 90 89 -
LMDZ5/branches/testing/libf/dyn3dpar/advtrac_p.F90
r2298 r2641 19 19 USE infotrac, ONLY: nqtot, iadv 20 20 USE control_mod, ONLY: iapp_tracvl, day_step, planet_type 21 USE comconst_mod, ONLY: dtvr 21 22 IMPLICIT NONE 22 23 ! 23 24 include "dimensions.h" 24 25 include "paramet.h" 25 include "comconst.h"26 include "comvert.h"27 26 include "comdissip.h" 28 27 include "comgeom2.h" 29 include "logic.h"30 include "temps.h"31 include "ener.h"32 28 include "description.h" 33 29 -
LMDZ5/branches/testing/libf/dyn3dpar/bernoui_p.F
r1910 r2641 25 25 #include "dimensions.h" 26 26 #include "paramet.h" 27 #include "logic.h"28 27 c 29 28 c Arguments: -
LMDZ5/branches/testing/libf/dyn3dpar/bilan_dyn_p.F
r1910 r2641 17 17 use misc_mod 18 18 use write_field_p 19 USE comconst_mod, ONLY: cpp, pi 20 USE comvert_mod, ONLY: presnivs 21 USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn 22 19 23 IMPLICIT NONE 20 24 21 25 #include "dimensions.h" 22 26 #include "paramet.h" 23 #include "comconst.h"24 #include "comvert.h"25 27 #include "comgeom2.h" 26 #include "temps.h"27 28 #include "iniprint.h" 28 29 … … 202 203 203 204 ! if (i_sortie.eq.1) then 204 ! 205 ! file='dynzon' 205 206 ! if (mpi_rank==0) then 206 ! 207 ! call inigrads(ifile,1 207 208 ! s ,0.,180./pi,0.,0.,jjm,rlatv,-90.,90.,180./pi 208 209 ! s ,llm,presnivs,1. … … 568 569 do l=1,llm 569 570 Q_cum(:,jjb:jje,l,iQ)=Q_cum(:,jjb:jje,l,iQ) 570 . 571 . /masse_cum(:,jjb:jje,l) 571 572 enddo 572 573 !$OMP ENDDO NOWAIT -
LMDZ5/branches/testing/libf/dyn3dpar/caladvtrac_p.F
r1910 r2641 10 10 USE infotrac, ONLY : nqtot 11 11 USE control_mod, ONLY : iapp_tracvl,planet_type 12 USE comconst_mod, ONLY: dtvr 12 13 c 13 14 IMPLICIT NONE … … 25 26 #include "dimensions.h" 26 27 #include "paramet.h" 27 #include "comconst.h"28 28 29 29 c Arguments: -
LMDZ5/branches/testing/libf/dyn3dpar/caldyn_p.F
r1999 r2641 10 10 USE parallel_lmdz 11 11 USE Write_Field_p 12 USE comvert_mod, ONLY: ap, bp 12 13 13 14 IMPLICIT NONE … … 31 32 #include "dimensions.h" 32 33 #include "paramet.h" 33 #include "comconst.h"34 #include "comvert.h"35 34 #include "comgeom.h" 36 35 -
LMDZ5/branches/testing/libf/dyn3dpar/conf_gcm.F90
r2471 r2641 16 16 USE infotrac, ONLY : type_trac 17 17 use assert_m, only: assert 18 USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, & 19 iflag_top_bound, mode_top_bound, tau_top_bound, & 20 ngroup 21 USE logic_mod, ONLY: fxyhypb, iflag_phys, ok_etat0, ok_gradsfile, & 22 ok_guide, ok_limit, ok_strato, purmats, read_start, & 23 ysinus 24 USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, & 25 alphax,alphay,taux,tauy 26 USE temps_mod, ONLY: calend 18 27 19 28 IMPLICIT NONE … … 34 43 include "dimensions.h" 35 44 include "paramet.h" 36 include "logic.h"37 include "serre.h"38 45 include "comdissnew.h" 39 include "temps.h"40 include "comconst.h"41 46 include "iniprint.h" 42 47 -
LMDZ5/branches/testing/libf/dyn3dpar/convmas1_p.F
r1910 r2641 32 32 #include "dimensions.h" 33 33 #include "paramet.h" 34 #include "comvert.h"35 #include "logic.h"36 34 37 35 REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ) -
LMDZ5/branches/testing/libf/dyn3dpar/convmas2_p.F
r1910 r2641 32 32 #include "dimensions.h" 33 33 #include "paramet.h" 34 #include "comvert.h"35 #include "logic.h"36 34 37 35 REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ) -
LMDZ5/branches/testing/libf/dyn3dpar/convmas_p.F
r1910 r2641 32 32 #include "dimensions.h" 33 33 #include "paramet.h" 34 #include "comvert.h"35 #include "logic.h"36 34 37 35 REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ) -
LMDZ5/branches/testing/libf/dyn3dpar/dissip_p.F
r1999 r2641 6 6 USE parallel_lmdz 7 7 USE write_field_p 8 USE comconst_mod, ONLY: dtdiss 8 9 IMPLICIT NONE 9 10 … … 29 30 #include "dimensions.h" 30 31 #include "paramet.h" 31 #include "comconst.h"32 32 #include "comgeom.h" 33 33 #include "comdissnew.h" -
LMDZ5/branches/testing/libf/dyn3dpar/divgrad_p.F
r1910 r2641 23 23 #include "comgeom.h" 24 24 #include "comdissipn.h" 25 #include "logic.h"26 25 c 27 26 INTEGER klevel -
LMDZ5/branches/testing/libf/dyn3dpar/dteta1_p.F
r1910 r2641 22 22 #include "dimensions.h" 23 23 #include "paramet.h" 24 #include "logic.h"25 24 26 25 REAL teta( ip1jmp1,llm ),pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm) -
LMDZ5/branches/testing/libf/dyn3dpar/dudv2_p.F
r1910 r2641 25 25 #include "dimensions.h" 26 26 #include "paramet.h" 27 #include "comvert.h"28 27 29 28 REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ), -
LMDZ5/branches/testing/libf/dyn3dpar/dynetat0.F
r1999 r2641 9 9 10 10 use control_mod, only : planet_type 11 USE comvert_mod, ONLY: pa,preff 12 USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa, 13 & lllm, omeg, rad 14 USE logic_mod, ONLY: fxyhypb, ysinus 15 USE serre_mod, ONLY: clon,clat,grossismx,grossismy 16 USE temps_mod, ONLY: annee_ref,day_ref,itau_dyn, 17 & start_time,day_ini,hour_ini 18 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 11 19 12 20 IMPLICIT NONE … … 29 37 #include "dimensions.h" 30 38 #include "paramet.h" 31 #include "temps.h"32 #include "comconst.h"33 #include "comvert.h"34 39 #include "comgeom2.h" 35 #include "ener.h"36 40 #include "netcdf.inc" 37 41 #include "description.h" 38 #include "serre.h"39 #include "logic.h"40 42 #include "iniprint.h" 41 43 -
LMDZ5/branches/testing/libf/dyn3dpar/dynredem.F
r1999 r2641 9 9 USE infotrac 10 10 use netcdf95, only: NF95_PUT_VAR 11 USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff, 12 & nivsig,nivsigs 13 USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad 14 USE logic_mod, ONLY: fxyhypb, ysinus 15 USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, 16 & taux,tauy 17 USE temps_mod, ONLY: annee_ref,day_ref,itau_dyn,itaufin, 18 & start_time,hour_ini 19 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 11 20 12 21 IMPLICIT NONE … … 18 27 #include "dimensions.h" 19 28 #include "paramet.h" 20 #include "comconst.h"21 #include "comvert.h"22 29 #include "comgeom2.h" 23 #include "temps.h"24 #include "ener.h"25 #include "logic.h"26 30 #include "netcdf.inc" 27 31 #include "description.h" 28 #include "serre.h"29 32 #include "iniprint.h" 30 33 … … 471 474 use netcdf, only: NF90_get_VAR 472 475 use netcdf95, only: NF95_PUT_VAR 476 USE temps_mod, ONLY: itau_dyn, itaufin 473 477 474 478 IMPLICIT NONE … … 480 484 #include "description.h" 481 485 #include "netcdf.inc" 482 #include "comvert.h"483 486 #include "comgeom.h" 484 #include "temps.h"485 487 #include "iniprint.h" 486 488 -
LMDZ5/branches/testing/libf/dyn3dpar/dynredem_p.F
r1910 r2641 10 10 USE infotrac 11 11 use netcdf95, only: NF95_PUT_VAR 12 USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff, 13 & nivsig,nivsigs 14 USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad 15 USE logic_mod, ONLY: fxyhypb, ysinus 16 USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, 17 & taux,tauy 18 USE temps_mod, ONLY: annee_ref,day_ref,itau_dyn,itaufin, 19 & start_time,hour_ini 20 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 12 21 13 22 IMPLICIT NONE … … 19 28 #include "dimensions.h" 20 29 #include "paramet.h" 21 #include "comconst.h"22 #include "comvert.h"23 30 #include "comgeom2.h" 24 #include "temps.h"25 #include "ener.h"26 #include "logic.h"27 31 #include "netcdf.inc" 28 32 #include "description.h" 29 #include "serre.h"30 33 31 34 c Arguments: … … 472 475 use netcdf, only: NF90_get_VAR 473 476 use netcdf95, only: NF95_PUT_VAR 477 USE temps_mod, ONLY: itau_dyn, itaufin 474 478 475 479 IMPLICIT NONE … … 481 485 #include "description.h" 482 486 #include "netcdf.inc" 483 #include "comvert.h"484 487 #include "comgeom.h" 485 #include "temps.h"486 488 487 489 -
LMDZ5/branches/testing/libf/dyn3dpar/exner_hyb_p_m.F90
r2056 r2641 33 33 ! 34 34 USE parallel_lmdz 35 USE comconst_mod, ONLY: cpp, kappa, r, jmp1 36 USE comvert_mod, ONLY: preff 35 37 ! 36 38 include "dimensions.h" 37 39 include "paramet.h" 38 include "comconst.h"39 40 include "comgeom.h" 40 include "comvert.h"41 include "serre.h"42 41 43 42 INTEGER ngrid -
LMDZ5/branches/testing/libf/dyn3dpar/exner_milieu_p_m.F90
r2056 r2641 30 30 ! 31 31 USE parallel_lmdz 32 USE comconst_mod, ONLY: cpp, kappa, r, jmp1 33 USE comvert_mod, ONLY: preff 32 34 ! 33 35 include "dimensions.h" 34 36 include "paramet.h" 35 include "comconst.h"36 37 include "comgeom.h" 37 include "comvert.h"38 include "serre.h"39 38 40 39 INTEGER ngrid -
LMDZ5/branches/testing/libf/dyn3dpar/fluxstokenc_p.F
r1910 r2641 21 21 #include "dimensions.h" 22 22 #include "paramet.h" 23 #include "comconst.h"24 #include "comvert.h"25 23 #include "comgeom.h" 26 24 #include "tracstoke.h" 27 #include "temps.h"28 25 #include "iniprint.h" 29 26 -
LMDZ5/branches/testing/libf/dyn3dpar/friction_p.F
r1910 r2641 12 12 USE ioipsl_getincom 13 13 #endif 14 USE comconst_mod, ONLY: pi 14 15 IMPLICIT NONE 15 16 … … 29 30 #include "paramet.h" 30 31 #include "comgeom2.h" 31 #include "comconst.h"32 32 #include "iniprint.h" 33 33 #include "academic.h" -
LMDZ5/branches/testing/libf/dyn3dpar/gcm.F
r2435 r2641 32 32 USE iniphysiq_mod, ONLY: iniphysiq 33 33 #endif 34 USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, r, rad 35 USE logic_mod ! all of it, because of copyin clause when calling leapfrog 36 USE temps_mod, ONLY: calend,start_time,annee_ref,day_ref, 37 & itau_dyn,itau_phy,day_ini,jD_ref,jH_ref,day_end, 38 & dt,hour_ini,itaufin 39 34 40 IMPLICIT NONE 35 41 … … 65 71 #include "dimensions.h" 66 72 #include "paramet.h" 67 #include "comconst.h"68 73 #include "comdissnew.h" 69 #include "comvert.h"70 74 #include "comgeom.h" 71 #include "logic.h"72 #include "temps.h"73 #include "ener.h"74 75 #include "description.h" 75 #include "serre.h"76 76 !#include "com_io_dyn.h" 77 77 #include "iniprint.h" … … 499 499 c write(78,*) 'q',q 500 500 501 c$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/,/logici/,/logicl/) 501 !c$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/,/logici/,/logicl/) 502 c$OMP PARALLEL DEFAULT(SHARED) 503 c Copy all threadprivate variables from temps_mod 504 c$OMP1 COPYIN(dt,jD_ref,jH_ref,start_time,hour_ini,day_ini,day_end) 505 c$OMP1 COPYIN(annee_ref,day_ref,itau_dyn,itau_phy,itaufin,calend) 506 c Copy all threadprivate variables from logic_mod 507 c$OMP1 COPYIN(purmats,forward,leapf,apphys,statcl,conser,apdiss,apdelq) 508 c$OMP1 COPYIN(saison,ecripar,fxyhypb,ysinus,read_start,ok_guide) 509 c$OMP1 COPYIN(ok_strato,ok_gradsfile,ok_limit,ok_etat0) 510 c$OMP1 COPYIN(iflag_phys,iflag_trac) 502 511 CALL leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,time_0) 503 512 c$OMP END PARALLEL -
LMDZ5/branches/testing/libf/dyn3dpar/geopot_p.F
r1910 r2641 28 28 #include "dimensions.h" 29 29 #include "paramet.h" 30 #include "comvert.h"31 30 32 31 c Arguments: -
LMDZ5/branches/testing/libf/dyn3dpar/gradiv_p.F
r1910 r2641 21 21 #include "paramet.h" 22 22 #include "comdissipn.h" 23 #include "logic.h"24 23 25 24 INTEGER klevel -
LMDZ5/branches/testing/libf/dyn3dpar/groupe_p.F
r2471 r2641 1 1 subroutine groupe_p(pext,pbaru,pbarv,pbarum,pbarvm,wm) 2 2 USE parallel_lmdz 3 USE comconst_mod, ONLY: ngroup 3 4 implicit none 4 5 … … 16 17 #include "dimensions.h" 17 18 #include "paramet.h" 18 #include "comconst.h"19 19 #include "comgeom2.h" 20 #include "comvert.h"21 20 22 21 ! integer ngroup -
LMDZ5/branches/testing/libf/dyn3dpar/groupeun_p.F
r2471 r2641 2 2 USE parallel_lmdz 3 3 USE Write_Field_p 4 USE comconst_mod, ONLY: ngroup 4 5 IMPLICIT NONE 5 6 6 7 #include "dimensions.h" 7 8 #include "paramet.h" 8 #include "comconst.h"9 9 #include "comgeom2.h" 10 10 … … 136 136 137 137 USE parallel_lmdz 138 USE comconst_mod, ONLY: ngroup 138 139 IMPLICIT NONE 139 140 140 141 #include "dimensions.h" 141 142 #include "paramet.h" 142 #include "comconst.h"143 143 #include "comgeom2.h" 144 144 -
LMDZ5/branches/testing/libf/dyn3dpar/guide_p_mod.F90
r2160 r2641 68 68 SUBROUTINE guide_init 69 69 70 USE control_mod 70 USE control_mod, ONLY: day_step 71 USE serre_mod, ONLY: grossismx 71 72 72 73 IMPLICIT NONE … … 75 76 INCLUDE "paramet.h" 76 77 INCLUDE "netcdf.inc" 77 78 ! For grossismx:79 include "serre.h"80 78 81 79 INTEGER :: error,ncidpl,rid,rcod … … 341 339 USE parallel_lmdz 342 340 USE control_mod 341 USE comconst_mod, ONLY: daysec, dtvr, cpp, kappa 342 USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner 343 343 344 344 IMPLICIT NONE … … 346 346 INCLUDE "dimensions.h" 347 347 INCLUDE "paramet.h" 348 INCLUDE "comconst.h"349 INCLUDE "comvert.h"350 348 351 349 ! Variables entree … … 619 617 SUBROUTINE guide_zonave(typ,hsize,vsize,field) 620 618 619 USE comconst_mod, ONLY: pi 620 621 621 IMPLICIT NONE 622 622 … … 624 624 INCLUDE "paramet.h" 625 625 INCLUDE "comgeom.h" 626 INCLUDE "comconst.h"627 626 628 627 ! input/output variables … … 707 706 USE mod_hallo 708 707 USE Bands 708 USE comconst_mod, ONLY: cpp, kappa 709 USE comvert_mod, ONLY: preff, pressure_exner, bp, ap 709 710 IMPLICIT NONE 710 711 711 712 include "dimensions.h" 712 713 include "paramet.h" 713 include "comvert.h"714 714 include "comgeom2.h" 715 include "comconst.h"716 715 717 716 REAL, DIMENSION (iip1,jjp1), INTENT(IN) :: psi ! Psol gcm … … 1099 1098 ! Calcul des constantes de rappel alpha (=1/tau) 1100 1099 1100 use comconst_mod, only: pi 1101 use serre_mod, only: clat, clon, grossismx, grossismy 1102 1101 1103 implicit none 1102 1104 1103 1105 include "dimensions.h" 1104 1106 include "paramet.h" 1105 include "comconst.h"1106 1107 include "comgeom2.h" 1107 include "serre.h"1108 1108 1109 1109 ! input arguments : … … 1813 1813 SUBROUTINE guide_out(varname,hsize,vsize,field,factt) 1814 1814 USE parallel_lmdz 1815 USE comconst_mod, ONLY: pi 1816 USE comvert_mod, ONLY: presnivs 1815 1817 IMPLICIT NONE 1816 1818 … … 1819 1821 INCLUDE "netcdf.inc" 1820 1822 INCLUDE "comgeom2.h" 1821 INCLUDE "comconst.h"1822 INCLUDE "comvert.h"1823 1823 1824 1824 ! Variables entree -
LMDZ5/branches/testing/libf/dyn3dpar/iniacademic.F90
r2160 r2641 16 16 use exner_hyb_m, only: exner_hyb 17 17 use exner_milieu_m, only: exner_milieu 18 USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm 19 USE logic_mod, ONLY: iflag_phys, read_start 20 USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner 21 USE temps_mod, ONLY: annee_ref, day_ini, day_ref 22 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 18 23 19 24 ! Author: Frederic Hourdin original: 15/01/93 … … 28 33 include "dimensions.h" 29 34 include "paramet.h" 30 include "comvert.h"31 include "comconst.h"32 35 include "comgeom.h" 33 36 include "academic.h" 34 include "ener.h"35 include "temps.h"36 37 include "iniprint.h" 37 include "logic.h"38 38 39 39 ! Arguments: -
LMDZ5/branches/testing/libf/dyn3dpar/initdynav_p.F
r1910 r2641 12 12 use misc_mod 13 13 USE infotrac 14 USE comconst_mod, ONLY: pi 15 USE comvert_mod, ONLY: nivsigs 16 USE temps_mod, ONLY: itau_dyn 14 17 15 18 implicit none … … 43 46 #include "dimensions.h" 44 47 #include "paramet.h" 45 #include "comconst.h"46 #include "comvert.h"47 48 #include "comgeom.h" 48 #include "temps.h"49 #include "ener.h"50 #include "logic.h"51 49 #include "description.h" 52 #include "serre.h"53 50 #include "iniprint.h" 54 51 -
LMDZ5/branches/testing/libf/dyn3dpar/initfluxsto_p.F
r1910 r2641 13 13 use Write_field 14 14 use misc_mod 15 USE comconst_mod, ONLY: pi 16 USE comvert_mod, ONLY: nivsigs 17 USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn 15 18 16 19 implicit none … … 45 48 #include "dimensions.h" 46 49 #include "paramet.h" 47 #include "comconst.h"48 #include "comvert.h"49 50 #include "comgeom.h" 50 #include "temps.h"51 #include "ener.h"52 #include "logic.h"53 51 #include "description.h" 54 #include "serre.h"55 52 #include "iniprint.h" 56 53 … … 106 103 CALL ymds2ju(zan, 1, idayref, 0.0, zjulian) 107 104 tau0 = itau_dyn 108 109 105 106 do jj = 1, jjp1 110 107 do ii = 1, iip1 111 108 rlong(ii,jj) = rlonu(ii) * 180. / pi … … 165 162 . 1, iip1, 1, jjn,tau0, zjulian, tstep, vhoriid, 166 163 . filevid,dynv_domain_id) 167 168 rl(1,1) = 1. 164 165 rl(1,1) = 1. 169 166 170 167 if (mpi_rank==0) then … … 191 188 call histhori(fileid, iip1, rlong(:,jjb:jje),jjn,rlat(:,jjb:jje), 192 189 . 'scalar','Grille points scalaires', thoriid) 193 190 194 191 C 195 192 C Appel a histvert pour la grille verticale … … 211 208 C 212 209 C Appels a histdef pour la definition des variables a sauvegarder 213 214 210 211 CALL histdef(fileid, "phis", "Surface geop. height", "-", 215 212 . iip1,jjn,thoriid, 1,1,1, -99, 32, 216 213 . "once", t_ops, t_wrt) … … 219 216 . iip1,jjn,thoriid, 1,1,1, -99, 32, 220 217 . "once", t_ops, t_wrt) 221 218 222 219 if (mpi_rank==0) then 223 224 220 221 CALL histdef(filedid, "dtvr", "tps dyn", "s", 225 222 . 1,1,dhoriid, 1,1,1, -99, 32, 226 223 . "once", t_ops, t_wrt) … … 289 286 if (mpi_rank==0) call histsync(filedid) 290 287 endif 291 288 292 289 #else 293 290 write(lunout,*)'initfluxsto_p: Needs IOIPSL to function' -
LMDZ5/branches/testing/libf/dyn3dpar/inithist_p.F
r1910 r2641 13 13 use misc_mod 14 14 USE infotrac 15 USE comconst_mod, ONLY: pi 16 USE comvert_mod, ONLY: nivsigs 17 USE temps_mod, ONLY: itau_dyn 15 18 16 19 implicit none … … 45 48 #include "dimensions.h" 46 49 #include "paramet.h" 47 #include "comconst.h"48 #include "comvert.h"49 50 #include "comgeom.h" 50 #include "temps.h"51 #include "ener.h"52 #include "logic.h"53 51 #include "description.h" 54 #include "serre.h"55 52 #include "iniprint.h" 56 53 -
LMDZ5/branches/testing/libf/dyn3dpar/integrd_p.F
r2160 r2641 7 7 USE parallel_lmdz 8 8 USE control_mod, only : planet_type 9 USE comconst_mod, ONLY: pi 10 USE logic_mod, ONLY: leapf 11 USE comvert_mod, ONLY: ap, bp 12 USE temps_mod, ONLY: dt 13 9 14 IMPLICIT NONE 10 15 … … 27 32 #include "dimensions.h" 28 33 #include "paramet.h" 29 #include "comconst.h"30 34 #include "comgeom.h" 31 #include "comvert.h"32 #include "logic.h"33 #include "temps.h"34 #include "serre.h"35 35 #include "iniprint.h" 36 36 … … 284 284 ije=ij_end 285 285 286 286 if (planet_type.eq."earth") then 287 287 ! Earth-specific treatment of first 2 tracers (water) 288 288 c$OMP BARRIER … … 385 385 DO l = 1, llm 386 386 massem1(ijb:ije,l)=massescr(ijb:ije,l) 387 388 c$OMP END DO NOWAIT 387 ENDDO 388 c$OMP END DO NOWAIT 389 389 END IF 390 390 c$OMP BARRIER -
LMDZ5/branches/testing/libf/dyn3dpar/leapfrog_p.F
r2408 r2641 28 28 & periodav, ok_dyn_ave, output_grads_dyn, 29 29 & iapp_tracvl 30 USE comvert_mod, ONLY: ap,bp,pressure_exner,presnivs 31 USE comconst_mod, ONLY: cpp, dtvr, ihf, dtphys, pi, jmp1 32 USE logic_mod, ONLY: iflag_phys,ok_guide,forward,leapf,apphys, 33 & statcl,conser,apdiss,purmats,ok_strato 34 USE temps_mod, ONLY: itaufin,jD_ref,jH_ref,day_ini, 35 & day_ref,start_time,dt 36 30 37 IMPLICIT NONE 31 38 … … 63 70 #include "dimensions.h" 64 71 #include "paramet.h" 65 #include "comconst.h"66 72 #include "comdissnew.h" 67 #include "comvert.h"68 73 #include "comgeom.h" 69 #include "logic.h"70 #include "temps.h"71 #include "ener.h"72 74 #include "description.h" 73 #include "serre.h"74 75 !#include "com_io_dyn.h" 75 76 #include "iniprint.h" … … 1455 1456 call Gather_Field(teta,ip1jmp1,llm,0) 1456 1457 call Gather_Field(pk,ip1jmp1,llm,0) 1457 1458 call Gather_Field(phi,ip1jmp1,llm,0) 1458 1459 do iq=1,nqtot 1459 1460 call Gather_Field(q(1,1,iq),ip1jmp1,llm,0) 1460 1461 enddo 1461 1462 call Gather_Field(masse,ip1jmp1,llm,0) 1462 1463 call Gather_Field(ps,ip1jmp1,1,0) 1463 1464 call Gather_Field(phis,ip1jmp1,1,0) 1464 1465 if (mpi_rank==0) then 1465 1466 CALL writedynav(itau,vcov, 1466 1467 & ucov,teta,pk,phi,q,masse,ps,phis) 1467 1468 endif 1468 1469 #endif 1469 1470 !$OMP END MASTER … … 1477 1478 IF( MOD(itau,iecri).EQ.0) THEN 1478 1479 ! Ehouarn: output only during LF or Backward Matsuno 1479 1480 if (leapf.or.(.not.leapf.and.(.not.forward))) then 1480 1481 c$OMP BARRIER 1481 1482 c$OMP MASTER … … 1515 1516 call Gather_Field(ucov,ip1jmp1,llm,0) 1516 1517 call Gather_Field(teta,ip1jmp1,llm,0) 1517 1518 call Gather_Field(phi,ip1jmp1,llm,0) 1518 1519 do iq=1,nqtot 1519 1520 call Gather_Field(q(1,1,iq),ip1jmp1,llm,0) 1520 1521 enddo 1521 1522 call Gather_Field(masse,ip1jmp1,llm,0) 1522 1523 call Gather_Field(ps,ip1jmp1,1,0) 1523 1524 call Gather_Field(phis,ip1jmp1,1,0) 1524 1525 if (mpi_rank==0) then 1525 1526 1526 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis) 1527 endif 1527 1528 ! CALL writehist_p(histid,histvid, itau,vcov, 1528 1529 ! & ucov,teta,phi,q,masse,ps,phis) … … 1650 1651 call Gather_Field(teta,ip1jmp1,llm,0) 1651 1652 call Gather_Field(pk,ip1jmp1,llm,0) 1652 1653 call Gather_Field(phi,ip1jmp1,llm,0) 1653 1654 do iq=1,nqtot 1654 1655 call Gather_Field(q(1,1,iq),ip1jmp1,llm,0) 1655 1656 enddo 1656 1657 call Gather_Field(masse,ip1jmp1,llm,0) 1657 1658 call Gather_Field(ps,ip1jmp1,1,0) 1658 1659 call Gather_Field(phis,ip1jmp1,1,0) 1659 1660 if (mpi_rank==0) then 1660 1661 CALL writedynav(itau,vcov, 1661 1662 & ucov,teta,pk,phi,q,masse,ps,phis) 1662 1663 endif 1663 1664 #endif 1664 1665 !$OMP END MASTER … … 1706 1707 call Gather_Field(ucov,ip1jmp1,llm,0) 1707 1708 call Gather_Field(teta,ip1jmp1,llm,0) 1708 1709 call Gather_Field(phi,ip1jmp1,llm,0) 1709 1710 do iq=1,nqtot 1710 1711 call Gather_Field(q(1,1,iq),ip1jmp1,llm,0) 1711 1712 enddo 1712 1713 call Gather_Field(masse,ip1jmp1,llm,0) 1713 1714 call Gather_Field(ps,ip1jmp1,1,0) 1714 1715 call Gather_Field(phis,ip1jmp1,1,0) 1715 1716 if (mpi_rank==0) then 1716 1717 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis) 1717 1718 endif 1718 1719 ! CALL writehist_p(histid, histvid, itau,vcov , 1719 1720 ! & ucov,teta,phi,q,masse,ps,phis) -
LMDZ5/branches/testing/libf/dyn3dpar/massbar_p.F
r1910 r2641 18 18 #include "dimensions.h" 19 19 #include "paramet.h" 20 #include "comconst.h"21 20 #include "comgeom.h" 22 21 c -
LMDZ5/branches/testing/libf/dyn3dpar/massbarxy_p.F
r1910 r2641 17 17 #include "dimensions.h" 18 18 #include "paramet.h" 19 #include "comconst.h"20 19 #include "comgeom.h" 21 20 c -
LMDZ5/branches/testing/libf/dyn3dpar/massdair_p.F
r1910 r2641 18 18 #include "dimensions.h" 19 19 #include "paramet.h" 20 #include "comconst.h"21 20 #include "comgeom.h" 22 21 c -
LMDZ5/branches/testing/libf/dyn3dpar/nxgrarot_p.F
r1910 r2641 22 22 #include "paramet.h" 23 23 #include "comdissipn.h" 24 #include "logic.h"25 24 c 26 25 INTEGER klevel -
LMDZ5/branches/testing/libf/dyn3dpar/qminimum_p.F
r1910 r2641 8 8 #include "dimensions.h" 9 9 #include "paramet.h" 10 #include "comvert.h"11 10 c 12 11 INTEGER nq -
LMDZ5/branches/testing/libf/dyn3dpar/sw_case_williamson91_6.F
r1910 r2641 26 26 c 27 27 c======================================================================= 28 USE comconst_mod, ONLY: cpp, omeg, rad 29 USE comvert_mod, ONLY: ap, bp, preff 30 28 31 IMPLICIT NONE 29 32 c----------------------------------------------------------------------- … … 33 36 #include "dimensions.h" 34 37 #include "paramet.h" 35 #include "comvert.h"36 #include "comconst.h"37 38 #include "comgeom.h" 38 39 #include "iniprint.h" -
LMDZ5/branches/testing/libf/dyn3dpar/top_bound_p.F
r1910 r2641 4 4 SUBROUTINE top_bound_p(vcov,ucov,teta,masse,dt) 5 5 USE parallel_lmdz 6 USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound, 7 & tau_top_bound 8 USE comvert_mod, ONLY: presnivs, preff, scaleheight 9 6 10 IMPLICIT NONE 7 11 c 8 12 #include "dimensions.h" 9 13 #include "paramet.h" 10 #include "comconst.h"11 #include "comvert.h"12 14 #include "comgeom2.h" 13 15 … … 40 42 ! NB: top_bound sponge is only called from leapfrog if ok_strato=.true. 41 43 42 ! sponge parameters: (loaded/set in conf_gcm.F ; stored in comconst .h)44 ! sponge parameters: (loaded/set in conf_gcm.F ; stored in comconst_mod) 43 45 ! iflag_top_bound=0 for no sponge 44 46 ! iflag_top_bound=1 for sponge over 4 topmost layers … … 69 71 REAL uzon(jjp1,llm),vzon(jjm,llm),tzon(jjp1,llm) 70 72 71 integer i 73 integer i 72 74 REAL,SAVE :: rdamp(llm) ! quenching coefficient 73 75 real,save :: lambda(llm) ! inverse or quenching time scale (Hz) -
LMDZ5/branches/testing/libf/dyn3dpar/tourpot_p.F
r1910 r2641 23 23 #include "paramet.h" 24 24 #include "comgeom.h" 25 #include "logic.h"26 25 27 26 REAL rot( ip1jm,llm ) -
LMDZ5/branches/testing/libf/dyn3dpar/vitvert_p.F
r1910 r2641 2 2 c 3 3 USE parallel_lmdz 4 USE comvert_mod, ONLY: bp 5 4 6 IMPLICIT NONE 5 7 … … 28 30 #include "dimensions.h" 29 31 #include "paramet.h" 30 #include "comvert.h"31 32 32 33 REAL w(ip1jmp1,llm),convm(ip1jmp1,llm) -
LMDZ5/branches/testing/libf/dyn3dpar/vlsplt_p.F
r1910 r2641 25 25 #include "dimensions.h" 26 26 #include "paramet.h" 27 #include "logic.h"28 #include "comvert.h"29 #include "comconst.h"30 27 31 28 c … … 208 205 #include "dimensions.h" 209 206 #include "paramet.h" 210 #include "logic.h"211 #include "comvert.h"212 #include "comconst.h"213 207 c 214 208 c … … 539 533 c -------------------------------------------------------------------- 540 534 USE parallel_lmdz 535 USE comconst_mod, ONLY: pi 541 536 IMPLICIT NONE 542 537 c 543 538 #include "dimensions.h" 544 539 #include "paramet.h" 545 #include "logic.h"546 #include "comvert.h"547 #include "comconst.h"548 540 #include "comgeom.h" 549 541 c … … 932 924 #include "dimensions.h" 933 925 #include "paramet.h" 934 #include "logic.h"935 #include "comvert.h"936 #include "comconst.h"937 926 c 938 927 c -
LMDZ5/branches/testing/libf/dyn3dpar/vlspltgen_p.F
r1910 r2641 27 27 USE VAMPIR 28 28 USE infotrac, ONLY : nqtot 29 USE comconst_mod, ONLY: cpp 29 30 IMPLICIT NONE 30 31 … … 32 33 #include "dimensions.h" 33 34 #include "paramet.h" 34 #include "logic.h"35 #include "comvert.h"36 #include "comconst.h"37 35 38 36 c … … 105 103 106 104 107 108 109 110 111 105 ijb=ij_begin-iip1 106 ije=ij_end+iip1 107 if (pole_nord) ijb=ij_begin 108 if (pole_sud) ije=ij_end 109 112 110 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 113 111 DO l = 1, llm 114 112 DO ij = ijb, ije 115 113 tempe(ij) = teta(ij,l) * pk(ij,l) /cpp … … 191 189 192 190 if(iadv(iq) == 0) then 193 194 195 196 191 192 cycle 193 194 else if (iadv(iq)==10) then 197 195 198 196 #ifdef _ADV_HALO 199 200 & 197 call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu, 198 & ij_begin,ij_begin+2*iip1-1) 201 199 call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu, 202 200 & ij_end-2*iip1+1,ij_end) 203 201 #else 204 205 & 202 call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu, 203 & ij_begin,ij_end) 206 204 #endif 207 205 … … 215 213 call VTe(VTHallo) 216 214 c$OMP END MASTER 217 215 else if (iadv(iq)==14) then 218 216 219 217 #ifdef _ADV_HALO … … 239 237 c$OMP END MASTER 240 238 else 241 242 239 240 stop 'vlspltgen_p : schema non parallelise' 243 241 244 242 endif … … 261 259 262 260 if(iadv(iq) == 0) then 263 264 265 266 261 262 cycle 263 264 else if (iadv(iq)==10) then 267 265 268 266 #ifdef _ADV_HALLO … … 270 268 & ij_begin+2*iip1,ij_end-2*iip1) 271 269 #endif 272 270 else if (iadv(iq)==14) then 273 271 #ifdef _ADV_HALLO 274 272 call vlxqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,qsat, … … 276 274 #endif 277 275 else 278 279 276 277 stop 'vlspltgen_p : schema non parallelise' 280 278 281 279 endif … … 301 299 302 300 if(iadv(iq) == 0) then 303 304 305 306 301 302 cycle 303 304 else if (iadv(iq)==10) then 307 305 308 306 call vly_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv) 309 307 310 308 else if (iadv(iq)==14) then 311 309 312 310 call vlyqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv,qsat) 313 311 314 312 else 315 316 313 314 stop 'vlspltgen_p : schema non parallelise' 317 315 318 316 endif … … 324 322 325 323 if(iadv(iq) == 0) then 326 327 328 329 324 325 cycle 326 327 else if (iadv(iq)==10 .or. iadv(iq)==14 ) then 330 328 331 329 c$OMP BARRIER … … 350 348 c$OMP MASTER 351 349 call VTe(VTHallo) 352 c$OMP END MASTER 353 c$OMP BARRIER 354 else 355 356 350 c$OMP END MASTER 351 c$OMP BARRIER 352 else 353 354 stop 'vlspltgen_p : schema non parallelise' 357 355 358 356 endif … … 369 367 c$OMP MASTER 370 368 call VTe(VTHallo) 371 c$OMP END MASTER 369 c$OMP END MASTER 372 370 373 371 c$OMP BARRIER … … 375 373 376 374 if(iadv(iq) == 0) then 377 378 379 380 375 376 cycle 377 378 else if (iadv(iq)==10 .or. iadv(iq)==14 ) then 381 379 c$OMP BARRIER 382 380 … … 388 386 c$OMP BARRIER 389 387 else 390 391 388 389 stop 'vlspltgen_p : schema non parallelise' 392 390 393 391 endif … … 414 412 415 413 if(iadv(iq) == 0) then 416 417 418 419 414 415 cycle 416 417 else if (iadv(iq)==10) then 420 418 421 419 call vly_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv) 422 420 423 421 else if (iadv(iq)==14) then 424 422 425 423 call vlyqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv,qsat) 426 424 427 425 else 428 429 426 427 stop 'vlspltgen_p : schema non parallelise' 430 428 431 429 endif … … 436 434 437 435 if(iadv(iq) == 0) then 438 439 440 441 436 437 cycle 438 439 else if (iadv(iq)==10) then 442 440 443 441 call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu, 444 442 & ij_begin,ij_end) 445 443 446 444 else if (iadv(iq)==14) then 447 445 448 446 call vlxqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,qsat, … … 450 448 451 449 else 452 450 453 451 stop 'vlspltgen_p : schema non parallelise' 454 452 … … 469 467 DO ij=ijb,ije 470 468 c print *,'zq-->',ij,l,iq,zq(ij,l,iq) 471 c 472 469 c print *,'q-->',ij,l,iq,q(ij,l,iq) 470 q(ij,l,iq)=zq(ij,l,iq) 473 471 ENDDO 474 472 ENDDO -
LMDZ5/branches/testing/libf/dyn3dpar/vlspltqs_p.F
r1910 r2641 25 25 USE mod_hallo 26 26 USE VAMPIR 27 USE comconst_mod, ONLY: cpp 27 28 IMPLICIT NONE 28 29 … … 30 31 #include "dimensions.h" 31 32 #include "paramet.h" 32 #include "logic.h"33 #include "comvert.h"34 #include "comconst.h"35 33 36 34 c … … 236 234 #include "dimensions.h" 237 235 #include "paramet.h" 238 #include "logic.h"239 #include "comvert.h"240 #include "comconst.h"241 236 c 242 237 c … … 584 579 c -------------------------------------------------------------------- 585 580 USE parallel_lmdz 581 USE comconst_mod, ONLY: pi 586 582 IMPLICIT NONE 587 583 c 588 584 #include "dimensions.h" 589 585 #include "paramet.h" 590 #include "logic.h"591 #include "comvert.h"592 #include "comconst.h"593 586 #include "comgeom.h" 594 587 c -
LMDZ5/branches/testing/libf/dyn3dpar/writedynav_p.F
r1910 r2641 12 12 USE misc_mod 13 13 USE infotrac 14 USE comconst_mod, ONLY: cpp 15 USE temps_mod, ONLY: itau_dyn 16 14 17 implicit none 15 18 … … 42 45 #include "dimensions.h" 43 46 #include "paramet.h" 44 #include "comconst.h"45 #include "comvert.h"46 47 #include "comgeom.h" 47 #include "temps.h"48 #include "ener.h"49 #include "logic.h"50 48 #include "description.h" 51 #include "serre.h"52 49 #include "iniprint.h" 53 50 -
LMDZ5/branches/testing/libf/dyn3dpar/writehist_p.F
r1910 r2641 12 12 USE misc_mod 13 13 USE infotrac 14 USE temps_mod, ONLY: itau_dyn 15 14 16 implicit none 15 17 … … 43 45 #include "dimensions.h" 44 46 #include "paramet.h" 45 #include "comconst.h"46 #include "comvert.h"47 47 #include "comgeom.h" 48 #include "temps.h"49 #include "ener.h"50 #include "logic.h"51 48 #include "description.h" 52 #include "serre.h"53 49 #include "iniprint.h" 54 50
Note: See TracChangeset
for help on using the changeset viewer.