Changeset 469
- Timestamp:
- 01/26/24 17:34:15 (4 months ago)
- Location:
- branches/GRISLIv3/SOURCES
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/GRISLIv3/SOURCES/3D-physique-gen_mod.f90
r468 r469 243 243 ! ===================== File id ========================================= 244 244 integer :: num_param = 10 !< Id of param_list.dat 245 integer :: num_coupe = 11 !< Id of coupe.ijz246 integer :: num_vitbil = 12 !< Id of vitbil247 integer :: num_newton = 15 !< Id of corrextion-newton.res248 245 integer :: num_file1 = 18 !< Id of 249 246 integer :: num_file2 = 19 !< Id of 250 247 integer :: num_forc = 20 !< Id of forcage 251 integer :: num_dat = 21 !< Id of dat252 248 integer :: num_cptr = 30 !< Id of cptr 253 249 integer :: num_debug = 33 !< Id of ritz debug 254 250 integer :: num_ritz = 35 !< Id of ritz 255 integer :: num_templlib = 36 !< Id of ritz templlib256 251 integer :: num_rep_42 = 42 !< Id of reponse 257 252 integer :: num_trace = 81 !< Id of trace-reboot 258 253 integer :: num_file = 88 !< Id of 259 integer :: num_ts_ritz = 89260 254 integer :: num_kelvin = 90 !< Id of kelvin.res 261 integer :: num_phil = 91262 integer :: num_templin = 98 !< Id of templin.res263 integer :: num_sealevel = 110264 integer :: num_mk_vitbil = 122 !< Id of masque_vitbil265 integer :: num_ic_vo = 171 !< vostok266 integer :: num_ic_by = 172 !< byrd267 integer :: num_ic_dm = 173 !< dml268 integer :: num_ic_dc = 174 !< dome C269 integer :: num_ic_df = 175 !< dome Fuji270 integer :: num_file3 = 993271 integer :: num_file4 = 995 !< Id of272 integer :: num_coor = 2004 !< Id of coord-Ant-40km273 255 274 256 ! Variables communes au main et aux subroutines du step -
branches/GRISLIv3/SOURCES/lineartemp-0.2.f90
r467 r469 40 40 ! T(I,J,K)=T(I,J,1)-(T(I,J,1)-TPMP(I,J,NZ))*(K-1)/(NZ-1) 41 41 T(I,J,K)=T(I,J,1)-(T(I,J,1)-T(I,J,NZ))*(K-1)/(NZ-1) 42 ! write(num_templin,*) T(I,J,K)43 42 end do 44 43 do K=NZ+1,NZ+NZM -
branches/GRISLIv3/SOURCES/main3D-0.4-40km.f90
r468 r469 131 131 subroutine grisli_init 132 132 133 USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen,num_sealevel, & 134 num_ts_ritz,num_ic_vo,num_ic_by,num_ic_dm,num_ic_dc,num_ic_df, & 133 USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen, & 135 134 s,h,b,bsoc,flot,mk,mk0,uxbar,uybar,hwater,time,timemax,ndebug,ndebug_max 136 135 use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,runname,itracebug … … 183 182 DTCPT=dtout 184 183 185 186 ! ************ OPEN FILES.RITZ ****************187 188 if ((geoplace.eq.'anteis1').or.(geoplace.eq.'ant20km')) then189 ! fichier de reference pour le niveau des mers190 open(num_sealevel,file=TRIM(DIRNAMEOUT)//'sealevel'//runname//'.ritz',position="append")191 open(num_ts_ritz,file=TRIM(DIRNAMEOUT)//'ts_'//runname//'.ritz',position="append")192 open(num_ic_vo,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'vo.ritz',position="append")193 open(num_ic_by,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'by.ritz',position="append")194 open(num_ic_dm,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'dm.ritz',position="append")195 open(num_ic_dc,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'dc.ritz',position="append")196 open(num_ic_df,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'df.ritz',position="append")197 endif198 199 184 !------------------------------ INITIALISATION ---------------------------- 200 185 ! -
branches/GRISLIv3/SOURCES/tab-litho-0.3.f90
r446 r469 88 88 WE(:,:) = WE(:,:)/(som*ROM*g) 89 89 90 91 92 ! ATTENTION DESSOUS VALABLE SEULEMENT EN 40 km93 94 ! la partie ci-dessous doit etre mise dans une routine specifique Philippe95 ! pour faire exactement la meme isostasie que Philippe96 !if ((icouple.eq.2).or.(icouple.eq.4)) then97 98 ! LBLOC=1099 100 ! do I=-LBLOC,LBLOC101 ! do J=-LBLOC,LBLOC102 ! We(i,j)=0. ! mise a 0 partout103 ! end do104 ! end do105 ! open(num_phil,file='../INPUT/iso-philippe')106 107 ! read(num_phil,*)108 ! read(num_phil,*)109 ! read(num_phil,*)110 ! read(num_phil,*)111 ! do J=-10,10112 ! read(num_phil,*)(We(i,j),i=-10,10)113 ! end do114 ! close(num_phil)115 ! do i=-LBLOC,LBLOC116 ! do j=-LBLOC,LBLOC117 ! we(i,j)=we(i,j)/11030/g/rom118 ! end do119 ! end do120 !endif121 122 90 return 123 91 end subroutine tab_litho
Note: See TracChangeset
for help on using the changeset viewer.