Changeset 3605 for LMDZ6/branches/Ocean_skin/libf/dyn3dmem
- Timestamp:
- Nov 21, 2019, 4:43:45 PM (5 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 12 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
-
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/caladvtrac_mod.F90
r1907 r3605 44 44 CALL allocate_u(massem,llm,d) 45 45 CALL allocate_u(pbaruc,llm,d) 46 pbaruc(:,:)=0 46 47 CALL allocate_v(pbarvc,llm,d) 48 pbarvc(:,:)=0 47 49 CALL allocate_u(pbarug,llm,d) 48 50 CALL allocate_v(pbarvg,llm,d) -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/call_calfis_mod.F90
r2603 r3605 52 52 CALL allocate_u(p,llmp1,d) 53 53 CALL allocate_u(pks,d) 54 pks(:)=0 54 55 CALL allocate_u(pk,llm,d) 56 pk(:,:)=0 55 57 CALL allocate_u(pkf,llm,d) 56 58 CALL allocate_u(phi,llm,d) -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/call_dissip_mod.F90
r1987 r3605 31 31 32 32 CALL allocate_u(ucov,llm,d) 33 ucov(:,:)=0 33 34 CALL allocate_v(vcov,llm,d) 35 vcov(:,:)=0 34 36 CALL allocate_u(teta,llm,d) 35 37 CALL allocate_u(p,llmp1,d) -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/conf_gcm.F90
r2665 r3605 25 25 USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, & 26 26 alphax,alphay,taux,tauy 27 USE temps_mod, ONLY: calend 27 USE temps_mod, ONLY: calend, year_len 28 28 29 29 IMPLICIT NONE … … 144 144 !Config 145 145 calend = 'earth_360d' 146 ! initialize year_len for aquaplanets and 1D 146 147 CALL getin('calend', calend) 148 if (calend == 'earth_360d') then 149 year_len=360 150 else if (calend == 'earth_365d') then 151 year_len=365 152 else if (calend == 'earth_366d') then 153 year_len=366 154 else 155 year_len=1 156 endif 157 147 158 148 159 !Config Key = dayref -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/gcm.F90
r2622 r3605 233 233 'GCM: AVANT iniacademic AVANT AVANT AVANT AVANT' 234 234 if (.not.read_start) then 235 start_time=0. 236 annee_ref=anneeref 235 237 CALL iniacademic_loc(vcov,ucov,teta,q,masse,ps,phis,time_0) 236 238 endif … … 368 370 369 371 !----------------------------------------------------------------------- 372 ! Initialisation des I/O : 373 ! ------------------------ 374 375 376 if (nday>=0) then 377 day_end = day_ini + nday 378 else 379 day_end = day_ini - nday/day_step 380 endif 381 382 WRITE(lunout,300)day_ini,day_end 383 300 FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//) 384 385 #ifdef CPP_IOIPSL 386 call ju2ymds(jD_ref + day_ini - day_ref, an, mois, jour, heure) 387 write (lunout,301)jour, mois, an 388 call ju2ymds(jD_ref + day_end - day_ref, an, mois, jour, heure) 389 write (lunout,302)jour, mois, an 390 301 FORMAT('1'/,15x,'run du ', i2,'/',i2,'/',i4) 391 302 FORMAT('1'/,15x,' au ', i2,'/',i2,'/',i4) 392 #endif 393 394 !----------------------------------------------------------------------- 370 395 ! Initialisation de la physique : 371 396 ! ------------------------------- … … 381 406 ENDIF ! of IF ((iflag_phys==1).or.(iflag_phys>=100)) 382 407 383 384 !-----------------------------------------------------------------------385 ! Initialisation des I/O :386 ! ------------------------387 388 389 if (nday>=0) then390 day_end = day_ini + nday391 else392 day_end = day_ini - nday/day_step393 endif394 395 WRITE(lunout,300)day_ini,day_end396 300 FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//)397 398 #ifdef CPP_IOIPSL399 call ju2ymds(jD_ref + day_ini - day_ref, an, mois, jour, heure)400 write (lunout,301)jour, mois, an401 call ju2ymds(jD_ref + day_end - day_ref, an, mois, jour, heure)402 write (lunout,302)jour, mois, an403 301 FORMAT('1'/,15x,'run du ', i2,'/',i2,'/',i4)404 302 FORMAT('1'/,15x,' au ', i2,'/',i2,'/',i4)405 #endif406 408 407 409 ! if (planet_type.eq."earth") then -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/guide_loc_mod.F90
r2740 r3605 1212 1212 enddo 1213 1213 endif 1214 if (pole_ nord) then1214 if (pole_sud) then 1215 1215 do i=1,iip1 1216 1216 qgui1(ip1jm+i,l)=qgui1(ip1jm+1,l) -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/iniacademic_loc.F90
r2622 r3605 101 101 time_0=0. 102 102 day_ref=1 103 annee_ref=0103 ! annee_ref=0 104 104 105 105 im = iim -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/integrd_mod.F90
r1907 r3605 23 23 CALL allocate_u(deltap,llm,d) 24 24 CALL allocate_u(ps,d) 25 ps(:)=0 25 26 26 27 -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/leapfrog_loc.F
r2622 r3605 29 29 USE call_dissip_mod, ONLY : call_dissip 30 30 USE call_calfis_mod, ONLY : call_calfis 31 USE leapfrog_mod 31 USE leapfrog_mod, ONLY : ucov,vcov,teta,ps,masse,phis,q,dq 32 & ,ucovm1,vcovm1,tetam1,massem1,psm1,p,pks,pk,pkf,flxw 33 & ,pbaru,pbarv,du,dv,dteta,phi,dp,w 34 & ,leapfrog_allocate,leapfrog_switch_caldyn,leapfrog_switch_dissip 35 32 36 use exner_hyb_loc_m, only: exner_hyb_loc 33 37 use exner_milieu_loc_m, only: exner_milieu_loc … … 1535 1539 1536 1540 #ifdef INCA 1537 call finalize_inca 1541 if (type_trac == 'inca') then 1542 call finalize_inca 1543 endif 1538 1544 #endif 1539 1545 … … 1583 1589 1584 1590 #ifdef INCA 1585 call finalize_inca 1591 if (type_trac == 'inca') then 1592 call finalize_inca 1593 endif 1586 1594 #endif 1587 1595 … … 1732 1740 1733 1741 #ifdef INCA 1734 call finalize_inca 1742 if (type_trac == 'inca') then 1743 call finalize_inca 1744 endif 1735 1745 #endif 1736 1746 … … 1820 1830 1821 1831 #ifdef INCA 1822 call finalize_inca 1832 if (type_trac == 'inca') then 1833 call finalize_inca 1834 endif 1823 1835 #endif 1824 1836 -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/temps_mod.F90
r2601 r3605 13 13 INTEGER annee_ref 14 14 INTEGER day_ref 15 INTEGER year_len 15 16 REAL dt ! (dynamics) time step (changes if doing Matsuno or LF step) 16 17 REAL jD_ref ! reference julian day date (beginning of experiment) -
LMDZ6/branches/Ocean_skin/libf/dyn3dmem/vlsplt_loc.F
r2765 r3605 19 19 include "dimensions.h" 20 20 include "paramet.h" 21 include "iniprint.h" 21 22 c 22 23 c … … 872 873 include "dimensions.h" 873 874 include "paramet.h" 875 include "iniprint.h" 874 876 c 875 877 c … … 1027 1029 ELSE ! countcfl>=1 1028 1030 1029 PRINT*,'vlz passage dans le non local' 1031 IF (prt_level>9) THEN 1032 WRITE(lunout,*)'vlz passage dans le non local' 1033 ENDIF 1030 1034 c --------------------------------------------------------------- 1031 1035 c Debut du traitement du cas ou on viole le CFL : w > masse … … 1059 1063 c le critère 1060 1064 DO WHILE (countcfl>=1) 1061 print*,'On viole le CFL Vertical sur ',countcfl,' pts' 1065 IF (prt_level>9) THEN 1066 WRITE(lunout,*)'On viole le CFL Vertical sur ',countcfl,' pts' 1067 ENDIF 1062 1068 countcfl=0 1063 1069
Note: See TracChangeset
for help on using the changeset viewer.